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The  problem  of  scattering  of  electromagnetic  waves  from  a 
perfectly  electrically  conducting  (PEC)  half -plane  is  a  classical 
problem  in  the  field  of  electromagnetics  and  is  the  basis  of  the 
GTD  (geometrical  theory  of  diffraction) .  The  application  of  this 
theory  to  modern  aerospace  structures  is  hampered  due  to  the 
presence  of  media  discontinuities  which  alter  the  scattering 
behavior.  Therefore,  the  study  of  diffraction  in  the  presence  of 
such  discontinuities  Is  an  important  prerequisite  to  applying  the 
GTD  to  practical  problems.  Two  related  problems  are  considered  in 
this  study  which  involve  planar  media  discontinuities.  The  first 
configuration  concerns  a  planar  Interface  between  two  semi-infinite, 
homogerleous  media  with  a  PEC  half-plane  covering  half  of  the 
interface.  The  second  configuration  concerns  a  dielectric  slab 
embedded  in  a  contrasting  medium.  On  the  upper  surface  of  the  slab 
resides  a  PEC  half-plane.  Both  structures  are  assumed  to  be  excited 
by  a  transverse  magnetic  plane  wave.  , 

CThe  boundary  value  problems  are  structured  as  integral  equations 
involving  the  electric  field  which  is  tangential  to  the  interface,  s 
The  integral  equations  are  Fourier  transformed  to  create  Wiener-Hopr 
equations  which  are  then  solved.  The  process  of  factorization  is 
performed  using  a  formal  Integral  factorization  which  is  converted 
into  a  form  that  is  amenable  to  numerical  evaluation.  The  solution 
is  structured  as  a  radiation  integral  and  is  evaluated  asymptotically 
to  arrive  at  the  far  field  ray-optic  contributors. 

Data  is  presented  for  the  interface  problem  to  exhibit 
distinctive  properties  of  the  solution.  Sample  cases  of  practical 
interest  are  given.  The  data  is  considered  in  the  light  of  present 
theory  concerning  wave  structure  about  a  dielectric  interface  and 
an  analysis  of  the  significant  features  of  the  data  is  conducted. 
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CHAPTER  I 


INTRODUCTION 

The  objective  of  this  study  is  the  solution  of  two  related  prob¬ 
lems  concerning  the  diffraction  of  time-harmonic,  electromagnetic,  plane 
waves.  The  first  problem  deals  with  the  geometry  given  in  Figure  1.1  and 
is  henceforth  referred  to  as  the  interface  problem.  As  seen  in  the  fig¬ 
ure,  two  half-spaces  are  separated  by  a  planar  interface  which  is  coin¬ 
cident  with  the  x'-z'  plane.  On  the  interface  is  a  perfectly  electrically 
conducting  (PEC),  infinitely  thin  half-plane.  The  half-plane  lies  in  the 
plane  y'  *  0  and  extends  along  the  positive  x'  direction  from  an  edge 
which  is  coincident  with  the  z'  axis.  An  incident  plane  wave  propagating 
strictly  in  the  x'-y'  plane  in  medium  1  (characterized  by  constitutive 
parameters  (y^, ep)  impinges  on  the  half-plane.  There  is  assumed  to  be  no 
variation  of  fields  in  the  z'  direction,  hence  the  problem  is  two-dimen¬ 
sional. 

The  second  problem,  referred  to  as  the  slab  problem,  consists  of  the 
above  half-plane  residing  on  a  dielectric  slab  characterized  by  (y^.e^) 
which  is  in  turn  imbedded  in  another  medium  (yg,£p  as  shown  in  Fig¬ 
ure  1.2. 

Aside  from  the  academic  interest  of  studying  such  configurations, 
they  resemble  or  model  junctions  involving  composite  materials  such  as 
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Figure  1.2.  Geometry  of  the  slab  problem 
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occur  in  modern  aerospace  structures.  Therefore  diffraction  from  such 
planar  composite  structures  is  of  practical  interest.  The  first  prob¬ 
lem  is  a  fundamental  building  block  toward  solving  and  interpreting  more 
complicated  diffraction  problems  involving  media  interfaces.  The  second 
structure  serves  as  a  canonical  problem  modeling  a  lapped- junction  between 
a  metallic  and  a  composite  sheet.  It  is  desirable  to  be  able  to  ana¬ 
lyze  these  structures  using  the  geometrical  theory  of  diffraction  (GTD), 
but  no  diffraction  coefficients  are  available  for  either  of  them.  As  a 
result,  scattering  problems  involving  these  canonical  structures  cannot 
be  modeled  using  the  GTD  at  present.  One  purpose  of  this  study  is  to 
solve  these  problems  in  an  asymptotic  limit  and  extract  the  ray  consti¬ 
tuents  of  the  solution  in  order  to  derive  the  correct  diffraction  coef¬ 
ficients  for  GTD  modeling  of  composite/metal  junctions. 

The  scope  of  study  of  the  problem  geometries  is  as  follows.  For 
the  interface  problem  the  equations  are  formulated  for  four  separate 
cases:  a  TM  (transverse  magnetic)  polarized  incident  plane  wave  with  the 
integral  equations  formulated  first  over  the  PEC  screen  and  then  over  the 
remaining  (aperture)  portion  of  the  interface  plane  and  a  TE  (transverse 
electric)  polarized  incident  wave  again  formulated  over  each  half  of  the 
interface  plane.  The  case  of  the  TM  wave  formulated  over  the  aperture  is 
developed  in  full  detail  from  the  integral  equation  formulation  to  the 
asymptotic  analysis  of  the  resulting  radiation  integrals.  The  other  cases 
can  be  developed  in  a  like  manner  and  the  similarity  of  the  analysis  is 
indicated.  The  interface  problem  is  handled  in  full  generality  with  the 


constitutive  parameters  of  the  media  allowed  to  be  independent  and  lossy 
(though  constrained  to  be  passive).  In  the  slab  problem,  only  the  TM  in¬ 
cident  wave  is  considered  with  the  integral  equation  formulated  over  the 
PEC  screen.  The  media  are  assumed  to  be  lossless  with  the  permeabilities 
equal  and  the  permittivity  of  the  slab  greater  than  that  of  the  surround¬ 
ing  medium. 

The  solution  of  diffraction  from  a  half-plane  can  be  traced  back  to 
the  problem  of  scattering  from  a  half-plane  in  a  homogeneous  medium  which 
was  solved  by  Sommerfeld  (1896).  Over  the  years  this  problem  has  been 
treated  from  a  number  of  points  of  view.  The  application  of  the  Wiener- 
Hopf  technique  to  the  half-plane  problem  is  described  in  Noble  (1958). 

The  Wiener-Hopf  approach  was  developed  to  deal  with  a  class  of  integral 
equations  on  a  semi-infinite  domain  and  is  a  special  case  of  the  Riemann- 
Hilbert  problem  (Carrier,  Krook  and  Pearson  (1966)).  It  was  applied  to 
diffraction  problems  by  J.  Schwinger  and  E.  T.  Copson.  Application  of 
the  Wiener-Hopf  technique  has  been  the  principal  subject  of  several  books 
(Noble  (1958),  Mittra  and  Lee  (1971),  Weinstein  (1969)).  In  the  Sommer¬ 
feld  half-plane  problem,  the  diffracted  fields  are  couched  in  the  form  of 
a  Fourier  inversion  integral  which  can  be  computed  as  Fresnel  integrals 
or  approximated  asymptotical  ly  as  rays.  This  ray-optic  result  was  uti¬ 
lized  by  Keller  (1957,1962)  to  formulate  a  high  frequency  electromagnetic 
scattering  theory  which  he  called  the  geometrical  theory  of  diffraction 
(GTD).  The  GTD  has  been  extensively  developed  by  others  (see  Hansen 
(1981)  for  a  survey).  The  solutions  of  various  canonical  problems  are 
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collected  in  the  GTD  and  one  subsequently  models  scattering  from  complex 
structures  by  isolating  effects  (localization)  and  synthesizing  an  approx¬ 
imate  response  from  the  canonical  building  blocks.  Each  new  solution  ex¬ 
tends  the  capability  of  the  theory,  and  the  solutions  presented  here  are 
intended  to  add  two  new  canonical  structures  to  those  available  for  such 
ray  solutions. 

The  Uiener-Hopf  technique  is  used  in  this  research  since  it  is  appli¬ 
cable  to  a  planar  geometry  which  leads  to  an  integral  equation  on  a  semi¬ 
infinite  domain,  and  it  yields  exact  results.  The  difficulty  in  Wiener- 
Uopf  analysis  arises  in  performing  the  required  factorization  and  decom¬ 
position  steps.  Therefore,  any  new  factorization  is  of  interest  as  it 
increases  the  collection  of  known  factorizations.  The  interface  problem 
has  been  addressed  by  several  authors  and  the  required  factorization  was 
not  found  to  be  available.  A  reduced  form  of  the  interface  problem  invol¬ 
ving  wave  propagation  across  a  seashore  was  solved  by  Clenmiow  (1966), 

Bazar  and  Karp  (1962),  Heins  and  Feshbach  (1954),  and  others.  Recogniz¬ 
ing  that  a  formal  factorization  integral  exists  in  certain  cases,  several 
authors  have  performed  an  analysis  of  problems  closely  related  to  the  in¬ 
terface  problem.  An  analysis  by  Heitman  and  van  den  Berg  (1975)  resulted 
in  a  numerical  evaluation  of  the  factorization  integral  and  did  not  attempt 
to  examine  the  asymptotic  field  behavior.  A  recent  analysis  by  Sunahara 
and  Sekiguchi  (1981)  considered  a  closely  related  problem.  However,  their 
solution  breaks  down  as  the  geometry  limits  into  the  interface  problem. 

Part  of  the  purpose  of  this  study  is  to  develop  an  efficient  integral 
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factorization  which  yields  all  the  analytical  properties  of  the  factor¬ 
ization  and  is  easy  to  evaluate  numerically.  This  represents  a  natural 
extention  of  the  application  of  the  Wiener-Hopf  procedure  from  problems 
where  Che  analytical  steps  required  can  be  performed  and  the  results  given 
in  terms  of  known  functions  to  problems  which  do  not  yield  such  closed 
form  or  determinate  solutions. 


CHAPTER  II 


ANALYTICAL  PRELIMINARIES 


2.1  Notation 


The  purpose  of  this  chapter  is  to  summarize  the  fundamental  mathe¬ 
matical  tools  required  for  the  solution  of  the  problems  of  interest.  It 
is  intended  as  a  reference  for  the  remainder  of  this  work.  Notation  and 
conventions  are  established  in  this  chapter  for  the  rest  of  the  volume. 
All  configurations  are  assumed  to  be  strictly  two-dimensional  (no  varia¬ 
tion  along  the  z  direction).  The  fundamental  time  convention  used  is 

jut  , 

e  ,  which  is  suppressed.  As  a  result.  Maxwell  s  equations  take  the 

form: 

^  •  e  *  q/e  (2. 1 . l.a) 

^  x  e  ■  -m  -  jwuh  (2.1.1.b) 

(yh)  ■  0  (2.1.1.c) 


±  ■*  -*• 
7  x  h  -  j  +  juee 


(2. 1 . 1 .d) 


where  e  is  the  electric  field,  h  is  the  magnetic  field,  q  is  the  electric 
charge  density,  j  is  the  electric  current  density  and  m  is  the  magnetic 
current  density. 

The  media  to  be  considered  are  passive  and,  in  general,  lossy.  They 
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can  be  characterized  by  complex  constitutive  parameters  ()j,e)  which  are 

complex  numbers  in  the  fourth  quadrant  (positive  real  part  and  nonpositive 

imaginary  part).  Also,  the  complex  wave  numbers  are  written  as: 

k  ■  k  -ik.  ,  where  k  >  0  and  k.>  0. 
r  i  r  i— 

Wiener-Hopf  analysis,  which  is  used  in  the  solution  of  the  problems 
to  be  considered,  is  performed  in  the  generalized  Fourier  transform  (spec¬ 
tral)  domain.  For  this  analysis  we  choose  the  following  Fourier  transform 
pair:* 


F(k  ) 

x 


F{f (x)} 


1 

/2tt 


E<V 


dx„ 


(2. 1.2. a) 


f(x) 


F"1{F(kx)> 


1 

flit 


so+jc  _1k  x 

F(kx)  e  X  dkx 


(2.1.2.b) 


k  °°+1  c 

X 


Here  the  transform  variable  is  allowed  to  be  complex.  The  constant  c 
is  suitably  chosen  for  convergence  of  the  Fourier  inversion  integral  in 
analogy  to  the  Bromwich  contour  of  the  Laplace  transform  (Mittra  and  Lee 
(1971a)). 

Various  methods  are  available  for  arriving  at  the  Wiener-Hopf  equa¬ 
tion  for  a  problem.  Because  of  the  pervasiveness  of  the  integral  equa¬ 
tion  formulation  for  electromagnetic  problems,  we  choose  to  use  this 
type  of  equation  as  the  starting  point  in  formulating  the  Wiener-Hopf 

^Throughout  this  work  we  use  capitalized  function  names  to  denote 
the  spectral  (Fourier  transform)  domain  image  of  a  spatial  domain  function. 
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equation.  The  integral  involved  will  be  a  convolution  of  a  desired 
unknown  in  the  problem  with  an  appropriate  Green's  function,  which,  for 
two-dimensional  problems,  is  of  the  form  H^(k  |x-Xq|).  When  the  inte¬ 
gral  equation  is  Fourier  transformed  with  respect  to  x,  according  to  the 
convolution  theorem,  the  integral  is  reduced  to  the  product  of  the  trans¬ 
form  of  the  unknown  times  the  transform  of  the  Green's  function.  For 
the  Green's  function  in  free  space,  we  have  (Jones  (1964)): 


F{g(p) }  ■  f|-  \  H<2)(k|p|)} 


(2.1.3) 


where  8  ■  nt^-k^  (the  details  of  the  definition  of  8  are  given 

below)  and  p  *  (x)x  +  (y)y. 

2.2  Specification  of  the  Value  of  8  *  ^k2  -k£ 

The  two  dimensional  Green's  function  is  the  solution  to  the  inhomo- 

-► 

geneous  wave  equation  with  a  line  source  at  P  *  0  where  we  require  waves 
emanating  from  the  line  source  to  be  outgoing  and  decaying.  This  means 
that  all  the  components  in  the  spectral  domain  along  the  Fourier  inver¬ 
sion  path  must  also  have  these  characteristics.  Observe  that  in  the  in- 
version  integral  of  (2.1.3)  for  g(  P  ),  we  must  have  Im(8)  £  0  to  give 
rise  to  decay  along  the  y  axis. 

Clearly  all  physical  waves  will  have  Re(S)  >  0  to  provide  outgoing 
wave  behavior  from  y  ■  0  and  will  satisfy  Im(8)  <0  so  that  the  outgoing 
waves  decay  (consistent  with  the  passiveness  of  the  medium).  Hence  any 
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physical  wave  will  lie  in  the  region  characterized  by  Im(6)£0,  so  the 
boundary  of  the  allowed  region  is  Im(B)*0,  that  is,  8  is  purely  real. 

In  the  0  plane  this  is  the  real  axis  and  in  the  82  plane  this  is  the 
positive  real  axis  so 

B2  -  [  (k2  - k2 )  -  (a2  -  T2)]  -  2j[krk±  +  aT]  (2. 2.1. a) 

where 

k  -  k  -  jk.  and  k  -  o  +  jt  . 

t  X  X 

So  in  the  82  plane,  the  boundary  is  given  by 

Im82  ■  -2(k  k  +ot)  -  0 

and 

Re  S2  -  t(k2 -k2)  -  (a2 -t2)]  >  0  , 

Figure  2.1  which  depicts  the  lines  specified  above  in  the 
kx  "  a+jT  Plane*  We  observe  that  the  condition  Im(B)  *  0  implies 
Im(B2)  *  0  and  Re(82)  0.  The  equality  obviously  gives  a  pair  of  hyper¬ 

bolas  in  the  second  and  fourth  quadrants.  The  branch  cuts  for  the  allow¬ 
able  region:  of  k^  (top  sheet  of  8  )  will  run  from  i  k  to  infinity  along 
these  hyperbolas.  To  decide  which  part  of  the  hyperbolas  contain  the 
branch  cuts  of  8,  we  observe  that  the  cuts  force  Im(8 )  <_  0  on  the  top  sheet 
which  implies  that  they  lie  along  Im(8)  ■  0  and  therefore  Re(B2)>^  0.  So 
the  branch  cuts  lie  along  the  lines  shown  in  Figure  2.1.  These  are  the 
conventional  hyperbolic  branch  cuts  and  preserve  decay  at  infinity.  Since 
our  analysis  will  take  us  off  the  real  axis  in  the  spectral  domain,  it  is 


(2.2. l.b) 

(2.2. l.c) 


aSfcw' » - 


45*  ASYMPTOTE 


BRANCH 


Kx=o-+jT 


Figure  2.1.  Domain  of  S  ■  /k2-k2 


important  to  understand  where  these  cuts  lie  and  to  recognize  that  moving 

them  may  cause  problems  in  closing  an  integration  contour  at  infinity. 

We  need  a  top  sheet  where  the  outgoing  wave  conditions  stated  above  are 

satisfied  on  the  Fourier  transform  path  (the  real  axis)  and  which  allows 

us  to  perform  any  semi-circular  closures  that  we  may  need  later.  The 

domain  in  Figure  2.2  shows  the  sign  of  Re(B)  for  points  in  the 

k  domain . 
x 

One  other  related  issue  must  be  dealt  with  before  we  leave  the  sub¬ 
ject  of  the  definition  of  8.  In  some  situations  subsequently,  we  eval¬ 
uate  the  factors  of  8  *  4c2-k2  *  A+k  *A-k  separately  and  must  impose 
the  above  conditions  in  so  doing.  The  complex  number  (k-kx)  can  be  re¬ 
presented  by  a  vector  pointing  from  +kx  to  k.  The  modulus  and  argument 
of  the  number  are  taken  as  the  length  of  the  vector  and  angle  with  res¬ 
pect  to  a  line  running  from  +kx  parallel  to  the  positive  real  axis  to¬ 
ward  infinity  as  shown  in  Figure  2.3.  The  complex  vector  representation 

for  (k+k  )  m  k  -(-k)  points  from  -k  to  k  and  we  utilize  a  similar  defi- 
XX  X 

nit  ion  for  its  argument  as  shown  in  Figure  2.3.  As  a  result 


8  -  A‘-k4  -  [(k  +  k  )  •  (k-k  )]** 

X  XX 


j 

e 


ei  +  02 


(2.2.2) 


is  found  to  obey  our  branch  definition  of  8  and  so  we  have  a  simple 


rule  for  establishing  the  arguments  of  A+k^  and  A-kT^  which  is  consis¬ 
tent  with  the  definition  of  8. 


Finally,  a  point  of  terminology  needs  to  be  established.  It  will 
prove  useful  to  consider  branch  definitions  of  8  which  have  branch  cuts 


Figure  2.2.  Upper  sheet  values  of  8. 
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that  are  different  than  the  hyperbolic  branch  cuts  used  in  this  section. 
These  alternative  branches  of  8  are  arrived  at  by  analytic  continuation 
of  the  function  from  the  original  branch  definition.  These  new  branch 
definitions  are  described  as  analytic  continuations  of  the  function,  8, 
arising  from  the  movement  of  the  hyperbolic  branch  cuts  of  6.  This 
language  is  used  to  preserve  the  intuitive  sense  of  connection  between 
successive  branch  definitions  and  an  apology  is  made  here  to  the  precise 
mathematician  for  the  use  of  this  phraseology. 

2.3  Edge  Condition 

In  order  to  uniquely  specify  the  field  solution  for  diffraction 
from  a  surface  which  contains  a  discontinuity  of  the  surface  normal 
vector,  it  is  necessary  to  specify  field  behavior  in  the  vicinity  of  the 
edge  since  unbounded  fields  arise.  This  problem  was  dealt  with  by  J. 
Meixner  (1954).  He  required  that  a  satisfactory  representation  for  the 
fields  in  a  small  volume  including  the  edge  must  be  limited  to  a  finite 
amount  of  energy.  For  a  problem  involving  scattering  from  a  wedge  in 
the  presence  of  different  media,  Mittra  and  Lee  (1971b)  have  given  the 
general  solution  for  three  media.  The  asymptotic  behavior  of  the  fields 
as  the  edge  of  the  half-plane  is  approached  from  large  distances  is  given 
below  for  the  problems  to  be  considered  in  this  study: 

E  transverse  H  transverse  Ex, Hz 

0(p_H)  0(p"S  0(pis)  .  (2.3.1) 
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2.4  Wiener-Hopf  Technique 

The  Wiener-Hopf  technique  applies  to  problems  formulated  on  a  semi¬ 
infinite  domain  in  a  separable  coordinate  system  and  its  use  is  ubiqui¬ 
tous  throughout  mathematical  physics  (see  Noble  (1958)).  For  example, 
in  a  two-dimensional  boundary  value  problem  where  the  boundary  values 
are  specified  along  a  line  so  that  along  one  semi-infinite  ray  a  Dirichlet 
condition  holds  and  along  the  remainder  of  the  line  a  Neumann  condition 
holds,  a  simple  solution  involving  the  Fourier  transform  of  the  differ¬ 
ential  equation  and  the  boundary  conditions  cannot  be  effected  since  the 

boundary  conditions  do  not  hold  along  the  entire  domain  of  the  trans¬ 
form.  We  do  not  repeat  the  development  leading  to  the  Wiener-Hopf  sol¬ 
ution  since  that  is  available  in  numerous  references  (see  Mittra  and  Lee 
(1971a)).  We  simply  state  the  results  of  the  development.  For  the 
convenience  of  the  reader,  we  follow  the  notation  of  Mittra  and  Lee. 

An  integral  equation  amenable  to  the  Wiener-Hopf  procedure  is  of 
the  form 

0 

— i—  f(x')  g(x,x')dx'  »h(x),  x  £  (-*,0)  ,  (2.4.1) 

The  semi-infinite  domain  involved  in  this  integral  equation  precludes 
Fourier  transformation  of  the  equation  as  it  stands.  In  order  to  trans¬ 
form  (2.4.1),  it  is  necessary  to  extend  the  equation  by  introducing  the 
definitions 


$  (x)  -  < 


0,  x  >  0 


f (x) ,  x  <  0 


and 
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(2.4.2) 


a  (x)  - 


0,  x  >  0 


h(x) ,  x  <  0 


(2.4.3) 


With  these  definitions  the  integral  equation  may  be  rewritten  as 


1 

fin 


$  (x)  g(x-x')dx'  -  a_(x)  +  b+(x),  xe  (-“>,*)  ,  (2.4.4) 


where  b+(x)  is  an  unknown  function  which  has  support  for  x  in  (0,°°  ) 
and  which  must  be  introduced  in  order  that  the  range  of  applicability  of 
the  original  equation  (2.4.1)  may  be  formally  extended  to  a  doubly- 
infinite  domain.  Since  b+(x)  is  not  a  part  of  the  original  problem  pre¬ 
scription,  it  is  a  new  unknown  which  should  be  recoverable  in  the  solution 
process . 

The  form  (2.4.4)  is  amenable  to  Fourier  transformation.  With  an  ap¬ 
peal  to  the  convolution  theorem  of  Fourier  transforms,  the  equation  be¬ 
comes 

$_(kx)  -G(kx)  -  A_(kx)  +  B+(kx)  .  (2.4.5) 

The  constraints  on  the  imaginary  part  of  k  are  dictated  by  the  exis¬ 
tence  of  the  Fourier  transforms  of  the  functions  in  a  particular  problem. 

The  Wiener-Hopf  procedure  leads  to  a  solution  for  the  unknown 


i 

l_  _ _ 
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functions  <t>+(kx)  and  B  (k^)  which  are  constructed  as  follows: 


S  (k  ) 

4-(V  ■  GlkJ 


VV  -  -  VV  •  G+<V 


(2.4.6) 


(2.4.7) 


The  functions  G  (k  )  and  G_(k  )  are  defined  through  the  factorization 
T  X  *  X 


of  the  transform  of  the  kernel,  GCk^),  as  follows: 


where 


G(kx)  -  G+(kx)  •  G_(kx) 

G+(k  )  is  regular  and  non-zero  for  Im(kx)  >  T_, 


(2.4.8) 


G  (k  )  is  regular  and  non-zero  for  Im(k  )  <  T.. 

-  x  x  + 

The  pair  of  functions  S  (k  )  and  S_(k  )  are  defined  through  the  decom- 
position  of  a  function  into  two  functions  which  are  each  individually 
regular  over  respective  halves  of  the  complex  plane  as  given  below: 


A.(k  ) 

GTOTT  -  W  +  S-<V 


(2.4.9) 


where 


S+(k^)  is  regular  for  Im(k J  _>  T_ 


Sjk^)  is  regular  for  ImC k  ^  <_  T+  . 

Thus,  in  summary,  we  can  recover  the  original  unknown  function, 
f(x),  through  its  transform  as  computed  in  (2.4.6).  The  other  unknown, 
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which  is  introduced  in  (2.4.4),  is  also  recovered.  To  construct  f(x),  we 

oust  find  functions  which  are  regular  and  non-zero  over  respective  halves 

of  the  complex  k  plane  according  to  the  factorization  definition  (2.4.8) 
x 

and  the  decomposition  definition  (2.4.9).  It  is  on  these  two  steps  that 
attention  is  focused  in  the  Wiener-Hopf  solution  process. 


2.5  Normalization 


Throughout  the  body  of  this  study,  a  normalization  has  been  used 
to  bring  the  scale  of  the  variables  to  a  consistent  and  convenient  level. 
Unnormalized  variables  in  the  initial  problem  formulations  are  primed. 

So  refers  to  the  actual  complex  permeability  of  medium  1  (physical 
property)  and  k^  refers  to  the  wave  number  of  medium  2.  Hence 
k^  *  a)  /  •  The  normalization  used  is  reflected  in  the  following 

equations  where  the  unprimed  variable  is  the  normalized  variable  and 
(Uq, e^)  are  the  constitutive  parameters  for  free  space: 


W' 

*1 


y  <"  v^o >  •  ui 


w e;  ■  veo  • 


The  effect  of  this  normalization  is  to  reference  values  of  any  variable 
to  the  value  which  holds  for  free  space  at  the  frequency  of  interest. 

Because  the  transform  variable,  kx>  is  normalized  to  the  free  space 
wavenumber,  kg  *  w  /  Mg*  £g,  i-c  is  necessary  to  normalize  distance  to  the 
free  space  wavelength: 

x  *  kg'x'*  (2ir»  (distance  in  wavelengths)) 

where  x'  is  the  unnormalized  distance.  This  ensures  that  k*x  *  k'*x'  , 
where  the  unprimed  quantities  are  normalized  and  the  primed  ones  are  not. 


CHAPTER  III 


WIENER-HOPF  ANALYSIS  OF  THE  INTERFACE  PROBLEM 

3.1  Statement  of  the  Problem 

The  diffracting  geometry  to  be  considered  is  shown  in  Figure  3.1. 

The  half  spaces  for  y'  >  0  and  y'  <  0  contain  media  characterized  by 
(y'.ep  and  respectively  and  their  interface  coincides  with  the 

x'-z’  plane.  The  media  are  general  and  characterized  by  constitutive 
parameters  which  are  complex  numbers  having  real  part  greater  than  zero 
and  imaginary  part  less  than  or  equal  to  zero.  An  infinitely  thin,  per¬ 
fectly  electrically  conducting  (PEC)  screen  lies  on  the  interface  for 
x'  >  0.  A  plane  wave  is  incident  from  medium  1  at  an  angle  9  to  the  PEC 
screen  (0  <  0  <  tt).  The  incident  wave  may  be  polarized  either  transverse 
magnetic  (TM)  or  transverse  electric  (TE)  relative  to  the  z'  axis  and  is 
invariant  with  respect  to  z'. 

We  seek  to  determine  the  total  fields  everywhere  in  space  in  the 
presence  of  the  medium  discontinuity  and  the  diffracting  PEC  half-plane. 
The  more  practical  result  developed,  however,  is  the  field  at  great 
distances  from  the  edge  of  the  half-plane,  that  is,  that  which  results 
from  applying  asymptotic  evaluation  techniques  to  the  inverse  Fourier 
transform  integrals  of  the  Wiener-Hopf  solutions  obtained. 

This  problem  may  be  formulated  in  terms  of  an  unknown  quantity  which 
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resides  on  the  interface  and  which  is  non-zero  only  over  half  of  the  do¬ 
main,  that  is,  for  x'  in  (-  «, «).  Either  the  component  of  the  total 
electric  field  tangential  to  the  interface  plane  is  considered  unknown 
for  x'  in  (-  °°,  0)  or  the  surface  current  density  on  the  PEC  half-plane 
for  x'  in  (O,00)  is  considered  as  unknown.  In  either  case  the  quantity 
in  question  vanishes  on  the  complementary  semi-infinite  domain  by  the 
physics  of  the  problem. 

We  proceed  in  this  chapter  to  state  four  integral  equations  result¬ 
ing  from  the  two  incident  field  polarizations  viewed  both  from  the  aper¬ 
ture  and  the  PEC  half-plane  formulations.  The  particular  case  of  TM  in¬ 
cidence  formulated  over  the  aperture  is  treated  in  detail.  The  other 
three  integral  equations  admit  to  similar  treatment.  We  conclude  the 
chapter  by  relating  them  to  the  TM  aperture  formulation. 

3.2  Integral  Equation  Formulations 

We  begin  by  considering  the  formulation  of  the  TM  scattering  problem 
as  an  integral  equation  with  the  tangential  electric  field  over  the  aper¬ 
ture  plane  (x  '  <  0)  as  the  unknown.  The  problem  conf igurat ion  is  shown 
in  Figure  3.1.  The  analysis  follows  the  procedure  described  by  Butler 
and  Umashankar  (1976).  The  problem  is  viewed  as  the  composition  of  two 
subsidiary  problems,  one  equivalent  problem  for  the  half  space  above  the 
interface  plane  (y  1  >0)  and  another  below  the  interface  (y  '  <0).  In 
each  equivalent  problem  the  PEC  half-plane  is  extended  to  cover  the  en¬ 
tire  interface  plane  and  an  unknown  surface  magnetic  current  is  intro- 
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duced  to  provide  the  correct  continuity  across  the  shorted  aperture  of 
the  tangential  components  of  the  fields.  The  resulting  integro-dif fer- 
ential  equation  is  given  below: 


hx”C<*') 


;  [>■ 


>2  I 


-jcue'  1  + 


(k p* 


3x'* 


eZ(x0) 


g1(x' 


xo)dxo 


+  jU)£^ 


for  x'  e  (-  °°,0)  , 


0 

ez(x0)  *  g2(x' .xQ)dx0  (3.2.1) 
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jk'x'cosS 
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where  klnc  *  - : - sin0  is  the  component  of  the  incident 

x 

magnetic  field  which  is  tangential  to  the  interface 
plane  (Eg  is  the  magnitude  of  the  electric  field  of 
the  incident  plane  wave  at  the  edge  of  the  PEC  screen 
and  ) , 

ez  is  the  total  electric  field  over  the  aperture, 

(e'  ,k^  )  ,(6^  ,1«2  )  ate  the  permittivity  and  wavenumber  for  the 
^  medium  in  the  upper/lower  half-spaces  respectively, 

and  g^,g  2are  Ch«  Green's  functions  appropriate  to  the  respective 
problems : 

8(1)  ^x' ,xo^“  "a- Ho  -xol>  * 

The  TM  problem  can  also  be  formulated  in  terms  of  the  1  -directed 
current  induced  on  the  PEC  half-plane  as  an  electric  field  integral 


equation.  The  resulting  integral  equation  is 


for  x'  £  (0,®), 


jz(x0)  -a(x\x0)dx0 


(3.2.2) 
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where  e  is  the  total  electric  field  present  with  the  PEC  half-plane 
removed,  that  is,  the  total  scattered  field  in  the  presence 
of  only  the  discontinuity  of  the  media, 

j z  is  the  unknown  induced  electric  current  on  the  PEC  half-plane, 

a  is  the  Green's  function  which  serves  as  the  kernel. 


The  kernel  is  usually  expressed  in  the  spectral  domain  as  it  involves  a 
Sonnnerfeld  integral.  For  this  problem  the  Fourier  transform  of  the 
kernel  a(x')  is 


A(kx) 


(3.2.3) 


Similar  procedures  can  be  formulated  for  the  TE  case.  The  problem 
can  be  formulated  as  an  integro-differential  equation  either  over  the 
aperture  or  over  the  PEC  half-plane.  The  resulting  equations  are  shown 
below: 


Aperture  Formulation 
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PEC  Formulation 
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for  x1  e  (0,»),  where  a(x',Xg)  is  the  kernel  of  the  integral  equation  in 


the  TM  case.  The  Fourier  transform  of  this  kernel  is 


A(k;) 


-J  £i 

m  +  if 


(3.2.6) 


It  is  noted  that  equations  (3.2.1)  and  (3.2.5)  are  integro-dif fei — 
ential  equations.  As  a  result,  the  subsequent  Fourier  transformation  of 
them  requires  integration  by  parts  which  results  in  the  formation  of  a 
term  commonly  referred  to  as  the  bilinear  concomitant.  A  detailed  exam¬ 
ination  of  the  conditions  for  the  vanishing  of  this  term  results  in  a 
boundary  condition  on  the  applicability  of  the  Wiener-Hopf  equation  in 
the  transform  domain.  The  remaining  two  integral  equations  are  seen  to 
involve  Sommerfeld  integrals  as  part  of  the  Green's  function  in  the  ker¬ 
nel.  On  examining  the  regions  of  validity  for  the  Fourier  transform  of 
the  integral  equations  resulting  from  the  two  formulations  in  a  parti¬ 
cular  polarization,  it  is  found  that  the  regions  are  complementary  in 
the  spectral  domain;  that  is,  the  union  of  the  domains  for  the  two  formu¬ 
lations  is  always  a  strip  in  the  transform  domain  which  is  bounded  by 
the  branch  points  for  that  wavenumber  in  the  problem  which  is  closest  to 
the  real  axis.  As  a  result  if  the  strip  of  analyticity  vanishes  in  one 
formulation  due  to  the  angle  of  the  incident  plane  wave,  it  is  found  that 


the  other  formulation  will  have  a  non-vanishing  Wiener-Hopf  strip. 
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The  unknown  function  b(x' )  is,  o-f  course,  not  required  by  the  original 
integral  equation.  It  is  an  artifice  in  the  Wiener-Hopf  procedure  and 
is  ultimately  eliminated. 


The  extended  domain  integral  equation  (3.3.1)  can  be  Fourier  trans¬ 


formed  according  to  (2.1.2)  to  obtain 


A’ (k')  +  B 
x 


f  k .  ]2  r 

'  (k')  -  -jioe  '  Jl rf  1-  $  '(k ' )  •  G.  (k ' ) 

X  X  ^  ^  X  X  X 


_  IV 

-JU3£'  m  1  -  4>’(k’)  •  G- (k ' ) 


(3.3.5) 


x  2  x 


where  we  have  used  the  convolution  theorem  of  Fourier  transforms.  The 
transformed  kernels  G  contain  the  Fourier  transforms  of  Hankel  func- 
tions  and  are  given  in  (2.1.3)  with  y  *  0. 

At  this  point  we  introduce  the  following  substitutions  and  normaliza- 
t ions : 


A_(V  l  1  En  sin  9 

-TT T"  *  A-(k  k  )  «  -X - 0 - 

inK-n  w  X  ink_ 


0  0  r2ir  H.  (k, cos  9  +  k  ) 
11  x 


(3.3.6) 


vv 

_  i, —  *  B'(k  k  )  ■  unknown  function  , 
n0  0  0  x 


(3.3.7) 
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With  these  normalizations,  the  result  is  that  for  1 

(-minOcli>  k2i)  -  t_)  <  T  <  (x+  *  minO<2i>  cos6))  (3.3.10) 

where  k  M  a+it  and  k  *  k  -ik , 
x  r  i 

the  following  equation  holds  in  the  Fourier  transform  domain: 


A_V  +  VV 


(3.3.11) 


This  is  the  Wiener-Hopf  equation  which  we  must  solve  for  ,  the  Fourier 
transform  of  the  electric  field  in  the  aperture  on  the  interface  between 
the  two  media. 

Requisite  to  the  Wiener-Hopf  procedure  is  the  determination  of  the 

common  strip  of  analyticity  in  the  k^  plane  for  the  quantities  appearing 

in  (3.3.11).  One  must  consider  the  region  of  existence  of  the  transform 

domain  functions  appearing  in  this  equation  in  the  light  of  the  Wiener- 

2 

Hopf  factorization  defined  in  (2.4.8). 

Figure  3.2  depicts  the  complex  k^  plane  and  the  boundaries  delimit¬ 
ing  the  half  planes  of  analyticity  for  the  functions  involved  and  the 
directions  from  each  boundary  in  which  the  regular  regions  extend.  (The 

*This  restriction  on  the  region  of  validity  of  the  Wiener-Hopf 
equation  arises  from  the  analyticity  and  existence  of  the  various  cons¬ 
tituent  transforms  and  is  discussed  subsequently. 

2 

We  consistently  use  the  (+)  subscript  to  denote  a  function  which 
is  regular  and  non-zero  in  the  upper  half  plane  and  (-)  to  denote  similar 
behavior  in  the  lower  half  plane. 
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subscript  notation  introduced  in  equations  (3.3.6)  through  (3.3.9)  refers 
to  the  regions  of  regularity  for  the  functions  defined  there  in  antici¬ 
pation  of  the  Wiener-Hopf  analysis  which  follows.)  The  common  region  of 
analyticity  for  all  the  components  is  observed  to  be  bounded  by  (3.3.10). 
We  also  note  that  the  Wiener-Hopf  equation,  when  formulated  over  the 
aperture,  does  not  apply  when  T+  <  T  which  occurs  if  k^cosQ  <  -k  . 
Observe  that  this  can  only  occur  for  9  >  tt/2  when  k^  >  k^. 

Our  Wiener-Hopf  system  is  (3.3.11)  and  applies  in  the  domain  defined 
by  (3.3.10).  To  use  the  Wiener-Hopf  procedure,  it  is  necessary  to  form 
the  factors  G+  and  G  defined  in  (2.4.8). 


3.4  Factorization  of  G(k  ) 

x 


We  consider  in  this  section  the  factorization  of  the  kernel  in  the 
Wiener-Hopf  equation  (3.3.11) 


G<V 


(3.4.1) 


Before  proceeding  to  formulate  a  formal  factorization  integral  for  G(kx) 
some  observations  about  the  function  are  in  order.  G(kx)  is  analytic 
except  on  the  branch  cuts  of  8^  ■  A|-k£  and  82  ■  /k|-k£ ,  and  is 
bounded  in  any  finite  region  of  the  complex  k^  plane.  It  is  a  multiple- 
sheeted  function  with  four  Riemann  sheets.  It  is  also  an  even  function. 
There  is  a  total  of  four  zeros  on  all  sheets  and  they  occur  in  pairs  on 
a  sheet.  Hence  two  sheets  will  have  no  zeros  and  two  will  have  two  zeros 
each. 
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It  is  useful  to  partially  factorize  G(kx>  by  identifying  explicit 
factors  which  are  regular  over  the  requisite  half  planes 


U1  +  U2 


s(klt)  '  <W  ’  “W  UjU2  “l  Vi  +  V2  “1  +  u2  B! 


u,  8- 

-r—  (3.4.2) 


The  first  factor  G,  (k  )  is  factorized  by  inspection: 
lx 


Gl(kx) 


yl  +  y2 


i,  +  k 
1  x 


ul  +  p2 


/k..  -  k  • 
1  x 


(3.4.3) 


The  G2  factor  is  now  the  function  which  remains  to  be  factorized.  We 
note  that  it  is  bounded  except  at  kx  *  ±k^.  In  particular  ■*  1  a s 

|k  It  is  a  two  sheeted  function  since  the  branch  points  are  intro¬ 
duced  only  in  the  ratio  of  ( ^  *-s  aPParent  if  one  considers 

the  value  of  the  function  on  a  closed  path  which  encircles  two  of  the  four 
branch  points.  It  is  an  even  function  (as  is  G(kx>)  and  is  bounded  every¬ 
where  except  at  k^  ■  ±kx  .  Since  there  are  two  branches  of  G2 ( kx ) ,  there 
will  only  be  two  roots  which  will  occur  and  they  will  be  located  at  the 


following  positions: 


k  *  ± 
x 


1.2  2  .2  2 
k2yl  ~  klU2 
2  2 
yl  -  y2 


*  yl  *  y2 


(3.4.4) 


If  y^  ■  y2»  there  are  no  roots  of  G2(kx)  uniess  £•]_  *  £2>i.e.  unless  the 
media  are  identical  signifying  the  degeneration  of  the  two  media  problem 
to  the  classical  half-plane  problem.  It  is  readily  determined  on  which 
sheet  of  G2(kx)  tke  r00ts  occur  by  substituting  the  roots  of  (3.4.4)  into 
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equation  (3.4.2)  for  G2(kx).  ^ey  can  on^y  occur  on  one  sheet. 

The  two  different  choices  of  branch  cuts  shown  in  Figure  3.3  prove 
useful  in  defining  G^  i-n  subsequent  developments.  The  cuts  shown  in 
Figure  3.3a  are  sectors  of  hyperbolas  which  are  asymptotic  to  the  coor¬ 
dinate  axes  and  pass  through  the  branch  points.  A  pair  of  these  cuts  is 
associated  with  the  respective  branch  points  at  +  k^  and  +  and  their 
choice  ensures  that 

la  B1  <  0  ,  Im  B2  <  0  (3.4.5) 

for  all  points  in  the  k^  plane  as  discussed  in  section  2.2.  With  6^  and 
8 2  defined  by  the  hyperbolic  branch  cuts,  it  is  difficult  to  discuss  the 
occurrence  of  the  roots  on  the  top  sheet  of  G2(kx).  It  is  desirable  to 
consider  an  analytic  continuation  of  the  definition  of  the  top  sheet  of 
^2^kx^  wh*-ch  manifests  the  fourth  quadrant  branch  cut  as  a  curved  line 
connecting  the  branch  points  k^  with  k^  (and,  similarly,  another  line 
connecting  -k^  with  -k2).  The  branch  cuts  for  this  analytic  continual 
tion  of  G2(kx)  are  shown  in  Figure  3.3b  and  account  for  the  two-sheeted 
nature  of  G2 .  That  this  branch  cut  deformation  is  permissible  is  a 
result  of  the  radicals  in  G2  appearing  as  a  quotient  only.  The  precise 
definition  of  this  branch  cut  and  the  properties  of  the  new  top  sheet  of 
G2(kx)  are  discussed  below. 

In  order  to  examine  the  functional  behavior  of  ^(k^)  with  the  new 
branch  cut  configuration,  the  functional  mappings  implied  by  G2  ( kx )  can 
be  viewed  as  a  series  of  conformal  mappings.  For  simplicity,  only  values 


Definition  of  G, (k  )  with  hyperbolic 

t  X 


branch  cuts. 


of  k  with  Im(k  )  <  0  will  be  considered.  (Similar  arguments  apply  for 
Im(k^)  >  0  also.)  The  mappings  are  depicted  in  Figure  3.4  along  with  the 
critical  points  k^  *  k^,  ,  «. 

The  definition  of  the  new  branch  cut  of  G.(k  )  is  shown  in  the  v 

2  x 

2  2 

plane  of  Figure  3.4.  The  branch  cut  connecting  k^  and  k^  is  a  straight 

line  which,  of  course,  goes  through  k  a  <*>.  The  transform  shown  as  B(k  ) 

x  x 

is  a  Mobius  mapping  and  hence  circle-preserving.  However  B(<»)  *  1  and 
B(k^)  ■  0  which  causes  the  branch  cut  to  map  into  the  negative  real  axis 
in  the  B  plane.  In  order  to  define  D  ■  /§,  it  is  necessary  to  specify 
both  the  branch  cut  and  a  point  on  the  top  sheet  of  D.  For  this  reason 
the  point  D(°°)  ■  i  is  used.  This  causes  G^Ck^)  *  as  I  ^  I"*  00  an<*  Pre“ 
serves  the  asymptotic  behavior  of  required  by  the  integral  factor¬ 

ization  theorem  which  we  apply  later.  Since  B(°°)  “  D(<=°)  *  1  and  D  *  /B  , 
it  is  clear  that  the  specified  branch  cut  maps  into  the  imaginary  axis  in 
the  D  plane.  Observe  that  the  root  of  G2^lcx'  occurs  at  D  *  — (^2 ^  * 
Since  both  an<*  are  >  ^-n  general,  complex  numbers  in  the  fourth 
quadrant,  their  ratio  is  a  complex  number  restricted  to  the  first  or 
fourth  quadrants.  Therefore,  roust  lie  in  the  left  half  of  com¬ 

plex  D  plane;  that  is,  the  root  of  D  must  lie  on  the  lower  sheet  of  D. 

This  proves  that  the  branch  definition  of  G2  which  is  specified  by  the 

2 

straight  cut  between  k^  and  (and  hence  -k^  and  -k^  )  in  the  Y  *  k 

plane  and  the  asymptotic  value  of  lim  G2^x^  *  *  39  I  I"*30  has  no  r00ts 
on  the  top  sheet.  This  cut  in  the  Y  plane  may  be  mapped  back  into  the 
k  plane,  of  course.  However  it  proves  less  complicated  and  no  more 
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Figure  3.4, 
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costly  to  transform  the  integral  in  kx  to  an  integral  in  Y  and  thereby 
obviate  the  need  to  explicitly  define  the  cut  in  the  k^  plane. 

We  proceed  to  perform  a  formal  integral  factorization  of  G^Ck^) 
defined  in  (3.4.2).  The  following  factorization  theorem  (Mittra  and  Lee 
(1971a),  Noble  (1958))  is  applicable  to  G2  as  we  have  constructed  it: 


Let  (»2(kx)  be  regular  and  nonzero  in  the 
strip  x_  <  t<t+.  Within  this  strip,  ■+■  1 

uniformly  as  [ O' f  ■*  °°.  Then  G£(kx)  can  be  fac¬ 
torized  such  that 


G2(kx)-  G2(k![)-G2(kJt)  , 


is  regular  and  nonzero  in  the  upper  half 
kx~plane  defined  by  T  >  T_  ;  and 


is  regular  and  nonzero  in  the  lower  half 
kjj-plane  defined  by  T  <  T+. 

Let  us  consider  the  integral  in  (3.4.6).  Figure  3.5  shows  the  inte- 
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gration  path  in  the  5  plane.  Since  we  are  considering  a  case  which  has 
no  roots  of  G2  on  the  top  sheet,  no  additional  cuts  are  required  for  the 
logarithm.  This  analytic  continuation  of  G2  has  precisely  the  same  val¬ 
ues  in  the  complex  6  plane  as  the  original  definition  of  the  top  sheet 
of  G2  at  all  points  except  those  in  the  region  bounded  by  the  hyperbolic 
cuts  and  the  new  branch  cut  connecting  k^  and  k2  (and  the  symmetric  image 
of  this  region  in  the  second  quadrant).  Since  no  roots  occur  in  this  de¬ 
finition  of  the  top  sheet  of  G2  and  the  function  is  bounded  at  all  points 
in  the  <$  plane  except  <$  *  ±  k^  ,  it  is  clear  that  no  branch  cut  of  the  lo¬ 
garithm  specified  in  equation  (3.4.6)  need  appear  on  the  top  sheet  of  the 
integrand  shown  in  Figure  3.5. 

We  observe  that  the  principal  value  of  the  logarithm  is  an  acceptable 
definition  for  the  integrand  in  (3.4.6).  Such  a  definition  can  be  seen 
to  have  the  logarithmic  branch  cut  on  the  lower  sheet  of  the  integrand. 
Considering  the  plot  of  the  D  plane  in  Figure  3.4,  the  logarithmic  cut 
will  run  from  the  root  of  G2  along  a  line  connecting  D(°°)  *  1  and  the 
root.  It  will  run  to  D(k^)  *  00  in  the  direction  opposite  to  D(co)  =  1. 

We  may  close  the  integration  contour  at  infinity  either  up  or  down 

with  no  change  in  the  value  of  the  integral  since  the  integrand  is 

uniformly  asymptotic  to  zero  along  any  arc  at  infinity.  Closure  down  is 

shown  in  Figure  3.6  where  we  show  the  integration  path  taken  in  the 

negative  direction  (for  convenience  later).  The  region  enclosed  by  the 

contour  is  analytic  and  therefore  the  contour  integral  equals  zero. 

Also  the  contributions  from  C  and  C  cancel  while  all  contributions 

1  2 
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from  are  zero. 


Therefore  we  have 


*  - 

=. 

(CT)  1 

(-0  C 

(3.4.8) 


We  have  converted  the  original  infinite  line  integral  into  a  finite  loop 

integral  encircling  the  branch  cut  from  ■  k^  to  5  *  k^. 

We  now  consider  the  loop  integral  around  the  branch  cut.  The  loop 

integral  may  be  considered  as  four  separate  integrals  as  shown  in  Figure 

3.7.  We  have  taken  the  branch  cut  as  the  line  connecting  k^  and  k^ 

2 

specified  implicitly  by  the  straight  cut  in  the  Y  ■  k^  plane.  Therefore 
we  have  the  following: 


In  G2(6) 
5  -  k 

-°°+j  c  x 


d6 


(3.4.9) 


Using  standard  analysis  one  can  show  that  the  contour  contributions 
vanish  on  the  two  semi-circular  arcs  and  C^.  Thus  the  original  inte¬ 
gral  reduces  to 


TiC  G2(<5) 


5  -  k 


d  6 


S,n  G2(6) 


5  -  k 


d  5 


-00+ jc 


Cl  +  C3 


(3.4.10) 


For  purpose  of  computation  it  is  necessary  to  establish  the  value  of 
G£  for  points  on  the  integration  path  C  This  is  straightforward  since 
the  path  specified  by  C^does  not  cross  the  original  hyperbolic  cuts  in 
Figure  3.7  and  the  original  hyperbolic  branch  values  of  the  square  roots 
may  be  used  along  C  This  can  be  shown  by  observing  that  the  hyperbolas 
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correspond  to  straight  lines  in  the  y  *  kx  plane  which  are  parallel  to  the 
real  axis  and  pass  through  the  branch  points.  Hence  the  hyperbolic  lines 
through  k^  and  k 2  in  the  kx  plane  of  a  problem  will  map  into  straight 
parallel  lines  in  the  y  plane  while  path  Carnaps  into  a  straight  line 
which  connects  the  branch  points.  So  the  value  of  on  a  path  will 
always  be  that  associated  with  the  hyperbolic  branch  cut  definition. 


As  depicted  in  Figure  3.7,  the  integral  along  C  ^  (path  from  k^  to  k^) 
is  on  the  original  top  sheet  of  the  hyperbolic  definition  of  and  the 
integral  along  C  is  on  the  original  lower  sheet  of  G^  (in  the  limit  as 
C  approaches  the  branch  cut).  This  is  due  to  the  fact  that  in  the  prob¬ 
lem  shown  in  Figure  3.7,  the  point  6  *  k  lies  on  the  original  hyperbolic 
branch  cut  which  was  nearer  the  coordinate  axes  in  the  6  plane.  Had  k^ 
lain  on  the  inside  hyperbola,  we  would  have  found  that  the  path  from  k^ 
to  k2  would  have  been  rather  than  C^.  This  suggests  an  algorithm  for 
establishing  values  of  lniG^)  on  the  line  integral. 

The  path  C^,  as  defined,  is  always  on  the  original  top  sheet  of  G2 
so  tftat  and  B2  are  clearly  and  simply  defined  on  C^(with  the  hyperbolic 
branch  cuts).  Also  we  observed  that  lies  just  across  the  square  root 
branch  cut  from  which  simply  means  a  change  of  sign  for  (  £,/  B^). 
Therefore  we  have 


Tic  In  G  (6) 

J  -rnr~dS 

—M-jc  X 


[£n(a+b  B^B^  -  4n(a-b  B^B^JdS 


(3.4.11) 
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where  a  =  U2/ 


and  b  -  Uj/Otj+I^) . 


We  have  reduced  the  original  infinite  line  integral  to  a  finite  line 
integral.  Since  the  ratio  of  the  logarithm  arguments  in  (3.4.11)  is 
less  than  tt  ,  it  is  permissible  to  express  the  integral  as 


<TJc  in  G2(6) 

6  -  k 


Jin 


dS 


a  +  b 


a  -  b  62/fij_i 
5  -  k 


Jin 


d6 


a8. 


be; 


L?6!  "  bS2-i 


5  -  k 


d6 


-°°+jc 


(3.4.12) 


For  k^  on  the  inner  hyperbola  (  branch  cut  nearest  the  real  and  imagi¬ 
nary  6  axes) : 


oH-jc 

a 

-oofjC 


Hn 


U26i  +  UjPj 


U?2®1  ” 

6  -  k 


dd 


(3.4.13) 


6^  and  maintain  the  hyperbolic  definition  and  In  is  the  principal 
value.  For  k2  on  the  inner  hyperbola: 


. 

-°°+jc 


in 


vih  +  visz 


~  Ul®2-J 
5  -  k 


d6 


(3.4.14) 


To  consolidate  these  expressions  we  define  a  function  DIR(k^J: 


DI*(k  ) 


+1,  when  k^  is  on  the  inner  hyperbola 

(3.4.15) 

-1,  when  k^  is  on  the  outer  hyperbola 
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This  leaves  us  with  the  final  result 


“M-jc 


-00+ jC 


in  G2(6) 
£-k 


in 


d<5  -  [DIR] 


^2^1  +  ^*1^2 

1_W2®1  ~  Ul®2-t 
5^k 


d6 


(3.4.16) 


Note  that  ln( ( *-*2 2 ^ ^ ^ ^2^1” ^1^2^ ^  a  weH“^eliaved  function  for  6 

on  a  line  between  and  k2»  As  5  ■*  k^  ,  In  -*■  ±jTT  and  as  6  -*■  k2>  In  ■+•  0. 

The  function  (3.4.16)  is  analytic  in  the  region  around  the  inte¬ 
gration  path  except  at  k.^  and  k£.  In  fact,  the  integral  is  a  Cauchy 
type  integral  since  the  logarithm  is  analytic  along  the  integration  path 
except  at  k^  and  k2(see  Markushevich  (1977),  Gakhov  (1966)).  That  is, 
the  logarithm  function  is  continuous  on  the  integration  path.  The  value 
of  the  integral  is  an  analytic  function  of  k^  for  all  points  off  of  the 
integration  path  and  tends  uniformly  to  zero  for  |k^|  -*■  °°.  We  can  charac¬ 
terize  the  jump  in  the  value  of  the  integral  as  kx  traverses  the  path  of 
integration  using  the  Plemelj  formulas.  The  important  point  is  that  the 
integral  constructed  has  well-behaved  properties  in  the  variable  kx*  At 
this  point  we  observe  that  the  expression  for  G^,  given  by  the  formal 
factorization  formula,  was  asserted  to  hold  just  for  kx  such  that 
T  <  T  <  T+.  We  now  have  a  well-behaved  form  which  is  the  analytic  con¬ 
tinuation  of  Gt  onto  the  entire  k  plane.  Therefore  we  have  the  factori- 
2  x 

zation  of  the  original  function  G(kx),  (3.4.1),  given  as  follows: 
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w 


in 

DIR  f 

y261(6)  +  y1B2(6)- 

« 

u,  + 

u2S1(5)  -  y132(c)_ 

d6  • 

y2  *  ^l’kx  *  exp  < 

2ttJ  j 

C1 

6  -  k 

X 

(3.4 

where  8.^,  82  and  are  the  original  values  defined  in  Chapter  2, 

In  is  the  principal  value  of  the  logarithm  and  is  the  path  from  k^ 
to  specified  by  62  *  (k2-k^)t+k^  where  t  is  a  real  parameter  and  t  is 
contained  in  [0,1]. 

It  is  a  straightforward  exercise  to  formulate  the  integral  in  terms 
of  the  real  parameter  t,  which  is  used  to  characterize  6  on  C^.  Direct 
substitution  results  in  the  following: 


VV 


/k  -  k 
1  x 


exp 


2  2 
(k2'kl) 


DIR 


4-rrj 


U2B1(6)  +  y102(6) 
yzBi(6)  -  y1S2(<5)_ 
(6  -  kx)5 


dt 


<5  (t) 


(3.4.18) 


where  6(t)  *  /(k^-k|)t+k£  .using  the  principal  value  of  the  square  root. 

We  now  establish  a  useful  relation  between  G+(kx)  and  G_(kx).  It  is 

clear  that  G(k  )*G  (k  ) *G  (k  )  is  an  even  function  of  k  .  Therefore, 

X  +  X  -  X  X 


w.w 

G.cy 


G-(V 


for  G_(-kx),  G_(kx)  +  0 


(3.4.19) 


Each  side  of  the  above  equation  is  an  entire  function  since  G+(kx)/G_(-kx) 


a 
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is  analytic  for  T  >  T_  and  G+( -k^) /G_(kx)  is  analytic  for  T  <  r+.  By  the 

concept  of  analytic  continuation  we  assert  that  the  two  sides  are  equal 

for  all  k  on  the  top  sheet  and  are  entire  functions.  As  )k  j-*- <*>,  we  have 
x  x 

by  construction  that  G+(-kx)  since  the  exponential  goes  to  1  and, 

similarly,  G_(  k)  ■*  Therefore,  G.(-k  )/G  (k  )  -*■  1  as  |k  I**  00 .  Hence 

G+(-kx)/G_(kx)  is  an  entire  function,  bounded  at  infinity.  By  Liou- 
ville 's  theorem  from  complex  analysis,  we  have  that  G+(-k  ) /G_(k  )  is  a 
constant  for  all  kx  and  hence  is  equal  to  unity.  So  we  have  that 
G+(-kx)“  G_(kx)  and  we  see  that  this  symmetry  will  be  true  of  the  factor¬ 
ization  of  any  even  function  G(kx)  which  can  be  used  in  a  Wiener-Hopf 
analysis . 

We  note  that  G_(k  )  is  analytic  in  the  lower  half  of  the  k  plane. 
Since  G^jCk^)*  GCk^/Gjk^**  G(kx>/G+(-kx) ,  the  integral  for  G+(-kx>  in 
this  expression  will  converge  more  quickly  for  k^  near  the  integration 
path  from  k^  to  than  will  the  original  expression  for  G+(kx>.  The 
judicious  use  of  these  two  representations  for  G+(kx>  ensures  conver¬ 
gence  of  the  integral  which  is  to  be  evaluated  numerically. 

In  a  similar  manner  to  the  analysis  leading  to  (3.4.18),  we  can 
derive  the  integral  representation  for  G_(kx)  from  (3.4.7).  It  can  also 
be  shown  using  these  two  expressions  for  G+(kx)  and  G_(kx)  that,  in  fact, 
G^k^)  *G_(kJ£)*  GCk^.  It  is  also  clear  from  their  form  that  the  func¬ 
tions  are  analytic  in  the  upper  and  lower  portions  of  the  complex  kx 
plane  as  required.  The  symmetry  of  the  two  functions  (G+(k  )■  G_(-k  )) 
is  also  apparent  from  the  factorization  integrals. 

Before  leaving  the  subject  of  the  factorization  of  G,  it  is 
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appropriate  to  examine  the  location  of  the  roots  of  G  (k  )  and  G  (k  ). 

We  observed  that  the  integrals  in  the  exponentials  are  analytic  functions 

except  on  the  integration  paths  (k^,  k2)  and  (-k^,  -k2 ) .  They  are  bounded 

except  possibly  in  the  vicinity  of  those  lines.  Therefore  G+  and  G_ 

must  be  bounded  except  possibly  in  the  immediate  vicinity  of  those  lines 

(which  are  actually  branch  cuts  of  G  and  G  ).  Now  we  know  that  G(k  ) 

+  x 

has  four  roots  which  occur  in  symmetric  pairs  on  two  of  its  four  sheets. 

If  two  of  these  roots  occur  on  the  top  sheet  of  G  then  they  also  exist 
for  (G+*G  ).  However  the  arguments  of  the  exponential  factors 
which  appear  in  the  expressions  for  the  functions  G+  and  G  (using  the  al¬ 
tered  branch  cut  definitions)  must  be  bounded  and  hence  G+  and  G_  cannot 
be  zero  (except  possibly  at  the  branch  points  kx  3  ±k^).  This  is  a  de¬ 
monstration  that  G+  and  G_  have  no  roots  on  the  top  sheet  except  possibly 
at  the  branch  points  themselves  if  G  has  no  roots  on  the  top  sheet  with 
the  deformed  branch  cuts.  Therefore  the  roots  of  G,  G+  and  G_  lie  on 
the  sheets  adjacent  to  the  top  sheet.  As  noted  above  G+  and  G_  are 
related  functions  and  therefore  will  each  contain  one  root  on  the 
improper  sheet. 

In  closing  this  section  on  the  factorization  of  G,  we  reiterate  that 
a  knowledge  of  the  position  of  the  roots  of  G  and  construction  of  a  cut 
which  placed  them  off  of  the  top  sheet  is  crucial  to  the  formation  of 
the  simplified  integral  expression  (3.4,17)  and  is  subsequently  seen  to 
be  important  in  forming  an  efficiently  computable  solution  to  a  Wiener- 


Hopf  problem. 
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3.5  Decomposition 

In  the  analysis  below  we  carry  only  the  formal  functional  expressions 

G  (k  )  and  G  (k  )  to  represent  the  functions  derived  in  the  preceeding 
+  x  -  x 

section.  As  pointed  out  there,  the  analytical  behavior  of  the  factorized 
functions  is  simply  inferred  from  their  form,  and  they  are  straightfor¬ 
ward  to  evaluate  numerically.  Therefore,  having  performed  the  factor¬ 
ization,  we  may  now  proceed  with  the  decomposition  operation  specified 
in  Section  2.4  (equation  (2.4.9)).  The  function  to  be  decomposed  is 


n 


S(kx)  -  A_/G+  -  (k  cos0  +  k  )  •  G  ('kj  “  Vkx 

X  X  T  X 


SL(kJ  +  S_(kx)  (3.5.1) 


where 


(2 


|  jEosin0] 


fix 


The  decomposition  is  seen  directly.  G+(kx)  with  its  branch  cuts  must  be 

part  of  S  (k  )  while  the  other  singularity  due  to  the  root  of  (k  +k1cos0) 
+  x  x  - 

will  clearly  lie  in  S  ,  so  that 


s.(V  * 


G+(-k1cos0) 


(kx+  k^  cos0) 


(3.5.2) 


and 


w  - a 


w 


G+(-k^cos0) 


(k  +  k.cosS) 
x  1 


(3.5.3) 


Note  that  the  pole  of  S+  due  to  the  zero  of  the  denominator  of  the 
second  factor  is  cancelled  by  a  zero  in  the  first  factor.  An  application 
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tion  of  L 'Hospital  's  rule  shows  that  the  function  S+  is  bounded  at  this 
point.  Having  performed  the  decomposition,  we  proceed  with  the  Wiener- 
Hopf  analysis  as  follows: 


B+ 

4>  •  G  -  S  -  7-  +  S . 

-  "  "  G+ 


(3.5.4) 


We  observe  that  the  left  side  of  the  equation  is  an  analytic  function  in 

the  upper  k^  half  plane  (r  >  r  )  and  that  the  right  side  is  an  analytic 

function  in  the  lower  k  half  plane.  The  equality  stated  in  (3.5.4) 

holds  in  the  domain  (t  <  T  <  T  ).  Therefore  by  analytic  continuation 

we  can  assert  that  the  two  functions  are  analytic  continuations  for  the 

upper  and  lower  half  plane  of  an  entire  function  P(k^). 

Consider  the  asymptotic  behavior  of  G  (k  )  and  G  ( k  )  as  |  k  |-*-  00 . 

T  X  “  X  X 

Clearly  the  integrals  in  the  exponents  vanish  since  kx  appears  only 

I  i*5 

in  the  denominator  of  the  finite  integral.  Hence  G+  ^  G_  ^  |k  |  , 
and  from  the  expressions  for  the  decomposition  functions,  we  see  that 
the  asymptotic  behavior  is  as  follows: 


S  (k  )  and  S  (k  )  ^  Ik  I  as  k  -*-00. 
+  x  —X  X  X 


(3.5.5) 


If  we  examine  *Kx)  ■  ez(x),  we  have  from  the  edge  condition  that 
U 

<J)(x)  ^  x  .  It  can  be  shown  (Mittra  and  Lee(  1971a))  that  this  implies 
-3/2 

<Hk  )  "V  k  .  Therefore  we  have 
X  x 

P(k  )  -  $  *  G  -S  -  [(k  "3/2)(k  ~1/2)  -  k  l]  -*0  .  (3.5.6) 

x  -  -  -  X  X  X 

We  now  have  need  of  the  asymptotic  behavior  of  B^Ck^)  which  is  essen¬ 
tially  the  Fourier  transform  of  the  difference  current  on  the  conducting 
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half-plane  due  to  the  edge.  It  can  be  shown  that  B  (k  )  ^  k  .  There- 

+  x  x 

fore  we  see  that  P(k  )  ^  k  ^  for  Ik  [-*■  ®  and  hence  P(k  )  *  0  for  jk  I  ■  08 , 

XX  X  X  X 

We  have  shown  that  upon  completion  of  the  decomposition  step  we  can 

create  a  bounded  entire  function  PCk^).  By  Liouville's  theorem  we  know 

that  P(k  )  is  equal  to  a  constant.  However  we  know  that  P(k  )  *  0,  hence 

X  X 

within  the  strip  of  analyticity  of  the  Wiener-Hopf  problem  we  arrive  at 
(2.4.6). 

Therefore  we  obtain  the  following  result  for  the  transform  of  the 
total  electric  field  over  the  aperture  of  the  interface  problem: 


E  (k  )  - 

2  X 


2 

S-<V 

/!¥  r 

l 

ko 

"5T  5 

(k  +k,  cos 9)  G  (k  ) 
x  1  -  x  _ 

,  (3.5.7) 


where 


-j EQsin6 


it  n1  G+(-kjCOS'9) 
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CHAPTER  IV 


ASYMPTOTIC  ANALYSIS  OF  THE  INTERFACE  PROBLEM 


4.1  Fields  in  the  Transmission  Region  (y  <  0) 


4.1.1  General  Formulation 


As  a  result  of  solving  the  Wiener-Hopf  problem  for  the  interface 
geometry  in  Chapter  III,  we  have  determined  the  Fourier  transform  of  the 
total  tangential  electric  field  over  the  plane  y  =  0.  From  aperture  the¬ 
ory  we  know  that  this  can  be  directly  related  to  ez  anywhere  in  the  tran¬ 
smitted  region  (see  Collin  and  Zucker  (1969)).  The  electric  field  has 
only  Che  z  component  since  the  excitation  is  transverse  magnetic  to  the 
z  direction. 

Because  the  transform  variable,  k^,  is  normalized  to  the  free  space 
wavenumber  k^  m  ,  it  is  necessary  to  normalize  distance  to  free 

space:  x  =  k^'x'  *  (2TT*(distance  in  wavelengths))  where  x'  is  the  unnor¬ 
malized  distance.  This  ensures  that  k»x  *  k'*x'  ,  where  the  unprimed 
quantities  are  normalized  and  the  primed  ones  are  not. 

The  total  electric  field  e zin  the  transmitted  region  may  be  con¬ 
structed  to  be 


ez(x\y') 


/ZtF 


E  (k') 
z  x 


■j  (k 


k '  ■-oc+j  c  ' 
x  J 


dk'  for  y  <  0. 
x 


(4.1.1) 
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where  k'  =  (k ' )x-(k ' )y ,  k'  *  31  =  /k' *-k'  1 ,  o'  =  (x '  )x+(y '  )y . 
x  y  y  z  z  x 

Reverting  to  normalized  parameters  and  substituting  in  E  (k  )  from 

Z  X 

(3.5.7),  we  have 

-j(B2y+kxx) 

j  (FSTtSIeiTTio  £«  y<0  ,  (4.1.2) 

k  *-aafic  x  x 

X 

where  c  places  the  inversion  path  in  the  strip  of  analyticity  fixed  in 
the  Wiener-Hopf  procedure,  i.e.  that  of  (3.3.10)  with  £  given  in  (3.5.7). 

Equation  (4.1.2)  is  an  exact  field  representation  of  the  solution  to 
the  interface  problem  in  the  transmitted  region.  It  is  largely  a  formal 
result,  however,  since  it  manifests  the  computational  foraidability 
typical  of  inversions  of  continuous  spectra.  For  diffraction  computa¬ 
tion,  the  asymptotic  evaluation  of  (4.1.2)  for  large  Pis  quite  useful, 
however,  and  we  turn  our  attention  to  this  task. 

The  choice  of  the  branch  cut  due  to  B-,  in  E  (k  )  will  not  affect  the 

x  z  x 

convergence  of  the  inversion  integral  (since  does  not  occur  in  the 
exponential).  The  factor  G_  contains  cuts  due  to  B^  and  8.,.  Moving  the 
semi-infinite  cut  for  8^  does  not  affect  the  convergence  of  the  Fourier 
inversion  integral  since  G_  does  not  have  an  exponential  behavior.  It 
is  useful  later  to  exercise  this  freedom  in  the  asymptotic  evaluation  of 
(4.1.2).  We  can  analytically  continue  G  (k^)  into  the  extended  form 
with  the  semi-infinite  hyperbolic  branch  cuts  from  -k^  and  -k2  without 
affecting  the  integral  or  the  convergence  behavior  at  |k^|  -*■  °°.  In  so 
doing  we  may  expose  a  root  of  G  which  resides  on  the  exposed  portion  of 
the  lower  sheet.  Also  the  method  for  evaluating  G_(kx)  in  this  region  of 
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the  kx  plane  is  to  use  G_=  G/G+  where  we  use  the  original  definition  of 
G(k^).  This  continuation  of  G  causes  the  integrand  to  have  the  hypei — 
bolic  branch  cuts  extending  from  -k^  and  -k^  to  00  in  the  second  quadrant. 

4.1.2  Angular  Spectral  Mapping 

We  now  make  use  of  the  trigonometric  substitution  =  ^2  s^nCt>  where 
a  is  a  complex  variable,  to  map  the  integral  onto  the  a  domain.  Such  a 
mapping  is  conventionally  referred  to  as  determining  the  angular  spec¬ 
tral  domain,  a.  The  cuts  become  the  lines  in  Figure  4.1  which  are  asymp¬ 
totically  vertical  and  pass  through  a  *  +  tt/2  and  -  tt/2.  The  lower 
sheet  of  the  original  plane  is  mapped  into  the  regions  indicated  with 
crosshatching  in  Figure  4.1.  The  former  two-sheeted  plane  of  S^ik^) 
thus  becomes  the  periodic  plane  shown  in  Figure  4.1  with  the  role  of 
the  branch  cuts  replaced  by  lines  which  are  their  images  .  The  effect 
of  this  mapping  is  to  open  up  the  function  62“  k£  *  and  a 

saddle  point  in  the  integral  which  may  be  used  in  a  steepest  descent 
analysis.  Substituting  kx  =  k2  sin  a,  x  =  p  cos  (p,  y  *  p  sinp  for 
-it  <  $  <  0  into  (4.1.2),  we  obtain 


ez(p,<J>) 


K 


•  cosa 

G_(k2sina)  •  (k^sina  +  k^osQ) 


-j  (k.p)  •  sin(a-4>) 
e  da, 

(4.1.3) 


where  F  is  the  mapping  of  the  integration  path  with  kx  on  (-°°+jc  ,+0Q+jc ) . 
In  Figure  4.1  is  shown  the  portion  of  the  complex  a  plane  which  is 
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of  interest  and  the  position  of  several  key  features  of  the  kx  plane. 

For  simplicity  and  definiteness  we  take  the  inversion  path  to  be  the  real 
kx  axis  (this  is  correct  for  0  <  9  <  tt/2).  A  discussion  of  the  equations 
for  various  a  plane  contours  of  interest  is  given  in  Appendix  A.  As  not¬ 
ed  above,  the  boundary  of  the  top  sheet  in  the  plane,  which  is  defined 
by  the  branch  cuts  of  g  corresponds  in  the  a  plane  to  the  curve  repre¬ 
senting  the  image  of  the  real  k  axis  moved  to  intersect  the  real  a  axis 

X 

at  a  *  ±  rr/2  (see  Figure  4.1). 

The  next  stage  in  the  asymptotic  analysis  of  e  is  to  deform  the 

z 

integration  path  to  the  steepest  descent  path  (SDP)  through  the  saddle 
point.  In  the  process  of  deforming  to  the  SDP,  singularities  of  the 
integrand  may  be  crossed  and  must  be  considered  separately.  The  saddle 
point  is  characterized  as  a  point  in  the  a  plane  at  which  the  derivative 
ot  the  exponent  vanishes. 

—  [-j  k  p  sin  (a  -<$>)]  ■  -j  k  pcos(a  -<f>)  =0  (4.1.4) 

aa  s  l  s 

s 

as  *  >n  an  integer  (only  n  =  1  is  of  interest). 

Note  that  as  the  angle  of  observation  moves  from  tp  m  0  (on  the  lower  edge 
of  the  PEC  half-plane)  to  <p  M  -n  (on  the  lower  edge  of  the  interface  aper¬ 
ture),  Che  saddle  point  moves  from  a  *  tt/2  to  a  =*  -tt/2.  The  shape  of  SDP 

is  a  function  of  (cr4>)  and  hence  of  ( ct-a  ),  therefore  it  is  unaffected 

s 

as  the  observation  angle  changes.  That  is,  the  SDP' s  shape  is  fixed  for 
a  specific  problem  but  its  position  moves  so  that  it  crosses  the  real  a 
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axis  at  a  .  There  fore, once  the  curve  of  the  SDP  is  found,  only  the  re¬ 
gion  of  the  a  plane  between  the  integration  path  and  the  SDP  need  be 
examined  for  singularities  of  the  integrand.  Hence  a  graphical  proce¬ 
dure  can  be  used  to  locate  the  dominant  contributors  to  the  far  field. 

The  total  field  is  represented,  after  deformation  of  T  to  the  SDP, 
as  the  sum  of  an  integral  along  the  SDP  and  the  singularity  contrib¬ 
utions  swept  during  the  deformation.  For  ll^Pl  °°,  the  SDP  contribution 
is  asymptotic  to  a  term  containing  the  value  of  the  integrand  of  (4.1.3) 
at  the  saddle  point.  The  contributors  to  the  far  field  are  the  follow¬ 
ing: 


1.  SDP  integral  contribution  —  Present  in  all 
cases  and  is  equivalent  to  the  diffracted  ray 
for  the  half-plane  in  a  homogeneous  region, 

2.  Pole  due  to  the  root  of  (k^sina+k-^cosB)  — 
Contributes  when  swept  and  is  the  counterpart 
to  the  geometrical  optics  field  in  the  homoge¬ 
neous  half-plane  problem, 

3.  Pole  due  to  the  root  of  G_  —  Contributes  if 
crossed  in  the  deformation  to  the  SDP, 

4.  Branch  point  of  G_  at  k^*-k^  —  Contributes  when 
swept  and  denotes  a  term  arising  from  the  pres¬ 
ence  of  the  semi-infinite  opposing  half  space 
which  is  filled  with  a  contrasting  medium. 
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All  the  singularity  contributions  are  functionally  characterized  by  the 
integrand  evaluated  at  the  singularity  position  in  the  a  plane  : 


^ez^saddle  +  ^ez^g.o. 


'S' 6 


point 


pole 


1 

branch 

cut 


(ez} root 
of  G 


(4.1.5) 


When  the  terms  of  (4.1.5)  are  evaluated  asymptotically,  they  exhib¬ 
it  an  exponential  behavior  which  is  determined  by  the  position  of  each 
contributor  in  the  a  plane.  It  is  useful  to  determine  what  region  of 
the  a  plane  will  correspond  to  outwardly  decaying  waves  (and  therefore, 
physical  ones).  Clearly  these  waves  have  the  real  part  of  the  exponent 
less  than  zero, 

Re[-j  k2  sin(cH>)]  *  -k2isin(u-4>)cosh  v  +  k^cos (u-<J>)sinh  v  (4.1.6) 


where  a  ■  u+jv  and  k2"  ^2r^2i*  Therefore  the  boundary  of  the  a  plane 
corresponding  to  physical  waves  has 


k2i  tan(u  -*}>)*  k2r  tanh(v).  (4.1.7) 

We  observe  that  the  mapping  of  the  real  kx  axis  into  the  ct  plane  is 
given  by:  k^  tan(u)  *  k2r  tanh(v).  Hence  the  boundaries  of  physical 
waves  are  given  by  the  real  kx  axis  mapping  shifted  to  the  position  u  *  <P 
and  u  *  <p  +  7T  (the  region  from  tt/2  to  the  left  of  the  saddle  point  over 
to  tt/2  to  the  right  of  the  saddle  point). 

If  the  inversion  path  T  is  the  mapping  of  the  real  kx  axis  (which 
is  permissible  if  0  <  0  <  tt/2),  then  we  are  guaranteed  that  only 
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singularities  for  decaying  waves  can  be  crossed  for  all  observation  angles 
-7T  <  <p<  0.  See  Figure  4.2  for  the  case  <p  *  -  tt  (maximum  negative  position 
of  the  saddle  point).  On  the  other  hand,  if  tt/2  <  9  <  tt,  then  the  inte¬ 
gration  path  in  the  kx  plane  must  be  below  the  real  axis  and  map  into  a 
curve  to  the  right  of  the  mapping  of  the  real  axis.  Because  the  path  lies 
to  the  right  of  the  mapped  real  axis,  as  -*■  -tt  it  is  possible  to  sweep 
a  region  between  the  real  axis  mapping  and  the  path  of  integration  which 
causes  non-decaying  components.  However  we  observe  that  there  are  no 
zeros  of  G(a)  in  this  region  of  the  a  plane.  Hence  the  only  possible 
singularity  is  the  pole  at  (k2sina+k^cos8)  *  0.  Though  it  appears  to  be 
anomalous  that  a  field  contributor  can  becom  e  unbounded  as  one  moves  to  an 
infinite  distance  from  the  diffracting  edge  of  the  half-plane,  this  can 
occur  and  is  interpreted  as  follows.  For  tt/2  <  9,  the  incident  field 
along  the  interface  becomes  unbounded  as  s*-00.  For  ^  <  (which  is  the 
angle  that  the  refracted  wavefront  makes  with  the  half-plane) ,  the  obser¬ 
vation  direction  cuts  across  planes  of  growing  amplitude  as  p-*».  There¬ 
fore  this  nonphysical  behavior  of  the  pole  is  due  to  the  fact  that  a  plane 
wave  in  a  lossy  medium  does,  in  fact,  correspond  to  a  wave  which  decays 
as  it  propagates  in  from  p  «* 00 .  We  note  that  the  pole  arising  from  the 
root  of  (k2sina+kjC0s9  -  0)  corresponds  to  the  geometrical  optics  pole  in 
the  half-plane  diffraction  problem  in  homogeneous  space  and  gives  rise  to 
the  simple  transmitted  field  for  the  unobscured  portion  of  the  half  space 
below  the  interface.  It  is  therefore  clear  that,  with  the  exception  of 
the  geometrical  optics  field,  this  problem  can  only  give  rise  to  decaying, 
outgoing  fields  in  the  transmitted  region. 
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4.1.3  Identification  of  Asymptotically  Dominant  Singularit ies 


In  diffraction  problems  where  the  media  are  lossless,  the  exponen¬ 
tial  behavior  of  dominant  field  contributors  is  that  of  a  propagating 
wave  and  the  asymptotically  dominant  field  terms  are  identified  based  on 
the  exponential  decay  of  the  non-propagating  factor.  The  presence  of 
material  loss  introduces  a  new  complexity  in  that  all  of  the  singulari¬ 
ties  are  complex-valued  and  exponential  decay  occurs  concomitantly  with 
propagation.  We  therefore  must  identify  relative  rates  of  this  exponen¬ 
tial  decay  for  the  various  contributors  and  identify  the  dominant  terms. 

Since  all  observation  angles  contain  the  diffracted  field  due  to 
the  steepest  descent  integral,  it  is  useful  in  the  spirit  of  the  above 
analysis  to  compare  the  exponential  decay  of  swept  singularity  waves  to 
that  of  the  SDP  diffraction  wave  as  one  moves  radially  away  from  the  edge 
of  the  PEC  screen.  Saddle  point  analysis  allows  an  asymptotic  approxi¬ 
mation  of  the  SDP  integral  to  be  made  for  P  00  (see  Appendix  B)  yield¬ 
ing 


saddle 

point 


Ertk  ,sin0sin 
0  ob 


•PV 


-3kob° 


n^G+(-k^cos0)G+(-kQ^cos<}))  (k^cos^  +  k^cos0) 


(4.1.8) 


where  /  is  the  principal  value,  and  kob  ■  k2«  Note  that  this 

is  the  leading  term  in  the  asymptotic  expansion  of  the  saddlepoint  inte¬ 
gral  and  represents  the  contribution  to  the  radiation  field  (varies  as 
-1/2 

p  )  of  that  integral.  For  observation  angles  in  the  vicinity  of  the 
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interface,  the  dominant  contribution  is  a  lateral  wave  type  field  (var- 
-3/2 

ies  as  P  )  and  the  relative  significance  of  these  two  fields  is  dis¬ 
cussed  in  Appendix  B. 

The  dominant  far  field  contributor  (or  contributors)  is  determined 
by  identifying  regions  in  the  a  plane  which  have  greater  or  lesser  expo¬ 
nential  decay  than  Re  (-jk^p).  For  $  *  0,  this  region  is  bounded  by  a 
contour  defined  by 


Re[-jk2sin(a-4>)]  -  -k2i 


-  k^sin  u  •  cosh  v  +  k^cos  u*sinh  v  *  ~^2i 


(4.1.9) 


where  cc  ■  u+j  v.  After  algebraic  manipulations  we  have 


sinh  v 


1  ±  cos  I u  *  I  •  ll  + 


k2il! 


k.. 


2r 


^.inV 

2i  2r 


sin  u' 


(4.1.10) 


where  u'  ■  (u-u  ),  u  is  the  position  of  the  saddle  point.  Note  the 
s  s 

following  properties: 

Use  (  +  )  for  the  (v  >  0)  leg  (  bounded  by  -  (tt-u  )<  u'<  (tt-u  )) 

res  —  —  res 

Use  (-)  for  the  (v  <  0)  leg  (  bounded  by  -u  <  u'<  u  where 

res—  —  res 

u  »  tan  ^(k„./k_  ),  principal  value), 
res  2i  2r 

The  double  sign  on  the  square  root  reflects  the  fact  that  there  are  two 
curves  in  the  a  plane  on  which  the  exponential  decay  is  the  same  as  that 
of  the  saddle  point  contribution. 
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A  typical  plot  of  the  equal  decay  contours  derived  above  is  shown 
in  Figure  4.3.  It  is  seen  that  these  contours  divide  the  ot  plane  into 
four  regions  in  which  the  exponential  decay  is  either  greater  than  that 
at  the  saddle  point  or  less  than  it.  Clearly  the  SDP  will  lie  in  the  two 
regions  of  greater  decay  and  therefore  identify  them.  Several  observa¬ 
tions  are  appropriate.  While  the  expression  for  the  equal  decay  curves 
was  derived  for  $  *  0,  it  applies  for  any  0  with  u  replaced  by  (u  -  <J>) 
in  (4.1.9).  These  curves  cross  the  real  ot  axis  at  the  saddle  point  just 
as  the  SDP  does  and  therefore  the  entire  set  of  curves  (SDP,  equal  decay 
contours)  can  be  considered  a  template  which  slides  across  the  a  plane 
as  moves  from  0  to  —rr.  We  refer  to  this  collection  of  curves  as  the 
dominance  template. 

The  foregoing  analysis  allows  the  determination  of  the  most  signif¬ 
icant  contributors  in  a  particular  problem  for  any  specific  angle  of 
observation.  Most  of  the  a  plane  lies  in  the  region  of  greater  decay 
than  the  saddle  point  value  and  hence  any  singularities  occurring  there 
can  be  ignored  since  the  greater  exponential  decay  implies  asymptotic 
insignificance  as  p  +  ”.  On  the  other  hand,  if  a  singularity  lies  in 
the  region  of  lesser  decay,  then  the  saddle  point  value  can  be  ignored. 
Only  the  dominant  contribution  needs  to  be  considered  asymptotically  in 
the  far  field,  therefore  this  analysis  of  the  significance  of  the  vari¬ 
ous  contributors  is  fundamental  to  the  problem  for  observation  points 
which  are  far  removed  from  the  edge  of  the  PEC  half-plane.  Note  that  the 
equal  decay  contours  are  functions  only  of  1^  and  do  not  involve  k^. 
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'ZERO  DECAY 
n CONTOUR v\ 


Region  containing  singularities  which  represent  waves 
with  exponential  growth,  non-physical 


Region  containing  singularities  which  represent  waves 
with  greater  exponential  decay  than  the  saddle  point 


11  Region  containing  singularities  which  represent  waves 
1  with  less  exponential  decay  than  the  saddle  point 


Figure  4.3. 


Demarcation  of  regions  in  the  a  plane  about  the 
saddle  point  characterizing  the  exponential  decay 
rate  of  waves  associated  with  singularity 
locations.  (Dominance  contours). 


| 
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As  k^  ■*"  0  (medium  2  becomes  lossless),  one  equal  decay  contour  co¬ 
incides  with  the  real  ct  axis  in  the  vicinity  of  the  saddle  point.  Thus, 
in  the  lossless  case,  the  geometrical  optics  pole  and  the  diffracted 
(saddle  point)  contribution  manifest  the  same  exponential  decay  rates, 
i.e.  no  decay.  For  lossy  media  the  two  contributions  will  not  be  of  the 
same  asymptotic  order  since  the  exponential  decay  rates  will  be  differ¬ 
ent.  For  k^  <<  k^r  the  region  of  exponential  dominance  over  the  saddle 
point  contribution  collapses  to  a  narrowly  defined  region.  Then  it  is 
approximately  true  that  any  singularity  off  the  real  a  axis  can  be  ig¬ 
nored.  However  as  k  increases,  this  is  not  valid  and  the  significance 
of  the  various  regions  of  the  a  plane  must  be  considered. 

Another  interesting  question  is  to  determine  what  portion  of  the 
plane  gives  rise  to  outgoing  waves  and  what  portion  gives  rise  to  incom¬ 
ing  waves.  These  regions  are  separated  by  lines  of  zero  phase  given  by 

Im[-jk2sin(a-<}))  ]  *  -k2rsin (u-0)cosh  v  -  k^cos^-^Osinh  v  *  0  (4.1.11) 
or  for  0*0 

k2r 

tanh  v  *  -  - -  tan  u  (4.1.12) 

k2i 

This  is  the  equation  for  the  mapping  of  the  imaginary  k^  axis  onto  the  a 

plane  and,  in  particular,  it  describes  two  loci  shifted  so  that  they 

intersect  the  real  a  axis  at  u  *  u  ±  7T  /2  where  u  is  the  saddle  point 

s  s 

location.  The  region  between  these  two  lines  has  a  negative  phase  change 
as  p  increases  (outgoing  wave)  while  the  region  outside  of  these  lines 


has  a  positive  phase  change  (incoming  wave)..  This  wave  interpretat ion 
can  be  incorporated  into  the  dominance  template  discussed  above  as  shown 
in  Figure  4.4. 

Diffracted  waves  must  manifest  the  phase  progression  of  outgoing 
waves  only.  Incoming  waves  are  excluded  in  the  solution  process  due  to 
the  fact  that  the  swept  region  characterized  by  incoming  waves  cannot 
contain  singularities.  For  example,  Figure  4.5  shows  the  largest  region 
swept  between  the  real  axis  and  the  SDP  for  <j>  *  0  that  violate  outgoing 
wave  behavior.  A  similar  plot  can  be  made  for  $  =*  -  tt  with  the  saddle 
point  at  ag=«-Tr/2.  It  is  clear  that  the  swept  regions  of  outgoing  waves 
in  these  extreme  cases  lie  in  the  mappings  of  the  first  and  third  quad¬ 
rants  of  the  original  kx  plane.  There  are  no  singularities  of  the  inter¬ 
face  problem  occurring  in  those  quadrants. 

Summarizing  these  results,  we  see  that  given  a  specific  comp¬ 

lex  medium,  k^,  one  can  construct  a  dominance  template  which  centers  on 
the  saddle  point  and  divides  the  ct  plane  into  quadrants  of  exponential 
growth  or  decay  relative  to  the  diffraction  term  of  the  saddle  point  con¬ 
tribution.  The  dominant  contribution  can  be  determined  graphically  (or 
by  simple  algorithms  on  the  computer).  Since  the  only  possible  scattered 
waves  are  decaying,  outgoing  waves,  one  can  infer  that  no  singularities 
arise  in  regions  of  the  a  plane  which  would  violate  this  property. 

The  contribution  of  the  geometrical  optics  pole  which  arises  from 
the  factor  (k2sina+k^cos0)  in  the  denominator  of  the  inversion  integral 
is  given  as  follows  (see  Appendix  B): 
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Region  containing  singularities  which  represent  outwardly 
propagating  waves 


Region  containing  singularities  which  represent  inwardly 
propagating  waves 


Figure  4.4.  Dominance  contours  supplemented  with  phase  behavior. 


Kv  AXIS 


Figure 


Region  containing  singularities  which  represent  inwardly- 
propagating  waves 


4.5.  Example  of  the  integration  path  deformation  to  the 
steepest  descent  path  through  u  =  *r/2. 
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2Ensin6  -j (k  p)  sin(a  +  4  )  (4.1.13) 

[e  ]  =  - H _  e  06  r  ’ 

2  geora.  n^G(-k  cosn) 
optics 
pole 

where  a^  is  the  pole  location  and  kQb  =  k2 ,  G(a)  =  (3.4.1). 

Also  the  contribution  at  the  root  of  G_(a)  is  given  as  follows  (see  Ap¬ 
pendix  B) : 


[ez] root 
of  G 


2EQsin6 


n1G+(-k1cos8) 


G ,  (k  ,  sina  ) 
+  ob  r 


(k  .sina  +  k,cos9)G'(k  .sina  ) 
ob  r  1  ob  v 


x  e 


-j(kQb0)  sin(ar  +  |<j>|) 


(4.1.14) 


where  a  is  the  root  location,  k  ,  *  k„ ,  u  ,  *tL ,  U  , 

r  ob  2  ob  2  opp  1 


and 


G'  (a r) 


u2  -  v2 
ob  opp 


U  „ 
ob  opp 


tan  a. 


Hate  that  when  —  y2  no  root  of  G  occurs  (see  (3.4.4) 


4.1.4  Branch  Cut  Contribution 


In  order  to  consider  the  evaluation  of  the  branch  cut  contribution, 
it  is  necessary  to  review  (4.1.8,13,14)  and  consider  how  the  function 
G+(  a)  is  to  be  evaluated.  From  (3.4.17),  we  observe  that  the  factor 
/k^-kx  exhibits  a  simple  square  root  branch  cut  running  from  kx  =  k^ , 
while  the  exponential  factor  has  a  branch  cut  running  from  k-j_  t*  .  As 
noted  in  section  3.4,  the  factor  G2(kx)  (which  factorizes  positively  into 
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the  exponential  in 
R.  cuts  are  chosen 


another,  resulting 


(3.4.17))  is  a  two-sheeted  function;  and  if  the  g  and 

1 

to  coincide,  they  can  be  viewed  as  annihilating  each 
in  a  finite  branch  cut.  If  the  8  and  8  branch  cuts 


of  G*  (k^)  are  allowed  to  separate,  we  can  view  the  process  as  evolving 

as  shown  in  Figure  4.6.  If  the  result  shown  in  Figure  4.6c  is  combined 

with  G+(k  ),  the  derived  value  of  G  (k  )  uses  the  simple  hyperbolic 
lx  +  x 

values  for  3^  and  62*  However  it  is  difficult  to  see  how  to  perform  this 

continuation  in  light  of  the  integral  factorization  available  in  (3.4.17). 

The  simplest  procedure  is  to  view  G  (k  )  as  follows: 

+  x 


G  (k  ) 
+  x 


(3.4.18);for  Im(k  )>0  or  if  Im(k  )=0,  then  Re(k  )<0 

X  X  X 

(4.1.15) 


G(k  )/G(k  )  =  G(k  )/G  (-k  ) ;otherwise 

X  “  X  X  T  X 

where  G(k  )=(3«4,1)  and  G  (~k  )“(3,4*18) 
x  *  x 


Here  we  are  assuming  that  the  branch  cuts  for  8^  and  8^  do  not  cross  the 
real  kx  axis.  Observe  that  G+(-k  )  is  unambiguously  given  by  the  inte¬ 
gral  factorization  formula  for  k  in  the  lower  half  of  the  complex  k 
plane.  The  branch  cut  behavior  is  contained  in  the  factor  G(kx)  and  can 
be  dealt  with  conveniently  in  this  form  since  the  effect  of  shifting  the 
branch  cuts  is  manifest  directly  in  Gik^  as  the  functions  8^  and 

It  is  useful  to  recognize  this  method  of  evaluating  G+(kx)  since 
in  the  process  of  evaluating  the  branch  cut  contribution  to  the  radi¬ 
ation  integral  (4.1.3),  it  is  necessary  to  deform  the  8^  branch  cut  and 
the  various  singularity  contributors  can  be  swept  by  the  branch  cut. 

When  a  contributor  is  swept,  the  evaluation  procedure  given  in  (4,1.15) 
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allows  the  factor  3^  given  in  G(k^)  to  be  easily  determined. 

The  development  of  the  asymptotic  formula  for  the  branch  cut  contri¬ 
bution  for  =  -k^  follows  the  analysis  of  Felsen  and  Marcuvitz  (1973a). 
It  requires  that  the  branch  cut  from  kx  =  -k^  be  deformed  in  the  ot  plane 
so  that  it  lies  along  the  steepest  descent  path  (SDPb)  from  the  branch 
point.  It  must  be  emphasized  that  the  SDP^  changes  as  the  saddle  point 
moves.  The  locus  along  which  the  SDP^  lies  is  a  line  of  constant  phase 
of  the  exponential  and  points  on  this  line  have  the  greatest  possible 
change  in  the  magnitude  of  the  exponential  factor.  Hence  these  lines  are 
contours  in  the  a  plane  where  (4.1.11)  equals  a  constant  y  £  0.  Figure 
4.7  shows  such  a  set  of  phase  contours  in  the  a  plane.  Clearly  these 
phase  contours  move  with  the  saddle  point  (since  y  =-k2r  corresPonds  t0 
the  SDP).  Hence  the  SDP^  ,  which  must  lie  on  the  phase  contour  running 
through  the  branch  point  corresponding  to  kx  *  and  be  asymptotic  to 
the  SDP  through  the  saddle  point,  changes  continuously  as  the  saddle 
point  moves  with  varying  observation  angles.  A  detailed  discussion  of 
the  formulas  describing  these  SDP^  contours  and  the  effect  on  the  other 
asymptotic  contributors  is  given  in  Appendix  C.  We  summarize  below  the 
conclusions  given  there. 

In  deforming  the  branch  cut  in  the  a  plane  for  the  branch  point  cor¬ 
responding  to  (-k^)  to  a  line  along  the  SDP^,  we  may  be  forced  to  cross  a 
singularity  contributor,  ax<  The  result  is  that  the  value  of  8^(0^)  which 
we  must  use  in  the  asymptotic  formulas  is  given  by 

(  ®l(ax))  swept*” (  ®l^ax^  original  hyperbolic  definition' 


Figure  4.7.  Contours  of  constant  phase  in  the  a  plane 
about  the  saddle  point  location. 
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Looking  at  the  asymptotic  formulas  (4.1.8,14)  we  see  that  G+(-k^cos6) 
must  be  evaluated.  One  of  the  results  of  Appendix  C  is  that  the  geomet¬ 
rical  optics  pole  for  an  observation  point  in  the  transmission  region 
(y  <  0)  is  never  swept  by  the  8^  branch  cut.  Hence  G+(-k^cos9)  may  be 
evaluated  in  all  cases  using  the  original  definition  of  Bisection  2.2). 
Also,  in  the  process  of  deforming  the  branch  cut  for  6^  to  the  SDP^ 
configuration,  a  portion  of  the  lower  sheet  of  G_(a)  is  brought  onto  the 
top  sheet.  This  may  result  in  the  root  of  G_(  a)  being  brought  onto  the 
top  sheet,  a  detail  which  is  dealt  with  in  Appendix  C. 

Having  performed  this  deformation  we  observe  that  when  the  branch 
point  is  crossed,  the  entire  cut  is  crossed.  The  branch  cut  contrib¬ 
ution  is  now  amenable  to  an  asymptotic  analysis  somewhat  similar  to  the 
usual  steepest  descent  analysis  about  the  saddle  point.  The  rather 
complicated  result  is  stated  below: 


Sz  branch 
cut , 

-k 

opp 


£  y2  G  (-k  )/2irk  exp 

ob  +  opp  opp 


PV 


36' 


y  (k  -k,cos0)A  ,cosa.  ^V|k  ,  pcos(c*+  |$| ) 
opp  opp  1  ob  b  ob  b  1  1 


3/2 


x  e 


-j (k^p)  sin(afe  +  |  <fi  ( )  , 


(4.1.16) 


where  (k  oV  U0J  -  (k^y^;  (k  Qpp,  yQpp)  -(k  r  y^  ; /_?V  is  the  principal 
value;  9+  is  the  angle  of  the  SDP  ^  from  (corresponds  to  k  x  =*  ”k0pp) 
in  the  a  plane  relative  to  a  ray  from  ct^  in  the  positive  real  direction, 
where  9+  >  [Arg(k  Qpp)-Arg(k  ^  cos  (  )  ] ,  modulo  2fT(that  is,  if 
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(  Arg( k  )  -  Arg( k  ,*cos(a,))]  =  0  ,  then  9  <  9  <  9  +  2tt). 

opp  ob  b  g  g  +  g 

4.1.5  Uniform  Asymptotic  Expansion 

In  general  the  wavenumbers  for  the  two  media  involved  in  the  inter¬ 
face  problem  have  different  loss  tangents  and  as  a  result  will  not  lie 
on  a  common  tine  in  the  kx  plane  running  through  the  origin.  One  expects 
that  the  geometrical  optics  pole  will  lie  on  a  curve  in  the  a  plane  which 
only  intersects  the  saddle  point  locus  at  the  origin.  As  a  result  the 
asymptotic  expansion  of  the  radiation  integral  given  in  (4.1.8)  is  an  ade¬ 
quate  representation  for  observation  point  sufficiently  removed  from  the 
edge  of  the  PEC  half-plane. 

In  practice  one  would  like  to  have  an  asymptotic  form  which  is  appli¬ 
cable  at  smaller  radii.  This  necessitates  utilizing  a  uniform  asymptotic 
expansion  of  the  steepest  descent  integral.  The  resulting  expression  is 
given  as  follows: 
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f  6  1  ■ 
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point , 

g.o.  pole 


.fr  Eosln6  -jkob: 
'  T, — 6 


kob-si„l»l 


PV  G+(-kQbcos0)  •G+(-k1cos0)  •  (kobcos9  +  k  cos0) 
J  ob  J 


(4.1.17) 


/-jkJJ  2  •  G(k^cos8)  •  sin  ft-  — ) 


t  J 


>2p  e 


j2k0bOsin2(-El-i) 


G(k^cos0) 


-PV  . 


for  Im  [/-j2kQ^  *sin(  (ap-ag) /2)  ]<  0,  where  k  ob  *  k  ^  as  is  the  saddle 
point  location  in  the  a  plane  and  is  the  geometrical  optics  pole  loca¬ 
tion,  Q(y)  is  the  complementary  error  function  defined  by 

00 

Q(y)  -  f  e_X  dx  (4. 1.17. a) 


The  first  term  in  the  braces  in  (4.1.17)  is  the  original  non-uniform  ex- 
panson  contribution  (4.1.8). 
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4.1.6  Remarks  Concerning  the  Asymptotic  Analysis 

It  is  clear  from  the  form  of  these  asymptotic  constituents  to  the 
far  field  that  their  exponential  factors  depend  on  the  location  of  the 
singularities  causing  them  as  stated  earlier  regarding  the  dominance 
analysis  of  the  various  contributors.  Extreme  care  must  be  exercised  in 
this  analysis  to  determine  whether  or  not  the  root  of  G  (a)  occurs  on 
the  top  sheet  of  the  integrand  used  for  the  asymptotic  analysis.  The 
expression  for  the  saddle  point  contribution  given  in  (4.1.8)  is  the  non- 
uniform  asymptotic  version  where  no  singularities  are  assumed  to  lie  in 
the  vicinity  of  the  saddle  point.  For  the  transmitted  region  (as  long  as 
k^/k^  not  Purely  reaD  this  is  the  case.  In  the  homogeneous  half-plane 
problem,  the  geometrical  optics  pole  (a  function  of  9,  the  angle  of  inci¬ 
dence)  lies  on  the  real  a  axis  between  (-it/ 2,tt/2).  Therefore  as  the  ob¬ 
servation  angle  <p  approaches  the  shadow  boundary  (which  is  in  the  trans¬ 
mission  region),  the  saddle  point  expansion  becomes  unbounded  due  to  the 
breakdown  of  the  asymptotic  analysis.  To  avoid  this  singular  behavior, 
several  uniform  asymptotic  expansions  are  available  for  a  pole  in  the  vi¬ 
cinity  of  the  saddle  point.  However,  in  the  interface  problem  when  k2/k^  , 
the  contrast  ratio  of  the  media,  has  a  significant  imaginary  part  rela¬ 
tive  to  the  real  part,  then  the  geometrical  optics  pole  follows  a  locus 
as  a  function  of  9(the  angle  of  incidence)  which  is  off  the  real  a  axis 
and  crosses  it  at  a  *  0.  Nonetheless,  a  uniform  asymptotic  expansion  is 
given  in  (4.1.17)  to  allow  a  more  accurate  evaluation  of  the  saddle  point 
contribution  for  small  observation  radius. 
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An  interesting  physical  interpretation  can  be  given  to  the  smooth 
behavior  of  the  steepest  descent  integral  when  the  observation  angle 
brings  the  saddle  point  near  the  geometrical  optics  pole  (FelsenC 1982) ) . 

In  the  well-understood  problem  of  diffraction  by  a  half-plane  residing 
in  a  homogeneous  medium,  the  singular  behavior  of  the  saddle  point  con¬ 
tribution  (diffracted  wave)  in  the  vicinity  of  the  boundary  between  the 
region  lit  by  the  incident  wave  and  that  shadowed  by  the  half-plane 
(hence,  the  shadow  boundary)  clearly  indicates  the  physical  importance 
of  the  shadow  boundary.  Therefore,  the  fact  that  the  contribution  of 
the  diffracted  term  is  well-behaved  for  distant  observation  points  near 
the  position  of  the  shadow  boundary  in  the  interface  problem  means  that 
the  complex  contrast  ratio  causes  a  "smearing''  of  the  shadow  boundary  in¬ 
to  a  smooth  transition  from  the  lit  space  to  the  shadowed  space.  This 
is  true  except  for  the  case  of  normal  incidence  when  the  geometrical  op¬ 
tics  pole  is  located  at  a  -  0.  Therefore  the  complexity  of  uniform  asymp¬ 
totic  analysis  is  not  warranted  in  the  transmission  region  for  distant 
ooints  of  observation.  Equations  describing  the  locus  of  the  geometrical 
optics  pole  are  given  in  Appendix  A. 

It  is  important  to  realize  that  in  some  problems  where  the  medium 
of  the  transmitted  region  is  significantly  lossy,  it  is  quite  possible 
for  the  diffraction  term  arising  from  the  saddle  point  to  be  dominated 
by  one  of  the  terms  arising  from  the  singularities  mentioned  for  some  ob¬ 
servation  angles.  It  is  not  possible  to  discard  these  field  constitu¬ 
ents  until  the  problem  is  specified  and  the  positions  of  the  singularities 
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is  determined.  The  saddle  point  contribution  is  a  homogeneous  cylindri¬ 
cal  wave  which  appears  to  emanate  from  the  edge  of  the  PEC  half-plane. 

The  other  waves  are,  in  fact,  inhomogeneous  plane  waves  which  exist  in  a 
region  bounded  by  a  ray  (extending  from  the  edge  of  the  PEC  at  a  critical 
angle)  and  the  interface  plane  and  do  not  appear  to  come  from  the  edge. 
This  can  be  seen  by  considering  the  exponential 

-j  (k^p)  3 in  (01-4) 
e 


where  p  is  the  normalized  distance  from  the  edge,  <j>  is  the  observation 
angle  (-it  <  4  <  0)  and  a  is  the  position  of  the  mode  singularity.  Consi¬ 
dering  the  constant  phase  planes,  we  have  the  following: 


Im[-j(k  p)  sin(a-4) ]  *  Im[  jk  y  -  jk  x]  =  [k  y  -  ax]  =  0  (4.1.18) 

l  y  x  yr 

where  *  k2sina  *  a+jt  and  k2  *  ^2r~^lc2i ’ k2i— Recognizing 
that  the  direction  which  is  normal  to  a  constant  phase  plane  (given  by  9^) 
is  characterized  by  tanG^  *  l/(-dy/dx)  *  -  k^/a  ,  we  have 


tan 


k^tan  u  •  tanh  v  -  k^ 
k^tanh  v  +  k^tan  u 


(4.1.19) 


where  a  ■  u+jv,  9^  is  the  apparent  angle  of  phase  progression  relative 

to  the  interface  (9  ■  0  is  the  positive  x  direction)  and  does  not  refer 

P 

to  the  edge  of  the  half-plane. 

A  similar  derivation  for  the  constant  amplitude  planes,  ignoring 
algebraic  decay,  leads  to  an  angle  9  of  apparent  maximum  decay: 

3 


* 


I 
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From  the  solution  of  the  Wiener-Hopf  equation  we  have  the  Fourier 
transform  for  the  total  electric  field  along  the  interface.  An  expedi¬ 
ent  way  of  constructing  the  solution  in  the  upper  half  space  is  to  consi 
der  the  total  field  as  composed  of  three  terms:  the  incident  field,  the 
short  circuit  reflected  field,  and  the  field  due  to  the  scattering  from 
the  interface  structures.  The  sum  of  the  incident  field  and  the  short 
circuit  reflected  field  has  a  null  at  the  interface  by  definition.  If 
we  now  inc'ude  that  field  which  gives  the  correct  tangential  electric 
field  on  the  interface,  then  the  solution  is  complete.  The  field  which 

we  must  construct  below  is  the  field  due  to  e  on  the  interface.  This 

z 

is  simply  the  field  in  the  aperture  radiating  into  the  upper  region. 
Therefore  the  inversion  integral  is  closely  related  to  that  for  the  tran 
smission  region, 

-j(S  y +k  x) 

e  dk  (4.2.1) 

x 
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for  y  >0  where  c  is  in  the  strip  of  analyticity  of  the  Wiener-Hopf  formu¬ 
lation  and  where 


At  this  point  the  mapping  k^  *  k^sina  is  used  and  the  resulting  equation 
is  found  to  be 


e2(p,$) 


k^cosa 

G_(k^sina)  •  (k^sina  +k^cos8)  e 


-j (k  p) sin(a-H)) 


da  ,  (4.2.2) 


where  0  <  0  <  it . 

The  analysis  fol’^s  in  the  same  way  as  in  the  transmitted  case. 

The  asymptotic  forms  there  are  valid  here  also  if  (k  are  taken  as 

(k.u,)  and  (k  ,u  )  are  taken  as  (k„,y_).  We  observe  that  the  saddle 
1  1  opp  opp  2  2 

point  is  now  given  by  a  “  tt/2 ~<p,  so  that  we  may  take  a  =  ir/2-|4>|  as 

s  s 

characterizing  the  saddle  point  in  both  half  space  inversion  integrals. 

One  significant  difference  between  the  two  solutions  is  that  the  geometri¬ 
cal  optics  pole  now  moves  along  the  real  a  axis  between  (-tt/2, tt/2)  re¬ 
gardless  of  the  medium  parameters.  Hence  the  reflection  boundary  will  be 
a  distinct  physical  feature  of  the  problem  as  in  homogeneous  problem.  The 
problem  of  a  uniform  asymptotic  expansion  for  the  saddle  point  integral 
in  the  vicinity  of  the  geometrical  optics  pole  is  an  issue  here.  There- 


fore  the  uniform  asymptotic  expansion  given  in  (4.1.17)  wit  h  kob'  kl 
must  be  used  when  the  saddle  point  is  in  the  vicinity  of  the  geometrical 
optics  pole  to  avoid  the  singularity  resulting  from  the  non-uniform  expa- 
sion  given  in  (4.1.8). 

4.3  Validation  and  Sample  Data 

It  is  clear  that  the  solution  of  the  interface  problem  should  reduce 
to  the  diffraction  from  a  half-plane  in  a  homogeneous  space  as  the  media 
parameters  limit  to  each  other.  It  is  also  clear  that  the  numerical 
evaluation  of  the  factorization  integral  in  (3.4.18)  will  become  unstable 

as  the  two  media  limit  to  the  homogeneous  case  since  the  integration  path 
vanishes  while  the  integrand  becomes  unbounded.  As  a  test  of  the  algo¬ 
rithm,  the  case  shown  below  was  evaluated  on  a  computer: 

u  =  l.-j.OOl,  £  =  1.  resulting  in  k  =  l.-j. 000500  ,  ,, 

x  1  i  (4.3.1, 

1^2  =  l-~j-0003,  ^2  =  l.-j.OOl  resulting  in  k7  =  l.-j. 000650 

Incidence  angle,  9  =  45°  and  observation  radius,  c  =  5 
in  free  space. 

The  resulting  evaluation  of  (4.1.8)  for  the  cylindrical  wave  of  the 
diffracted  field  was  compared  with  the  Keller  diffraction  coefficient. 
Plots  of  the  two  diffracted  field  computations  were  found  to  overlay  one 
another  except  in  the  vicinity  of  the  observation  angle,  *>  =  130°.  The 
modulus  of  the  diffracted  field  for  the  parameters  given  above  is  plotted 
in  Figure  4.8  as  a  function  of  the  observation  angle,  taken  with  res¬ 
pect  to  the  illuminated  side  of  the  PEC  half-plane  where  0o<i<360°. 

The  anomaly  in  the  vicinity  of  i  »  180°  arises  from  the  vanishing 
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Figure  4.8.  Comparison  of  diffraction  coefficients 
from  the  interface  problem  with  the 
Keller  values.  The  two  agree  closely 
except  in  the  vicinity  of  <p  •  180°. 
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of  the  leading  term  in  the  asymptotic  expansion  of  the  steepest  descent 
integral  about  the  saddle  point  as  discussed  in  Appendix  B.  As  indicated 
there,  this  failure  of  the  leading  term  to  suffice  in  the  asymptotic  exp¬ 
ansion  is  localized  to  roughly  a  ±  2°  sector  about  the  interface  in  the 
present  example.  As  the  contrast  between  the  media  is  increased,  the  de¬ 
parture  from  the  Keller  diffraction  coefficient  around  <j>  *  180°  is  en¬ 
hanced.  This  behavior  is  seen  in  Figure  4.9. 

In  Figure  4.9  another  feature  of  the  diffracted  field  is  apparent, 
namely,  the  spur  occuring  at  <p  *  156°.  For  this  configuration  the  angle 
of  total  internal  reflection  is  4>  ■  24°.  The  spur  occurs  at  this  criti¬ 
cal  angle  and  always  lies  in  the  more  dense  medium  (i.e.,  that  medium 
having  the  larger  magnitude  of  the  real  part  of  the  complex  wave 
number).  Figure  4.10  is  the  same  configuration  as  in  Figure  4.9  but  view¬ 
ed  at  an  observation  radius  of  p  »  50  A  from  the  edge  of  the  PEC  half¬ 
plane.  It  is  evident  that  the  spur  is  diminished  as  the  observation  ra¬ 
dius  increases. 

The  spur  arises  when  the  two  media  have  almost  the  same  loss  tangent. 
In  that  case  the  branch  point  for  the  less  dense  medium  (-k2  for  the  ex¬ 
ample  shown  in  Figure  4.10)  approaches  the  line  in  the  kx  plane  which  con¬ 
nects  the  branch  point  of  the  denser  medium  (-k^)  to  the  origin.  This 
line  segment  maps  into  the  line  segment  between  a  *  -tt/2  and  a  ■  0  which 
is  the  locus  of  the  saddle  point.  In  essence,  the  spur  represents  the  ef¬ 
fect  of  the  branch  point  due  to  (— k 2^  on  the  saddle  point  integral  as  the 
saddle  point  passes  near  the  branch  point. 


•••  •  ' 


Modulus  of  diffraction  coefficient 
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The  validity  of  the  asymptotic  expansion  is  questionable  for  obser¬ 
vation  angles  near  the  spur  unless  P  is  large.  In  order  to  determine 
under  what  conditions  the  asymptotic  saddle  point  expansion  is  valid, 
it  is  necessary  to  establish  a  minimum  limit  for  the  observation  radius, 
p.  This  can  be  done  by  considering  the  value  of  the  exponential  enve¬ 
lope  of  the  integrand  in  the  saddle  point  integral.  We  may  consider  a 
ci  -cular  region  in  the  a  plane  which  is  centered  on  the  saddle  point 
and  has  the  branch  point  on  its  boundary  (the  region  of  validity  for  the 
Taylor's  series  expansion  of  the  integrand).  If  we  require  the  magnitude 
of  the  exponent  in  the  integrand  to  drop  to  one  tenth  of  its  value  at  the 
saddle  point  then  it  can  be  shown  that  the  smallest  observation  radius 
for  which  the  simple  asymptotic  value  of  the  saddle  point  integral  is 
valid  is  given  by 

p  %  . 73/r2  ,  (4.3.2) 

where  p  is  the  observation  radius  in  wavelengths  in  the  medium  containing 
the  observation  point,  r  is  the  separation  of  the  branch  point  and  the 
saddle  point  in  the  a  plane. 

It  must  be  said,  however,  that  the  spur  is  a  valid  phenomenon  in 
the  far  field  of  the  diffraction  pattern  and  represents  an  interaction 
between  the  radiation  portion  of  the  saddle  point  contribution  and  the 
branch  point  of  the  opposing  medium.  If  the  loss  tangents  of  the  two 
media  are  fixed  and  the  contrast  is  decreased,  the  spur  is  seen  to  move 
toward  4>  *  180°. 

A  typical  situation  is  shown  in  Figure  4.11,  where  all  the  terms  con- 


Observation  angle,  $  (degrees) 


-  Saddle  point  contribution  (radiation  field) 

-  Geometrical  optics  pole  contribution 

-  Branch  cut  contribution  (lateral  wave) 


Figure  4.11.  Typical  result  of  scattering  from  a  PEC 

half-plane  on  an  interface  between  dissimilar 
media.  (Parameters  —  ■  Uj  “  1*  ^  *  1.-J.0005, 

k,  -  1. 2248- j. 0041,  p  -  5  X  (free  space),  6  -  45°) 
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tributing  to  the  scattered  field  are  shown.  The  geometrical  optics  field 
contributes  in  the  region  bounded  by  the  reflection  and  shadow  boundaries 
of  geometrical  optics.  In  the  reflection  region  (y  >  0)  the  total  geome¬ 
trical  optics  field  is  composed  of  the  sum  of  the  short  circuit  field 
(that  due  to  the  incident  plane  wave  reflecting  off  of  the  interface  plane 
where  the  PEC  half-plane  has  been  extended  to  cover  the  aperture)  and  the 
scattered  field  of  the  geometrical  optics  pole,  in  the  portion  of  the 

reflection  region  which  lies  beyond  the  reflection  boundary  (tt— ©  <  <  9) 

these  two  constituents  must  sum  to  equal  the  field  that  would  exist  with 

no  half-plane  present.  As  a  result  the  field  of  the  geometrical  optics 
pole  can  have  a  magnitude  greater  than  one  in  a  lossless  configuration 
due  to  phasing  between  the  short  circuit  and  geometrical  optics  pole 
fields . 

In  the  transmitted  region  (y  <  0)  the  geometrical  optics  pole  is  the 
total  geometrical  optics  field.  Since  we  require  the  total  tangential 
fields  to  be  continuous  across  the  interface  and  the  short  circuit  field 
vanishes  at  the  interface,  it  is  seen  that  constituents  of  different 
radial  behavior  must  be  continuous  across  the  interface. 

The  decrease  in  the  geometrical  optics  field  as  $  moves  away  from 
the  interface  is  a  manifestation  of  the  greater  loss  in  medium  2  since 
the  ray  must  attenuate  more  quickly  as  the  ray  path  in  medium  2  increa¬ 
ses.  If  medium  1  is  much  more  lossy  than  medium  2,  it  is  possible  that 
the  geometrical  optics  pole  field  will  actually  increase  in  magnitude 
with  increasing  $  from  the  interface.  It  is  noted  that  when  medium  1  is 
more  dense  than  medium  2  and  the  angle  of  incidence  exceeds  the  critical 


angle  of  total  internal  reflection,  a  cutoff  effect  of  the  diffracted 
field  structure  is  seen  in  medium  2. 

The  lateral  wave  field  arises  from  the  branch  cut  integral  and  is 
seen  to  lie  in  the  more  dense  medium  with  its  peak  magnitude  occurring 
near  the  critical  angle  for  internal  reflection.  As  is  seen  in  Figure 
4.11,  it  seems  probable  that  the  mechanism  which  creates  the  spur  on  the 
saddle  point  diffraction  field  represents  some  interaction  with  the  la¬ 
teral  wave  as  both  occur  at  approximately  the  same  observation  angle. 

In  Figure  4.12  the  configuration  of  Figure  4.8  is  shown  versus 
increasing  radius.  It  is  clear  that  the  lateral  wave  suffers  a  larger 
algebraic  and  exponential  attenuation  than  the  other  contributors  and 
therefore  decays  more  quickly  with  distance. 

An  interpretation  of  the  phenomena  may  be  offered  to  suggest  the 
physical  mechanism  which  gives  rise  to  them.  The  lateral  wave  has  been 
described  in  numerous  places  (Felsen  and  Marcuvitz  ( 1973b, c),  Felsen 
(1967))  as  the  result  of  a  ray  which  propagates  along  the  media 
interface  in  the  less  dense  medium  and  leaks  energy  into  the  more  dense 
medium  at  the  critical  angle  due  to  refraction  of  the  ray.  In  these  in¬ 
terface  problems  one  may  visualize  the  initial  scattering  from  the  PEC 
half-plane  as  launching  a  homogeneous  cylindrical  wave  as  in  the  diffrac¬ 
tion  from  a  PEC  half-plane  in  a  homogeneous  medium.  However  the  ray  bun¬ 
dle  which  propagates  along  the  interface  experiences  precisely  the  same 
environment  as  a  lateral  wave  source  ray.  The  differing  phase  velocity 
of  the  two  media  refracts  energy  into  the  more  dense  medium  as  the  source 
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ray  propagates  and  depletes  the  energy  of  the  cylindrical  diffraction 
wave  in  the  vicinity  of  the  interface.  The  greater  the  contrast,  the 
greater  the  energy  loss  of  the  diffraction  field,  thereby  minimizing  the 
scattered  energy  in  the  direction  of  the  interface. 

This  viewpoint  also  explains  the  mechanism  by  which  the  lateral  wave 
is  launched.  One  asks  the  following  question:  how  can  the  diffracted 
field  propagate  with  e  J  1  in  the  upper  region  and  e  2  in  the  lower 
region  along  the  interface  and  represent  a  valid  solution  to  the  boundary 
value  problem?  Surely  even  an  approximate  solution  will  manifest 
continuity  of  the  tangential  electric  field  across  the  interface.  The 
solution  to  this  paradox  lies  in  the  realization  that  the  radiation  field 

Ac 

(which  varies  as  P  arising  from  the  asymptotic  evaluation  of  the  steep¬ 
est  descent  integral  at  the  saddle  point  vanishes  along  the  interface 

and  the  second  term  in  the  expansion  is  the  leading  term.  Recognizing 

-3/2 

that  both  the  saddle  point  and  branch  cut  contributions  vary  as  p  but 
with  a  propagation  behavior  of  e  J  1  and  e  J  2  respectively,  one  suspects 
that  the  saddle  point  contribution  for  region  1  will  equal  the  branch  cut 
contribution  for  region  2.  This  can  be  shown  to  be  the  case.  Therefore 
we  see  that  the  lateral  wave  contribution  is  essential  along  the  inter¬ 
face  to  match  the  saddle  point  contribution  for  the  other  medium  and  to 
provide  continuity  of  the  total  diffracted  field  across  the  interface. 

This  viewpoint  would  seem  to  explain  the  mathematical  necessity  of  the 
existence  of  the  lateral  ray  to  provide  continuity  of  phase  and  attenua¬ 
tion  across  the  interface  when  combined  with  the  saddle  point  contribu- 
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tion.  It  also  provides  a  useful  device  for  estimating  the  magnitude  of 
the  second  term  in  the  asymptotic  expansion  for  the  saddle  point  integral. 
That  is,  the  second  term  will  equal  the  lateral  wave  term.  Clearly  for 
p  *  5  X  this  contribution  is  small  as  seen  in  the  data  presented  below. 

In  Figures  4.13-4.18  sample  data  is  presented  for  a  number  of  possi¬ 
ble  media  combinations.  All  media  parameters  are  normalized  to  free  space 
In  all  cases  *  1>  the  permeability  of  free  space.  Plots  to  the 

left  are  associated  with  higher  loss  in  medium  2,  while  those  further 
down  are  associated  with  increasing  contrast  (that  is,  the  modulus  of 
£2^£1^  between  the  two  media  .  The  data  in  Figure  4.13  serve  as  a 
reference  case  where  the  angle  of  incidence  is  9  *  45°,  the  radius  of  ob¬ 
servation  is  p  *  5  Xin  free  space.  The  contrast  ratio  increases  by  a  fac¬ 
tor  of  three  as  one  moves  down  and  the  ratio  is  one  on  the  top  row.  It 
is  clear  that  increasing  contrast  causes  the  angle  of  refraction  of  the 
geometrical  optics  ray  to  increase  and  causes  the  lateral  wave  peak  to 
shift  away  from  the  interface.  Also  a  general  flattening  of  the  diffrac¬ 
tion  contribution  from  the  saddle  point  (particularly  about  the  interface) 
is  seen  with  increasing  contrast  .  Increasing  the  loss  of  medium  2  is  seen 
primarily  to  attenuate  the  fields  in  the  transmitted  region. 

Figure  4.14  is  identical  to  Figure  4,13  but  with  the  radius  of 
observation  increased  to  p  ■  20  wavelengths  in  free  space.  The  magni¬ 
tude  of  the  lateral  wave  in  the  case  of  no  contrast  (top  row)  is  ex¬ 
plained  by  noting  that  the  middle  picture  depicts  a  case  of  almost  equal 
loss  tangent  in  the  two  media.  Clearly  the  branch  point  which  gives  rise 


Relative  Loss:  Tan (6 


Relative  Loss:  Tan(6„)  =  0  Tan(6„)  =  Tan(6.)  Tan(<S0)  >  Tan(6.) 


case  as  Figure  4.13  with  p  =  20A,  free  space 


Relative  Loss:  Tan(60)  =  0  Tan(<$0)  =  Tan(6.)  Tan(6_)  >  Tan(6.) 


13  with  angle  of  Incidence 


Relative  Loss:  Tan(60)  =  0  Tan(6„)  =  Tan(6.)  Tan(6„)  >  Tan(6.) 


3  with  angle  of  incidence 
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to  the  lateral  wave  is  nearer  the  real  axis  of  the  a  plane  than  in  the 
two  surrounding  pictures  on  the  top  row  of  Figure  4.13.  Therefore  the 
lateral  wave- contribution  should  be  greater  at  a  given  observation  ra¬ 
dius.  Also  noticeable  is  the  fact  that  the  peak  of  the  saddle  point  term 
is  greater  on  the  reflection  boundary  than  on  the  shadow  boundary.  In 
fact,  the  peak  on  the  reflection  boundary  is  truly  a  singularity  (of  the 
non-uniform  asymptotic  expansion)  since  the  geometrical  optics  pole  re¬ 
sides  on  the  real  a  axis  for  observation  points  in  the  incident  medium. 

The  peak  on  the  shadow  boundary  is  not  truly  singular  since  the  geometri¬ 
cal  optics  pole  will  lie  on  the  a  axis  only  at  a  *  0  (<j>  *  90°)  for  pro¬ 
blems  involving  media  with  unequal  loss  tangents. 

Figures  4.15  and  4.16  are  given  for  comparison  with  Figure  4.13  and 
depict  the  same  case  but  with  the  incidence  angle  shifted  to  9  ■  30°  and 
60®  respectively.  As  anticipated  the  shift  in  the  reflection  and  shadow 
boundaries  is  apparent.  Figure  4.17  is  included  as  an  example  where 
medium  1  is  essentially  lossless  (k^  *  l.-j. 00005).  The  second  row  in 
this  figure  corresponds  to  a  value  of  such  that  Re(e2)  ■  3  which  is 
typical  of  many  plastic  dielectric  materials.  The  third  row  character¬ 
izes  medium  2  with  Re (£2)  ■  9  which  is  near  the  value  for  ceramic  dielec¬ 
trics.  In  contrast  to  the  earlier  cases,  Figure  4.18  depicts  a  case  where 

medium  1  is  relatively  dense  (k^  ■  1 . 732- j . 000029) .  The  top  row  of  this 
figure  represents  a  contrast  ratio  which  is  less  than  1  (le^/e^J  ■  1/3) 

as  opposed  to  equality  as  in  all  of  the  earlier  figures.  This  implies  a 

critical  angle  of  total  internal  reflection  of  about  55®  and  therefore  the 
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plot  of  this  figure  exhibits  the  phenomenon  of  total  reflection. 
Figure  4.19  depicts  the  case  given  by 


e  -  l.-j. 00001,  kx  -  l.-j. 000005 
e2  -  3. -j. 00001,  k2  -  1. 732- j. 000003 


p  *  5X  in  freespace 

Unlike  the  earlier  figures,  this  one  shows  the  diffraction  as  the  angle 

of  incidence  increases  from  15°  to  the  maximum  incidence  angle  for 

the  validity  of  the  Wiener-Hopf  analysis  (-Im(k  cos0  )  ■  Im(k_)) 

JL  max  z 

which  for  this  case  is  9  *  137°. 

max 

It  must  be  pointed  out  that  in  none  of  the  data  computed  has  the  pole 
due  to  the  root  of  G_  significantly  contributed  to  the  field  structure 
in  a  problem.  This  pole  is  associated  with  the  root  of  (3.4.1)  and  seems 
to  correspond  to  the  Zenneck  pole  contributed  since  they  both  are  roots 
of  the  same  equation.  This  pole  contribution  arises  only  when  the  branch 
cut  in  the  a  plane  exposes  the  pole.  No  clear  determination  has  been 
possible  to  extablish  the  actual  existence  of  the  contribution  of  this 
term  to  the  far  field  structure  in  a  problem. 


4.4.  Interpretation  of  Results 

The  intent  of  the  analysis  in  the  study  of  the  interface  problem  is 
to  provide  a  means  for  readily  computing  the  various  ray  optic 
contributors  which  arise  and  thus  a  means  to  apply  the  GTD  to 
configurations  that  contain  features  for  which  the  interface  problem  is 


me** 


Scattered  fields  for  the  case  described  by  (4.3.3)  with 
increasing  angle  of  incidence,  8.  p«5A,  free  space. 
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a  suitable  canonical  geometry.  It  is  desirable,  therefore,  to  discuss 
the  scattered  fields  in  physical  configuration  space  (that  is,  the 
spatial  domain)  and  describe  the  relative  magnitudes  of  the  ray 
contributors . 

A  ray  representation  of  the  contributing  fields  given  by  (4.1.8), 
(4.1.13)  and  (4.1.16)  is  shown  in  Figure  4.20  where  both  media  are 
assumed  to  be  lossless,  and  the  lower  medium  is  assumed  to  be  the  more 
dense.  That  the  Zenneck  pole  contribution  (4.1.14)  arises  is  a 
possibility,  too.  As  discussed  above,  however,  no  cases  have  been  found 
where  it  does  arise,  and  no  further  address  is  made  to  it  here. 

As  shown  in  Figure  4.20  the  geometrical  optics  field  in  the  upper 
half  space  is  given  by  the  simple  reflection  off  the  PEC  half-plane  in 
the  region  lying  between  the  half-plane  (at  <J>  *  0)  and  the  reflection 
boundary  (at  $  »  ^or  remaining  portion  of  the  upper  half 

space  the  geometrical  optics  field  is  simply  the  reflection  of  the 
incident  wave  off  the  media  interface.  In  the  lower  medium  the 
transmitted  wave  is  seen  to  bend  to  an  angle  $refr  and  represents  the 
refraction  of  the  incident  wave  at  the  media  interface. 

The  presence  of  the  edge  of  the  half-plane  launches  two  different 
waves:  an  edge-diffracted  wave  and  a  lateral  wave.  The  diffracted  wave, 
which  appears  as  a  cylindrical  wave  propagating  away  from  the  edge  of 
the  PEC,  is  seen  in  all  observation  directions.  When  the  media 
parameters  are  equal,  this  wave  is  described  by  the  Keller  diffraction 
coefficients  (see  Figure  4.8).  However,  as  the  contrast  is  increased, 
the  diffracted  wave  is  quite  different  from  the  "Keller"  wave  (for 
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example,  the  diffracted  wave  decreases  to  near  zero  as  the  observation 
angle  approaches  the  interface  while  the  "Keller"  wave  does  not,  see 
Figure  4.9)  . 

The  lateral  wave  is  launched  by  the  PEC  half-plane.  This  wave  is  an 

inhomogeneous  plane  wave  which  may  be  represented  by  a  ray  bundle 

emanating  from  the  edge  and  propagating  along  the  interface  in  the  less 

dense  medium,  shedding  energy  into  the  more  dense  medium  through 

refraction.  This  energy  thus  propagates  away  from  the  interface  at  the 

critical  angle  of  total  internal  reflection,  4>cr"  This  field  is 

present  in  the  region  of  space  between  the  boundary  ray  emanating  from 

the  edge  of  the  PEC  half-plane  at  an  angle  of  and  the  interface. 

The  lateral  wave  is  a  plane  wave  which  propagates  in  the  direction  <f>cr 

and  is  most  intense  along  the  boundary  ray  at  <)>cr.  At  a  fixed  radius 

the  magnitude  of  the  lateral  wave  decays  rapidly  as  the  angle  of 

observation  rotates  from  <t>  toward  the  interface.  This  lateral  wave 

cr 

is  a  unique  manifestation  in  problems  involving  an  interface  between 
media  and  does  not  occur  unless  dissimilar  media  are  present  in  a 
scattering  configuration.  Note  that  in  Figure  4.20  the  lateral  wave  is 
shown  propagating  into  the  lower  medium.  This  is  due  to  the  assumption 
that  the  lower  medium  is  more  dense  than  the  upper  medium.  Had  the 
reverse  been  assumed  then  the  lateral  wave  would  propagate  into  the 
upper  medium. 

Figure  4.20  shows  all  the  possible  ray  contributors  at  an 
observation  point.  It  is  clear  that  when  the  upper  medium  is  less 


dense,  any  observation  point  in  that  medium  will  receive  scattered 
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radiation  from  the  geometrical  optics  field  and  the  cylindrical  edge- 
diffracted  field.  The  geometrical  optics  field  will  dominate  in  the 
upper  half  space  when  lossless  media  are  involved  except  for  observation 
points  near  the  reflection  boundary  ^  where  the  geometrical  optics 
field  and  the  edge-diffracted  field  will  be  of  the  same  magnitude.  It 
is  noted  that  since  the  lower  medium  is  assumed  more  dense,  the  lateral 
wave  contribution  in  the  upper  medium  is  quite  small  and  evanescent  in 
nature,  decaying  away  from  the  interface. 

In  the  lower  medium,  it  is  seen  that  observation  points  in  the 
shadow  region  will  only  be  illuminated  by  the  edge-diffracted  field. 

For  observation  angles  between  the  interface  and  the  shadow  boundary, 
the  geometrical  optics  field  contributes  as  does  the  edge- 
diffracted  field  and  the  lateral  wave  contributes  in  the  sector  between 
the  interface  and  the  ray  at  the  critical  angle  of  total  internal 

reflection.  Throughout  most  of  the  "lit"  region  of  the  lower  medium  the 
dominant  contributor  is  the  geometrical  optics  field  and  the  least 
important  contributor  is  the  lateral  wave  field.  However,  the  lateral 
wave  dominates  the  others  in  a  small  angular  sector  bounded  by  the  ray 
at  ^cr*  I**  essence  a  pencil  beam  of  lateral  rays  appears  to  be 
emanating  from  the  edge  roughly  in  the  direction  4^  and  is  the 
dominant  field  contributor  there.  The  hierarchy  of  contributions  to  the 
total  field  is  summarized  in  the  table  in  Figure  4.21. 

It  is  emphasized  here  that  when  the  two  media  are  different,  the 
edge-diffracted  field  is  substantially  different  from  that  in  the 
homogeneous  problem.  The  lateral  wave  also  is  a  manifestation  of  the 


-o 

TO 

T3 

TJ 

TO 

l“l 

i-4 

r—4 

i—4 

0> 

<U 

<U 

0) 

TJ 

•  «4 

TO 

•*4 

TO 

T3 

•i4 

T3 

*44 

1—4 

U-l 

•-4 

*44 

1-4 

*44 

*44 

»-4 

Si 

a; 

0) 

a; 

CO 

•*4 

w 

•  -4 

CO 

•*4 

CO 

•  -J 

CO 

*  «4 

o 

*44 

u 

<44 

o 

*44 

O 

<44 

o 

<44 

•  —4 

«*4 

* 4 

•  -4 

4J 

TO 

4J 

A4 

T3 

44 

T3 

4J 

T3 

a 

0) 

CL 

0) 

a 

0) 

a 

0) 

a. 

V 

o 

44 

o 

<u 

o 

44 

o 

44 

o 

44 

U 

u 

o 

<D 

V 

o 

u 

*-4 

to 

*— 1 

to 

^4 

CO 

> 

> 

f—4 

CO 

i—4 

CO 

n j 

Li 

u 

to 

M 

C0 

to 

to 

L> 

CO 

Li 

u 

*44 

u 

*44 

*44 

3t 

o 

*44 

o 

*44 

■i4 

*44 

•-4 

*4-1 

•  ^ 

*44 

*  ^ 

*44 

*44 

u 

•  »4 

L 

•  >4 

Ll 

•*4 

f*4 

i— 4 

Li 

•H 

u 

'*4 

44 

to 

44 

T? 

44 

T3 

CO 

to 

44 

TO 

44 

TO 

a> 

1 

<y 

) 

a> 

l 

L» 

u 

<u 

1 

a) 

1 

gov 

O  00  J-I 

33  2 


u  e  <u 

•u  O  OX) 
0)  T3 
J  O 


CS  CO 


-O 

rO  <f 


\  * 
\ 

\ 


Significance  of  the  ray  contributors  in  the  scattered 
field  for  the  interface  problem.  The  lower  medium  is  more 
dense  in  (a)  and  the  upper  medium  is  more  dense  in  (b). 


Ill 


presence  of  the  interface.  Hence  the  field  structure  manifests  markedly 
different  features  from  the  classical  problem  of  scattering  from  a  half¬ 
plane  in  a  homogeneous  medium. 

When  losses  are  introduced  several  interesting  changes  occur.  Most 
marked  is  the  possibility  that  the  dominance  of  the  geometrical  optics 
field  over  the  edge-diffracted  field  in  the  lower  medium  can  reverse 
when  the  constrast  ratio  of  the  media  wave  numbers  (k^/k^)  attains  a 
significant  imaginary  part.  In  this  circumstance  the  locus  of  the 
geometrical  optics  pole  leaves  the  real  axis  in  the  angular  spectral 
plane  and  induces  the  exponential  dominance  of  the  edge-diffracted  field 
as  discussed  in  Section  4.1.3.  A  second  manifestation  of  the  presence 
of  media  loss  is  that  the  shadow  boundary  is  seen  to  "smear"  and  lose 
definition.  Further,  the  significance  of  the  lateral  wave  contribution 
is  strongly  dependent  on  the  contrast  ratio  (k^k^)  between  the  two 
media.  In  some  cases  the  lateral  wave  will  dominate  both  the 
geometrical  optics  wave  and  the  edge-diffracted  wave  for  observation 
points  sufficiently  close  to  the  interface. 


CHAPTER  V 


THE  SLAB  PROBLEM 


5.1  Statement  of  the  Problem 

The  structure  considered  in  this  chapter  is  shown  in  Figure  5.1.  A 
dielectric  slab  of  thickness  t’  with  constitutive  parameters  (y^.e^)  is 
embedded  in  a  surrounding  medium  with  constitutive  parameters  (y^.e^) 
where  and  are  real  and  6^  >  •  The  upper  surface  of  the  slab  lies 

in  the  plane  y'  =  0.  Residing  on  the  slab's  upper  surface  is  a  PEC  screen 
covering  the  half-plane  (x'  >_  0)  with  its  edge  at  x'  *  0.  A  plane  wave  is 
incident  from  the  upper  half-space  (y'  >  0)  and  propagates  toward  the  slab 
at  an  angle  0  with  respect  to  the  PEC  half-plane.  The  incident  wave  is 
assumed  to  be  polarized  transverse  magnetic  (TM)  to  the  z'  axis.  As  in 
the  interface  problem,  we  assume  no  variation  in  the  z'  direction  and  so 
this  is  a  strictly  two-dimensional  problem.  It  is  noted  that  only  loss¬ 
less  media  are  admitted.  This  stands  in  contrast  to  the  general  media 
considered  in  the  slab  problem. 

5.2  Wiener-Hopf  Formulat ion  o ver  the  PEC  Half-Plane 

For  this  problem  we  structure  the  integral  equation  over  the  PEC 
half-plane  rather  than  over  the  aperture  (as  in  the  interface  problem) 
to  demonstrate  the  features  of  this  formulation  and  how  they  contrast  to 
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INCIDENT 
PLANE  WAVE 


Figure  5.1.  Geometry  of  the  slab  problem. 
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the  features  of  the  aperture  field  formulation  used  to  solve  the  inter¬ 
face  problem.  The  resulting  integral  equation  is  an  electric  field 
integral  equation  (EFIE)  which  enforces  that  the  total  tangential  elec¬ 
tric  field  over  the  PEC  screen  be  zero  and  has  the  electric  current  on 
the  screen  as  the  unknown  function.  Since  we  are  considering  the  TM 
case  we  see  that  the  surface  current  on  the  plane  is  purely  z'  directed 
and  therefore  induces  only  the  z'  component  of  the  magnetic  vector  po¬ 
tential,  A2 .  To  formulate  the  integral  equation,  we  must  solve  two  asso¬ 
ciated  problems.  First,  the  problem  of  a  z'  directed  filamentary  source 
of  z'  directed  electric  current  located  on  the  upper  surface  of  the  slab 
(at  x'  *  0)  is  used  to  establish  the  Green's  function  for  the  integral 
equation.  Second,  the  problem  of  the  plane  wave  incident  on  the  slab  with 
the  PEC  screen  removed  is  used  to  establish  the  source  field  for  the  inte¬ 
gral  equation.  Having  established  these  two  functions,  we  form  the  EFIE, 


e^Cx’)  -  -€®*d(x') 

z  z 


j(x')k(x\x’)dx’ 

yy?  J  0  0  u 


(5.2.1) 


for  x'  e  (0,°°)  , 


NP  ,  ,  x  . 

where  ez  (x  )  is  the  total  field  present  due  to  a  plane  wave  impinging 
on  the  slab  in  the  absence  of  the  screen  and  k(x'  ,Xq  )  is  the  Green's 
function  relating  the  z' -directed  electric  field  on  the  interface  elec¬ 
tric  field  on  the  interface  to  a  z' -directed  current  filament  residing 


there.  It  is  convenient  to  consider  this  equation  directly  in  the  Fourier 


] 
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transform  domain,  therefore  we  extend  the  definition  of  the  functions  in 
the  manner  described  in  section  2.4  to  arrive  at 


a“(x’)  +  b“(x') 


w 

—  {  ♦i(*6)k(x,,*6)dJC0  for  x '£(—.•) 

/2tt  / 


(5.2.2) 


*  k 


NP, 


where  <)>,  (x')  *  j  (x‘)  and  a,  (x1)  M  e  (x')  for  x’  >  0.  Fourier  trans- 

T  J  Z  +  Z 

forming  this  equation  with  respect  to  x' ,  we  arrive  at 


A"(k')  +  Bu(k')  -  *"(k')  *  K(k’) 

T  X  *  X  t  X  X 


(5.2.3) 


for  k^  within  a  strip  where  all  the  functions  are  analytic.  We  now 
establish  the  Fourier  transform  of  the  known  functions,  K  and  A^. 

We  require  the  kernel  k(x' )  in  (5.2.2),  that  is,  the  Green's  func¬ 
tion  relating  the  field  in  the  >*'*  0  plane  to  a  z'-directed  filament  of 
current  in  that  plane.  Actually,  since  a  Wiener-Hopf  analysis  is  subse¬ 
quently  used,  we  need  only  the  Fourier  transform  of  this  kernel.  The 
analysis  for  the  transform  is  given  in  Whitmer  (1948)  and  results  in 

ax  ,  y’  >  0 


a(k;) 


a,  »  -t'  <  y'  <  0 


,  y  £  -t' 


(5.2.4) 


where 
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-j8^t'  ejB2tr 

-jS.'y’ 

p  -  "*■ 

•  \ 

-j 

e  r 

J 

^  6 

2/2? 

.  J 

(5. 2. 5. a) 


-jS’<y'+t')  e  J  (y '+t ' ) 

c 

-j 

e  r 

2/2? 

7 

(5.2.5.b) 


1  'I 


j± 


2/2? 


jB.’Cy’+c') 


(5.2.5-c) 


and 


6l-62 

8’  +  S’  * 


-jb:c  js  :t* 

(8*  +  8p  (r  e  2  -  j  e  2  ) 

The  Green's  function  which  we  require  for  (5.2.3)  is  given  by  a 
above  with  y'  *  0.  Hence  we  have 


wy' 

K(k;>  -  28[  F(k;> 


(5.2.6) 


where 


-0u>  ' 
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FOc;)  -  26J 


-jB't’  JB't* 

(SJ  -  Spe  -  (3{  +  B’)e 

2  -jS't'  2  jB't 

(6{  -  S')  e  1  -  (Bj_  +  S’)  e  1 


(5. 2. 7. a) 


or  equivalently, 


jB't’  +  Blt'cotCBlt') 

F(k')  -  28' t '  - = - - - 5 — - - - -  (5.2. 7. b) 

X  1  j[(BJt’r  +  (B^t’T]  +  2(B’t,)(B^t,)cot(B^t') 


An  interesting  observation  can  be  made  at  this  point  regarding 
(5.2.6).  The  poles  of  the  Green's  function  characterize  the  surface 
waves  for  the  dielectric  slab  while  the  roots  characterize  the  surface 
waves  of  the  problem  of  a  slab  with  the  upper  surface  covered  by  a  PEC 
screen.  The  problem  treated  here  is  a  composite  of  these  two  problems 
for  the  two  regions  x'  <  0.  One  expects  the  solution  for  these  two  re¬ 
gions  to  manifest  the  wave  structure  appropriate  to  each  region  in  the 
vicinity  of  the  slab.  As  is  seen  later,  one  result  of  the  Wiener-Hopf 
analysis  is  that  both  roots  and  poles  of  (5.2.6)  are  singularities  of 
the  integrand  of  the  radiation  integral  and  contribute  in  the  two  regions 
x  ^  0  as  expected. 

To  formulate  Au  (k' ) ,  which  is  the  Fourier  transform  of  the  incident 
+  x 

plane  wave  in  the  presence  of  the  slab  but  with  the  PEC  screen  removed, 
we  must  solve  the  requisite  boundary  value  problem.  The  analysis  is  rela¬ 
tively  straightforward  and  the  resulting  fields  are  given  below: 


9 


•  jk'(x'cos0  +  y'sin0)  jk'(x'cos8  -  y'sin9) 

EQ  e  1  +  A  e  1  ,  y '  0 

jk’(x'cos9  +  y'sinS)  jk'(x'cos 9  -  y'sinS) 

e  ■  is  e  +  C  e  ,  -t  <  y '  <  0 

2  —  — 


,  y'l-t' 

(5.2.8) 


where  k^cos(0)«k2  cosO^)  and  9,02  are  the  ray  directions  for  the  reflec 
ted  and  transmitted  ray  respectively  from  the  plane  of  the  slab.  The 
constants  in  (5.2.8)  are 

A  -  -EQ  +  B  +  C  (5. 2. 9. a) 


jk* (x'cos0  +  y'sin9) 
D  e  1 
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2B 

sine  SIn<V 


"i 


T 


sin02  j(k^sln0  +  k^sinO^t' 

~*r  * 


(5.2.9.d) 


where 


A  - 


sine  .  ,is*2 

2  jk2t'sin02 

Sine  sl”e2 

-jk^t'sin@2 

"i  "2  1 

6 

ni  n2 

6 

. 

s  - 

- 

ni  "  ’^rJ7iI  •  n2  “  \/z'l 


If  y’  is  set  to  zero  in  (5.2.8)  and  the  resulting  function  is  then  Four¬ 
ier  transformed  (leading  to  a  semi-infinite  integral  from  x'  *  0  to 
x'  ■  <*•),  the  resulting  transform  is 

_ Lfils _ 

+  x  /2tt  (k^  +  k^cos0) 


(5.2.10) 


where  Eq  is  the  magnitude  of  the  incident  plane  wave 

jk..'(x'cos8  +  y’sin0) 

einc  «  E.  e  1  and  7 1  -  F(k’  -  k'cos0). 

z  0  oxi 

Applying  the  normalizations  defined  in  section  2.5,  we  arrive  at  the 
following  form  for  the  Wiener-Hopf  equation: 

A+(kx)  +  B_(kx)  -  *+(kx)  G(kx)  (5.2.11) 


where 
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W 


J  Fe  Eo 


(k  +  k  cos6)  *  F9  *  F(kx  “  ^i^ose , c) 


B  (k^)  ■  unknown  function 


Vkx>  “  ^6  Jl  Jz(k0kx)t  WV  m  F{jz(x’)} 

c(kx)  -  F(kx)/e1 


5.3  Factorization  of  G(k ^ 

The  next  step  in  the  Wiener-Ropf  analysis  is  the  factorization  of 
G(k J  given  in  (5.2.11).  This  function  is  composed  of  two  factors  which 
are  considered  separately,  that  is,  1/  8^  and  FCk^.  Me  see  that  (1/8^) 
factorizes  by  inspection  as  in  (3.4.3).  Therefore  we  concentrate  on  the 
factorization  of  F  through  the  formal  integral  factorization  formula 
given  in  (3.4.6). 

To  meet  the  hypothesis  of  the  integral  factorization,  we  require 

that  F(k  J  -*■  1  uniformly  within  the  Wiener-Hopf  strip  of  analyticity  as 

| a |  -*■  °°,  where  k  *  a  +j  T.  We  observe  from  (5.2. 7. a)  that  as  |k  |  +  00  on 
x  x 

the  top  sheet  of  F  (where  8^  and  8 2  assume  the  original  definition  of  Sec¬ 
tion  2.2), 

Re(-j  g2t )  -*■  -  « 

and  therefore  F  ~  28^ 


0  -  (Sj+B^e 


+jS2t‘ 


2  +^2C 

0  -  (Bj+Bj)  e 


28, 


Sl+62 
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However  6,  ^  3„  ^  Jv?  *  ±  k  such  Chat  Ira(  S)  <  0  as  Ik  (  ■+<*>.  There - 
1  2  X  x  x 

fore  F  %  1,  as  |kx|-*-  “•  We  observe  that  even  if  we  move  to  along 

the  hyperbolic  82  branch  cut  (and  therefore  ReC-jg^t)  *  0),  again  F  +  1  as 

the  coefficient  (0  ~8_)  -*•  0.  It  is  clear  that  F  -*■  1  as  Ik  I-*-  °°  anywhere 
12  1  x1 

on  the  top  sheet. 

We  may  therefore  perform  the  integral  factorization  of  F  and  it  will 
apply  along  the  real  k  axis  since  the  Wiener-Hopf  strip  collapses  to  the 
real  axis.  It  is  necessary  to  consider  the  definition  of  ln(F)  which  ap¬ 
pears  in  the  factorization  formula  (3.4.6).  In  order  to  establish  an 

acceptable  definition,  we  must  consider  in  detail  the  behavior  of  F(k  ) 

x 

in  the  complex  k^  plane. 

The  details  of  the  behavior  of  F  are  essentially  the  same  as  those 
for  the  Green's  function  of  the  slab  problem  and  have  been  considered  by 
several  authors  (e.g.  Whitmer  (1948),  Barone  (1956),  Collin  (I960)).  The 
important  facts  concerning  F  are  presented  below  and  in  Appendix  D. 

The  points  k^  ■  ±  k^  will  be  branch  points  of  F  due  to  the  function 
8^.  However,  kx  ■  ±  k2  are  not  branch  points.  Observe  that  F  is,  in  fact, 
an  even  function  of  82*  If  we  consider  a  path  encircling  k2  in  the  kx 
plane,  we  observe  that  82  ”82  but  F  -*•  F.  One  can  make  the  observation 

here  that  the  Fourier  transforms  of  the  Green's  functions  for  the  slab 
problem  and  the  interface  problem  exhibit  branch  points  at  the  wave  num¬ 
bers  for  the  media  which  span  a  half-space  but  do  not  exhibit  branch 
points  for  the  media  whose  extent  in  the  y  direction  is  finite  (that  is, 
the  slab  medium).  Therefore  the  only  branch  points  of  F  lie  at  k  ■  ±  k^. 


In  addition  to  the  branch  point  singularities,  F  also  has  poles 
which  occur  at  the  roots  of  the  denominator  of  F.  These  roots  are  consi¬ 
dered  in  Appendix  D  and  they  are  found  to  lie  on  the  real  axis  between 
and  k£  (and  therefore,  by  symmetry,  between  -k2  and  -k^).  Interspersed 
alternately  between  the  poles  of  F  are  roots  of  F  due  to  roots  of  the  nu¬ 
merator  of  F  (also  constrained  to  the  real  kx  axis  between  (+k^,+k2)  and 
(-k2»~k^)).  The  number  of  poles  and  zeros  is  a  monotonically  increasing 
function  of  (a  t, normalized  values).  No  roots  or  poles  of  F  exist 

off  of  the  real  axis  between  (-k2  ,-k^)  and  between  (+k^  ,+k2).  If  the 
interval  (k^,  k^)  on  the  positive  real  k^  axis  is  considered,  the  typical 
distribution  of  the  singularities  of  ln[F]  is  shown  in  Figure  5.2. 

Consider  F  as  given  in  (5.2.7.b).  Since  k^  and  k^  lie  on  the  posi¬ 
tive  real  k  axis  with  k.  >  k,,  we  observe  that  values  of  k  in  the  inter- 
x  2  1  x 

val  (k^,  k^)  must  have  6^  31  -j  4>(kx)  and  S2  *  +  ty(kx)  where  <J>,  ip  are  po¬ 
sitive  real  functions.  Substituting  for  B^and  B2  in  (5.2.7.b),  we  see 
that  F  is  purely  real.  Therefore  if  we  consider  the  principal  value  of 
the  logarithm  for  ln(F),  the  branch  cuts  required  will  connect  respective 
roots  with  poles  as  shown  in  Figure  5.2  and  lie  on  the  real  kx  axis.  The 
logarithmic  branch  cuts  are  restricted  to  lie  on  (kpk2)  and  similarly 
on  (-k2,-k^).  We  recall  that  F  +  1  as  | kx | °°,  so  there  are  no  logarith¬ 
mic  branch  cuts  running  to  k  *  °°.  It  can  be  shown  that  a  path  in  the  k 
plane  which  encircles  a  pole  and  its  adjacent  root  will  not  suffer  an  in¬ 
crementing  of  the  argument  of  F  as  the  value  of  kx  moves  around  the  path. 
Therefore  we  have  adequately  described  ln(F)  with  the  branch  cuts 
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K|  K2 

FOR  a<ir/2  :  o**~*  □ 

ir/2<a<7r :  □  □ 

ir<a<3ir/2:  o— x  □ 

3ir/2<a<2ir :  o  o~*x  O'**  □ 


Branch  cut  of  the  logarithm 
0  Root  of  F  (and  branch  point  of  ln[F]) 

X  Pole  of  F  (and  branch  point  of  ln[F]) 


Figure  5.2.  Depiction  of  the  branch  cuts  of  ln[F]  due  to 
the  root  and  pole  progression  of  the  function 
F  in  (5.2.7)  as  the  thickness  of  the  slab 

increases  (a  ■  / t2) 
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described  above  using  the  principal  value  of  the  logarithm. 
The  integral  factorization  of  F  is  given  as 


W  =  exp 


air  1 


&n(F) 

6-k 


d  6 


\  for  T  =  Im(k  )  =  0 
x 


(5.3.1) 


where  the  path  of  integration  is  indented  below  6  =  as  shown  in  Figure 
5.3.  As  noted  above  FCk^)  -*■  1  uniformly  for  i  *>,  hence  we  can  deform 

Che  path  of  integration  of  (5.3.1)  to  that  shown  in  Figure  5.4.  It  is 
clear  that  wich  this  deformed  path 


d<5  -  i 27T  £n[F(k  )]  •  u(a)  +  j2irH(k  ) 
6-k  J  x  x 

x 


where  u(a) 


1  ,  a  >  0 

,  a  *  0  with  a»Re(kx) 

0  ,  a  <  0 


(5.3.2) 


and 


Therefore , 


H<V  ■  T 


wililldj 

6-k 

x 


F  (k  )  -  [ F (k  )]U(0)  •  exp{H(k  )}  for  Im(k  )  -  0.  (5.3.3) 

r  X  X  X  X 

Let  us  consider  the  behavior  of  (5.3.2)  for  kx  off  of  the  real 
axis.  Since  F  is  an  even  function,  the  integral  can  be  reduced  to 

^  ie  ■  2“x  j  d5  •  <5-3-4> 

-J-  X  0  X 

Since  ln(F(6))  ■+■  0  as  |6|-r  00  and  the  integration  path  is  removed  from 


81- 

INTEGRATION 
'  -PATH 


K> 


128 


any  singularities  of  the  integrand  (as  long  as  k^  is  off  of  the  imaginary 
5  axis),  it  is  seen  that 

<  M/|62-k2| 

for  a  suitably  chosen  constant  M.  Now  the  right  hand  side  of  (5.3.4)  can 

be  partitioned  as  /-j[2Im(k  )] 

X  + 

0  j[2Im(kx)] 

The  first  integral  is  clearly  convergent  for  any  k  and  for  the  second. 


Therefore  (5.3.4)  is  uniformly  convergent  and  is  a  regular  function  for 
k^  off  of  the  imaginary  6  axis  (see  Mittra  and  Lee(1971),  section  3.3). 
Hence  if  k^  is  any  point  in  the  complex  plane,  then  (5.3.2)  applies  and 
is  the  analytic  continuation  of  the  exponential  factor  of  F+(kx).  Also 
(5.3.3)  is  the  valid  expression  for  the  factorizat ion  of  F  for  all  points 


in  the  kx  plane. 

Observe  that  with  the  expression  given  in  (5.3.3),  it  is  unneces¬ 
sary  to  be  concerned  with  the  detailed  behavior  of  the  roots  and  poles  of 
the  function  F(kx)  in  evaluating  the  integral  H(kx).  It  is  a  straight¬ 
forward  operation  to  parameterize  (5.3.4)  using  <5  =*  j  cot(^t/2),  for  t  on 


the  interval  (0,  1),  to  arrive  at 


2k*  J  2k*  ] 


(5.3.5) 
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Observe  that  (5.3.5)  is  an  odd  function  of  kx  implying  that 
exp{(5.3.5)}  ■  1  for  kx  =  0. 

Some  interesting  observations  can  be  made  concerning  the  factoriza¬ 
tion  process  and  how  it  works.  In  electromagnetics  problems  one  often 
encounters  a  function  G(e)  which  must  be  factorized  and  satisfies  the 
following  conditions: 

1.  G  is  an  even  function  of  e. 

2.  The  Wiener-Hopf  strip  includes  the  real  e  axis. 

3.  G  meets  the  conditions  of  the  integral  factorization 
theorem. 


Then 


G+(£) 


exp  \ 


_i_  f  to  Gffl  { 

2nj  1  e-c 


(5.3.6) 


for  £  in  the  Wiener-Hopf  strip  and  Im(£)>  0.  When  £  is  extended  into  the 
rest  of  the  complex  plane,  we  have 


G+(e) 


exp 


1 

2TTj 


f  In  G  (6) 


S-£ 


l—  -as 


,  Im(£)  >  0 


(5.3.7) 


G(£)  •  exp 


1 

00 

In  G  (6) 

2TTj  , 

_  — <3 

6-e 

O  _ 

,  Im(£)  <  0 
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If  we  consider  only  the  integral  in  the  exponential,  then  we  observe  that 
it  must  be  an  odd  function  of  E  . 


H(-e) 


1 

2*j 


in  G(6) 


2-j  J 
(-«)■ 


in[G(-(-6) )  ] 
-[ (—6 ) — e ] 


[-d(-e)] 


1  in  G(-Y) 
2ffj  Y-e 

ym~ot> 


dy 


-  H(e) 


H(-e)  -  -  H(e) 


(5.3.8) 


Therefore  if  we  consider  an  £  such  that  Im(£)>  0,  then  we  have 


G+(-£)  -  G(-e)  •  exp [H(-e) ]  -  G(e)  -exp[-H(e)]  -  G(e)/G+(e).  (5.3.9) 

Therefore  if  we  accept  G  (£)  *  G+(-  (as  can  seen  directly  from  the 
formula  for  the  function  G  ),  then 

G+(e)  •  G_(e)  -  G(e) 

where  G+  is  analytic  in  the  upper  half-plane  and  G  is  analytic  in  the 
lower  half-plane.  Since  H(e;)  is  analytic  and  therefore  bounded,  G+  and 
G_  will  be  non-zero  in  their  respective  half-planes  of  analyticity. 

We  collect  below  the  results  for  the  factorization  of  the  slab 
problem: 

w  - 

G  <k  )  -  F  (k  ) //k  +  k  (5.3.10) 

-  X  -  X  lx 

where  the  square  root  functions  are  interpreted  as  in  section  2.2  and 
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FCk^)  *exp[H(lO],  Re(kx)>0 


W 


(5.3. 10. a) 


exp[H(k  )]  ,  Re(k  )<0 

X  X  7 


where 


H(kx)  -  -  -f 


£n{F[ 5 (t) ] } 


ti0  cos2(f)  +  fkxsin(^)]2 


(5.3. 10. b) 


6(0  -  j  cot  (~)  } 

F(k  )  -  (5.2. 7.b)  with  (k',t')+(k  ,t), 
and  x  xx 

F.(kx)  -  F+(-kx). 

To  complete  the  factorization  analysis,  it  is  necessary  to  determine 
the  asymptotic  behavior  of  G+(kx)  as  |kxl’*‘  ”•  w«  consider  first  H(kx)  as 
00  along  a  ray,  then 


H(V  -  -  T 


ci0  cos2  ^  +  fkxsintr^ 


- - - 


J  k2sin2  f— ) 
t-0  xSln  1  2  > 


(5.3.11) 


since  ln(F)  +  0  as  t  +  0.  Therefore 


F+(kx)  -  F(kx)-exp(H(kx))  -1*1 


(5.3.12) 


and  so 


W  ~  lkx' 


as  |k 


(5.3.13) 
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A  similar  result  applies  for  G_(kx). 


5.4  Decomposition 


Having  performed  the  factorization  of  the  function  G(kx)  in  (5.3.10), 
we  proceed  to  the  decomposition  of  the  function  A+/G_.  As  in  the  inter¬ 
face  problem,  this  step  can  be  carried  out  by  inspection.  The  resulting 
functions  are 

'F 


S  (k  ) 


-  x;  G  (kx)  •  (kx+k^cos9)  G  (-k^cos0)  •  (k^k^cos 


6) 


(5.4.1) 


and 


S.(k  ) 


+v  x;  G  (-k1cos0)  (k^k^cosO) 


(5.4.2) 


where  ^  *  jFq'Eq,  defined  in  (5.2.7)  and  (5.2.11).  The  formal  steps  to 
complete  the  Wiener-Hopf  analysis  are  similar  to  those  given  in  section 
3.5.  Therefore  the  Fourier  transform  of  the  current  on  the  PEC  screen 
is 


VkoV 


Htf-  $+(kx)  "  .rrr  |g 


|G  (-k^cos0) 


Gt (k  ) • (k  +k, cos0) 
4*  x  X  1 


(5.4.3) 


To  establish  the  radiation  integral  for  this  problem,  it  is  neces¬ 
sary  to  reconsider  the  initial  formulation.  The  procedure  was  to  parti¬ 
tion  the  fields  of  the  problem  into  two  components:  the  fields  due  to 
diffraction  of  the  incident  plane  wave  from  the  slab  inhomogeneity  with 
no  PEC  screen  present  and  the  fields  radiated  by  the  induced  electric 
currents  on  the  screen.  Therefore 


v  ta-  -  -  - 


133 


134 


When  converted  to  normalized  coordinates,  we  have 


_Rad ..  N 
Ez 


/H  k„ 


-jS.y 

$  (k  )*G(k  )e  where  y  *  k..y'.  (5.5.2) 

T  X  X 


Now  the  radiation  integral  itself  is  the  inverse  Fourier  transform  of 
(5.5.1), 


Rad  *  \  ,  v  .”1 

e  (x  »y  )  -  F 


-Radi 


1 

/2tt 


r  T5  j  -jk'x' 

'  •  e  x  dk'  .  (5.5.3) 

Z  X 


Again  reverting  to  normalized  coordinates, 


-f? 


-jCB.y+kx) 

W*G(Ve  dkx 


2tt 


QO 

• 

00 

k 

G  (-k^cosO) 


(k  +k.cos6)*G  (k  ) 
x  1  +  x 


-j(B  y  +k  x) 

G(k  )e  2  x  dk  . 
x  x 

(5.5.4) 


Substituting  for  ^  from  the  definition  following  (5.4.2),  we  arrive  at 

(5.5.5) 


Rad .  .  _ 

e  (x.y)  -  C 

z 


G  (k  )  -j  (B.y  +  k  x) 

Z  * _  a  L  X 

(k  +k,cos9)  x’ 


x  1 


where 


-j  Egk^sinQ  G+(-k^cos0) 


*  2tt 

and  (G+,  G  )  are  the  factorizations  discussed  in  section  5.3. 

The  integral  in  (5.5.5)  can  now  be  transformed  into  one  in  the 
angular  spectral  domain  using  the  substitution,  kx  ■  k^  sin(a),  to  give 
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Rad .  . 

e2  (P>4>) 


•  cos(a)  •G_(k1sina)  -jk  psin(ot-HJ)) 
(sina  +  coso)  e 

r 


(5.5.6) 


where  is  the  angle  of  observation  (0  <  <p  <  it).  The  saddle  point  for 

(5.5.6)  is  given  by  a  ■  (tt/2)-<P  and  the  equation  for  the  SDP  is  given 
by  (A. 5.2)  with  <J>  replaced  by  -<p  and  k^O;  that  is, 

cos(u-ug)  cosh(v)  ■  1  (5.5.7) 

where  ug  is  the  saddle  point  location  and  a  =  u+jv.  The  integrand  of 

(5.5.6)  contains  no  branch  cuts  since  the  mapping  removes  the  8^  branch 
cuts  and  S2  has  none.  Therefore  the  singularities  are  as  follows: 


1.  Geometrical  Optics  Pole  -  Root  of  (sina  +  cos0).  It  is 
clear  that  this  pole  lies  on  the  real  a  axis  between  — (it / 2) 
and  +(tr/2). 

2.  Poles  of  G_  -  Lie  on  the  line  u  *  -(tt/2),  v  <  0  between 
a  a  -(7r/2)  and  the  mapping  of  -k2. 


Forming  the  asymptotic  form  of  the  steepest  descent  integral  for 

(5.5.6)  through  the  saddle  point  proceeds  as  in  the  foregoing  analysis 
and  is  given  by 


sin9sin4>  G  (k^cosS)  G  (k^cos4>) 
(cos9  +  cos<|0 


-j^P 

(5.5.8) 


Similarly  the  analysis  leading  to  the  residue  contribution  of  the  geome¬ 
trical  optics  pole  is  straightforward  and  results  in 


r  Rad\ 

(e  ) 

z  g.o. 


jk  pcos(9+<J>) 

- j 2tt  C  G  (-k^cos9)e 


(5.5.9) 
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Observe  chat  the  geometrical  optics  contribution  suffers  no  decay  as  p 
increases  and  will  therefore  dominate  the  saddle  point  contribution  if  it 
is  present.  As  in  the  classical  Sommerfeld  half-plane  diffraction  prob¬ 
lem,  the  diffraction  term  is  only  important  in  the  shadow  region  of  the 
problem. 

To  formulate  the  residue  contribution  for  the  poles  of  G  (k^sina) , 
it  is  necessary  to  determine 


Lim  {(a-a  )  G  (k.sina)}  -  H(a  ) 
p  -  x  p 


(5.5.10) 


To  evaluate  (5.5.10),  G  is  reconstructed  as  G/G+  and  it  is  recognized 

that  G  is  analytic  and  well-behaved  at  a  • 

+  P 


G(k^sina) 

H(a  )  «  Lim  (a-a  ) 

P  +1 


sinaT  Ltm  { (ct*ap)  G(klsina)}  (5.5.11) 


G  (k.  sina_) 

+  1  P  a-*a 


Using  (5.2. 7. b)  for  G*F/8^,  we  have 


H(ap)  - 


2 ( j Si (ap)  +  82(ap)*cot(82(a p)*t) 


G, (k, sina  ) 

+  1  p 


x  Lim 


cr*a  j[B*(a)  +  B|(a)]  +  2  •  B1(a)  *  S2(a)  •  cot(B2(a)*t) 


(5.5.12) 
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Applying  L'Hospital's  rule  to  (5.5.12)  and  recognizing  that  the  denomi¬ 
nator  in  the  braces  in  (5.5.12)  is  zero  for  ct  ,  we  arrive  at 

P 


Lim 

tt+ct 


a-a 


[  1 


-J28  (2+j3  t)sina 

lp  Ip _ 2. 

sin2 (8.  t) 

2p 


(5.5.13) 


where  Blp  *  and  82p=  82  (a  ).  Using  (5.5.12)  and  (5.5.13),  the 

residue  contribution  for  the  surface  wave  poles  is 
b-j  cos(a  )  *H(a  )  -jk  Osin  (a  -K>) 


-  -2n 


5 


-jk  psin(a  H>) 

sin2(82pt)[jBlp+62pcot(82pt)]  e _ P 

k10  (2+jB.  t)  G  (k  sina  )*tan(a  )*(sin(a  )+cos0) 

1  Xp  Xp  +  1  p  p  p 


(5.5.14) 

Note  that  (5.5.14)  exhibits  the  exponential  decay  factor, 

-k,p  sin  <j>  sinh  tv  I 

where  oip  *  -(Tr/2)+jvp.  Therefore  at  large  P  the  surface  wave  contrib¬ 
utions  are  dominated  by  both  of  the  other  two  components.  This  is  not 
surprising  since  large  p  implies  an  observation  point  far  from  the  slab 
while  the  surface  waves  are  confined  to  the  region  near  the  slab's 
surface . 

Since  the  geometrical  optics  pole  is  present  in  every  problem  and 
lies  on  the  real  a  axis  between  -(tt/2)  and  (tt/2),  it  is  necessary  to  for¬ 
mulate  a  uniform  asymptotic  expansion  for  the  SDP  contribution  to  the  rad¬ 
iated  field  when  the  saddle  point  is  near  the  geometrical  optics  pole. 

This  result  is  given  below: 
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/  Rad'. 

(ez  SDP 
unif 


/tt  E,  e 


'j^P 


j2kxpcos2  (~) 


sgn(0-h}>-7T)G_(-k1cos0)e 


' 

jir/4 

/2k^D 

cos(— } 

eJ 

j 

J»/4 


/k^p 


/2  sin0  G_(kxcos<fc)  G  (-k^cosQ) 


(cos<t  +  COS0) 


SI  cos(^) 


(5.5.15) 


where  Q(y)  is  given  in  (4. 1.17. a). 

Before  moving  on  Co  the  fields  for  the  transmitted  region  (y  <  -  t), 
the  following  observation  can  be  made.  Recall  from  (5.4.4)  the  parti¬ 
tioning  of  the  fields  in  the  problem.  For  y  >  0  the  fields,  which  arise 
from  the  diffraction  of  the  incident  plane  wave  off  of  the  slab  with  the 
PEC  screen  removed,  can  be  written  as 


HP  Inc  ref  1  ,  JV*  coaS  +  y  sl»9> 

+  er  ‘  E0  e 


(5.5. 16. a) 


+  W1)e 


jkx(xcos0  -  ysin0) 


where  F  *  F(k  cos0)  with  F  defined  in  (5.2.7)  or  can  be  written  as 
0  1 
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jk  pcos?cos$ 

eNP  *  E_  e  [ j2sin(k  psincsinP) ] 

2  U  i 

jk  pcos  (0-hJ>) 

+  Egk^sin9*G(k^cos9)  e 


(5.5. 16. b) 


Observe  that  the  second  term  of  (5.5. 16. b)  is  canceled  by  the  geomet¬ 
rical  optics  contribution  (5.5.9)  in  the  shadow  region  for  observation 
angles  such  that  the  geometrical  optics  pole  is  crossed,  tp  >  (  tt— © ) . 

See  Figure  5.5.  The  term  in  (5.5. 16. b)  which  is  canceled  is  the  reflec¬ 
tion  which  would  be  seen  from  a  PEC  screen  covering  the  entire  upper  in¬ 
terface  of  the  slab.  Since  the  actual  screen  is  only  a  half-plane,  this 
geometri.cal  optics  field  is  forced  to  zero  in  the  shadow  region,  but  the 
correction  term  (compensating  between  the  slab  and  a  PEC  sheet)  repres¬ 
ented  by  Fg  remains. 


5.6  Asymptotic  Evaluation  of  the  Radiation  Integral  for  the  Transmission 
Region  (y  <  -t) 


To  begin  the  asymptotic  evaluation  of  the  radiation  integral,  it  is 
necessary  to  form  the  Fourier  transform  of  the  Green's  function  for  the 
transmission  region.  This  is  done  using  (5.2.4)  and  (5.2.6)  to  arrive 


tran ..  . 

G  (kx,y) 


[7  j8lC 

k  no  e 


<81"S2) 2e 


-jS2t 


-  (B^)  2e 


js2t 


JBjy 


(5.6.1) 


where  normalized  coordinates  are  used  and  m  /y^/ Eg  .  Now  the  radiated 
field  is  given  by  the  inverse  Fourier  transform  (in  normalized  coordinates^ 
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eRad(x,y)  =  F-1{J (k  )  •  [ T’H  Gtran(k  )]} 
2  x  x  ,y 


u 

00 

f 

jW 

GCran(k  ) 
x,y 

o  J 

! 

nr 

-j  k  x  ,, 

J  x  dk 


-2  5 


G_(kx)  82  e 


(k^+k^cosQ) [ (S1-82)e 


-i62t 


j82t 


-  (Sj+B^e  ] 


j(6  y-k  x) 

X  e  1  x  dk 

x 


(5.6.2) 


where  £  is  given  in  (5.5.5). 

Before  proceeding  we  note  the  singularities  of  the  integrand  of 
(5.6.2).  As  in  the  radiation  integral  for  the  reflection  region,  this 
integral  has  a  geometrical  optics  pole  at  the  root  of  (kx  +kjCOs9  ).  For 
values  of  k  with  X  =  Im(k  )  <  0  (or  if  X  =*  0,  then  -  =  Re(k,)  >  0), 

G_(kx)  is  analytic  and  non-zero.  Therefore  the  roots  of  the  bracketed  ex¬ 
pression  in  the  denominator  of  (5.6.2)  with  o  >  0  (which  lie  on  the  posi¬ 
tive  real  kx  axis  between  k^  and  k2)  are  poles  of  the  integrand.  While 
the  symmetric  roots  with  o  <  0  are  canceled  by  the  corresponding  roots  of 
G  .  However  for  o  <  0,  G_(kx)  has  poles  of  its  own.  Therefore  exclusive 
of  the  geometrical  optics  pole,  the  integrand  has  poles  at  the  poles  of 
G  (kx)  which  lie  between  -k2  and  -k^  and  poles  at  the  roots  of  the 
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denominator  factor  which  lie  between  k^  and  k ^ .  These  poles  correspond 
to  the  surface  waves  for  the  slab  problem. 

It  is  interesting  to  note  that  for  a  <  0  the  pole  locations  corre¬ 
spond  to  the  surface  waves  for  a  slab  of  thickness  t  and  must  propagate 
in  the  negative  x  direction.  Clearly  the  field,  which  is  launched  by 
the  presence  of  the  PEC  screen,  that  propagates  in  the  negative  x  direc¬ 
tion  from  the  edge  of  the  PEC  screen,  must  be  a  surface  wave  for  the 
slab.  On  the  other  hand  the  waves  moving  in  the  positive  x  direction 
must  characterize  the  surface  waves  for  a  slab  of  thickness  (2t)  due  to 
the  imaging  of  the  slab  through  the  PEC  screen.  In  fact,  the  pole  reso¬ 
nances  for  a  >  0,  which  characterize  the  waves  propagating  in  positive  x 
direction,  show  just  this  behavior  since  they  are  roots  of  G  and  hence 
roots  of  the  numerator  of  G.  These  details  are  discussed  in  Appendix  D. 

To  conveniently  evaluate  the  poles  of  (5.6.2)  exclusive  of  the 
geometrical  optics  pole,  two  alternative  forms  can  be  used. 


00  J6lc 

6Z  (x,y)  (k^k^cose)  B(kx)e 


G  (k  ) 

-  X 


j  (B.y  -  k  x) 


(5. 6. 3. a) 


-J8,11  JB-t 

[  (B1-82)e  -  (B^B^e  Z  ] 


,  a  >  0 


B(kx) 


Li 


♦  (5.6.3.b) 
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'jS„t 


,  ~j02C  JB-t 

[G+(kx)[(S1-62)2e  -  (B1+82)2e  ] 


,  a  <  0 
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Moving  to  the  spectral  domain  using  the  mapping,  k 


k^sin  n  ,  the 


radiation  integral  becomes 


e^d(p,P)  -  -2  £ 


jk^t  cosa 

r  8?'cosa*e 

- r— - - - 5-T —  B(k  sina) 

(sina  +  cos 0)  1 


-jk  psin(a-p) 
x  e  da 


(5.6.4) 


where  4>  is  the  observation  angle  with  respect  to  the  PEC  half-plane 
(-ir  <  <p  <  0).  The  integration  in  the  a  plane  is  shown  in  Figure  5.6 
(where  the  integration  path  is  offset  for  clarity). 

The  asymptotic  form  of  the  saddle  point  integration  is 


,  Rad.  f.  I  2 it 

(ez  W25(k^ 


^2^1C0S^^  ‘sin$*e 


(-jk^tsinp+j-j-) 


(co  sp  +  COS0) 


-B(k^cosp) 


e-Jkr  . 

(5.6.5) 


The  residue  contribution  for  the  geometrical  optics  pole  is 


,  jk  tcos0 

(e*3  )  *  j  (A.ttO  [g  (-k  cos0)  *B(-k  cos0)e  ] 

z  g.o.  l  1  1 


pole 


x  e 


jk1pcos(0-P) 


(5.6.6) 


The  geometrical  optics  pole  is  captured  in  the  deformation  of  the  inte¬ 
gration  path  for  p  >  (0-tt). 

Finally  the  residue  contribution  for  the  surface  wave  poles  must  be 
given.  Clearly  the  integration  path  must  indent  to  the  right  of  the 
poles  when  Re(ap)  <  0  and  to  the  left  of  the  poles  when  Re(a^)  >  0.  When 
the  residue  is  formed,  the  result  is 


Figure  5.6.  Integration  path  for  the  radiation  integral 
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(5.6.7) 


where 


SIGN  =  1  for  Re  (a  )  *  ± 

P  2 


The  limit  in  (5.6.7)  must  now  be  evaluated.  Consider  the  case  when 

Re(a  )*-(tt/2),  then 
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1 

(5.6.8) 

Applying  L'Hospital's  rule  to  the  limit  expression  and  simplifying, 
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(5.6.9) 


for  Re(a  )  ■  -(tt/2),  where  8,  *  8, (a  )  and  8-  *  6, (a  ).  A  similar  anal- 

p  lp  1  p  2p  2  p 

ysis  to  that  used  in  (5.6.8)  and  (5.6.9)  but  using  the  other  expansion  for 
B(k^  sin(a)),  given  following  (5.6.3),  results  in 
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(5.6.10) 


for  Re(  a  )  *  +(  tt/2)  . 

P 

As  in  the  asymptotic  analysis  for  the  reflection  region,  it  is 
necessary  to  give  the  uniform  asymptotic  expansion  for  the  saddle  point 
contribution  since  the  geometric  optics  pole  will  approach  the  saddle 
point  as  <t>  approaches  (-tt+0).  The  result  is  as  follows: 
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(5.6.11) 


where  Q(y)  is  given  in  (4. 2. 17a), g  -  g.(a  )  and  g_  ■  B„(a  ). 

2g  2  go  2s  2  s 

Finally  for  convenience,  the  portion  of  the  total  field  in  the 
transmitted  region  due  to  the  incident  plane  wave  with  the  PEC  screen 


removed  is 
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where  8^  *  k^sin(8),  82  *  k2sin(S2>  *  /k^-k|  cos2 (8) 
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CHAPTER  VI 


CONCLUSIONS 

The  objective  of  this  study  is  the  solution  for  the  scattering  of  an 
electromagnetic  plane  wave  from  a  PEC  half-plane  in  the  presence  of  pla¬ 
nar  media  discontinuities.  Two  specific  geometries  are  addressed.  The 
interface  problem  represents  a  basic  canonical  configuration  in  this 
class  of  problems  while  the  slab  problem  is  one  of  practical  interest. 

The  case  of  TM  incident  polarization  was  chosen  as  manifesting  all  the 
primary  characteristics  of  interest.  It  is  the  author's  belief  that  the 
TE  polarization  will  not  manifest  any  new  phenomena.  In  particular, 
the  functions  to  be  factorized  in  the  interface  problem  are  very  similar 
to  that  of  the  case  considered  in  Chapter  III.  Based  on  this  it  is  as¬ 
serted  that  the  analysis  presented  here  can  be  generalized  to  the  other 
polarizat ion. 

In  the  interface  problem  the  analysis  leading  to  the  radiation  inte¬ 
gral  for  the  TM  polarized  wave  is  given  along  with  an  efficient  form  of 
the  integral  factorization.  The  asymptotic  formulation  of  the  various 
far  field  contributors  is  given.  Numerical  data  is  presented  at  the  end 
of  Chapter  IV  to  demonstrate  the  behavior  of  the  various  field  contribu¬ 
tors  as  parameters  are  changed.  It  is  seen  that  a  diffracted  cylindri¬ 
cal  wave  is  launched  from  the  edge  of  the  half-plane  which  is  similar  to 
the  wave  described  by  the  Keller  diffraction  coefficient  as  the  media 
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contrast  is  eliminated.  An  interesting  feature  of  the  diffracted  wave 
for  the  interface  problem  is  the  decrease  in  the  radiation  field  as  the 
observation  angle  approaches  the  interface.  This  condition  is  seen  to  be 
necessary  for  the  total  solution  to  provide  a  continuous  tangential  elec¬ 
tric  field  across  the  interface. 

The  typical  geometrical  optics  field  for  the  interface  problem  is 
seen  to  reflect  and  refract  according  to  the  theory  of  geometrical  optics. 
A  peculiar  feature  of  the  interface  problem  is  the  launching  of  a  lateral 
wave  which  is  seen  to  have  its  peak  amplitude  in  the  more  dense  medium 
(which  is  defined  as  the  medium  having  the  greater  real  part  of  its  wave- 
number).  Finally,  the  theory  indicates  the  possibility  of  a  pole  term 
which  appears  to  be  similar  to  the  Zenneck  pole.  An  example  of  this  pole 
has  not  been  found  which  contributes  significantly  to  the  asymptotic 
field . 

The  analysis  of  the  slab  problem  is  presented  in  Chapter  V.  In  chis 
case  no  lateral  wave  is  present  but  surface  waves  are  launched  instead. 
The  asymptotic  analysis  has  been  performed  for  the  media  surrounding  the 
slab  but  is  not  presented  for  the  slab  itself  as  a  self-consistent  ex¬ 
pansion  has  not  yet  been  found.  No  data  is  presented  for  this  case. 

While  this  investigation  has  laid  the  foundations  for  study  of  the 
interface  and  slab  geometries,  it  has  by  no  means  exhausted  even  the 
theoretical  topics  of  interest.  In  the  interface  problem  the  issue  re¬ 
mains  concerning  the  possibility  of  a  Zenneck-like  pole  contribution. 
While  this  pole  occurs  in  the  TM  case  only  if  the  media  have  dissimiliar 
permeabi lit ies ,  it  may  occur  in  the  TE  case  for  media  which  have 


151 


differing  permittivities  and  is  therefore  of  importance  in  the  TE  case. 
Another  area  of  future  work  would  be  the  development  of  the  uniform 
asymptotic  expansion  of  the  saddle  point  integral  in  the  presence  of  a 
branch  point.  The  data  given  in  Chapter  IV  seems  to  imply  that  the  local 
disturbance  in  the  diffracted  cylindrical  wave  due  to  the  lateral  wave 
is  representative  of  the  field  structure  at  small  observation  radii  even 
though  the  formal  mathematics  of  asymptotic  theory  breaks  down.  It 
would  be  useful  to  confirm  this  conjecture.  Clearly,  another  avenue  of 
possible  work  would  be  to  perform  the  analysis  for  the  TE  case.  In  the 
slab  problem  numerical  implementation  remains  to  be  performed.  The  a- 
symptotic  expansion  of  the  field  structure  in  the  slab  also  must  be  per¬ 
formed  and  the  work  completed  indicates  that  an  explicit  modal  field 
will  result. 

The  form  of  the  solutions  is  structured  to  generate  the  ray-type 
diffraction  coefficients  needed  for  GTD  analysis.  While  the  diffracted 
cylindrical  wave  and  the  geometrical  optics  field  behave  as  simple  homo¬ 
geneous  waves,  the  lateral  wave  and  Zenneck-type  pole  wave  are  inhomogen¬ 
eous  plane  waves  and  will  require  an  extension  of  the  GTD  to  include 
scattering  of  complex  rays  in  order  to  be  utilized.  It  is  hoped  that 
this  work  provides  a  clear  demonstration  that  the  numerical  evaluation 
of  the  formal  Wiener-Hopf  factorization  integrals  is  a  feasible  process 
for  generating  ray-optic  solutions  to  complex  scattering  geometries. 

The  application  of  the  GTD  to  model  lap  junctions  between  metallic  and 
composite  dielectric  sheets  is  a  relatively  straightforward  result  of 


this  work. 


APPENDIX  A 


CHARACTERISTIC  EQUATIONS  OF  THE 

MAPPING  k.  -  k  sin  a 
x 


A.l  Introduction 


The  purpose  of  this  appendix  is  to  establish  the  principal  results 
which  determine  the  mapping  of  significant  curves  in  the  complex  plane 
onto  the  a  plane.  In  the  mapping  we  take  k  =  k^-jk^  as  the  complex  wave 
number  of  the  medium  of  observation  and  define  a  =  u  +  jv.  The  mapping 
is  written  as  follows: 

a  +  jt  *  (kr-jki)sin(u+jv)  -  (kr-jki)  (sin  u  cosh  v  +  jcosusinhv) 

CJ  *  k  sin  u  cosh  v  +  k^cos  u  sinh  v  ;  T  *  k^cos  u  sinh  v  -  k^sin  u  cosh  v  . 

(A. 1.1) 


Note  that  (A. 1.1)  maps  the  entire  plane  into  a  finite-width  strip  in 
the  a  plane  and  that  the  a  plane  is  periodic  in  its  real  part  with  peri¬ 
od  2tt.  These  results  are  based  on  the  assumption  that  both  media  are 
lossy  in  the  problem  treated. 

It  is  also  of  interest  to  provide  an  inverse  mapping  from  the  top 
sheet  of  the  k^  plane  onto  the  a  plane.  This  is  given  by 
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(A. 1.2) 

where  In  is  the  principal  value  of  the  logarithm,  i.e.  if  Arg(ln)  =  9 

then  -tt  <  9  <  7T  and  8  =  /k2-k2  as  defined  in  Section  2.2. 

A. 2  Real  and  Imaginary  kx  Axes 

By  setting  the  real  and  imaginary  parts  of  (A. 1.1)  to  zero,  we 
arrive  at  mappings  of  the  axes  in  the  k^  plane  onto  the  a  plane.  The 
real  axis  is  mapped  into  the  contour  defined  by 

tanh(v)  =  (k^yk^ )tan(u) .  (A. 2.1) 

From  this  formula  it  is  clear  that  the  mapping  of  the  real  kx  axis  (about 

a  ■  0)  lies  in  the  regions  of  the  a  plane  characterized  by  (u,v)  >  0  and 

(u,v)  <  0  for  any  media  to  be  considered.  Also  the  origin  in  [he  k 

x 

plane  corresponds  to  the  origin  in  the  a  plane. 

The  imaginary  k^  axis  maps  into 

tanh(v)  *  -(kj./k^)  tan(u)  (A. 2. 2) 

Therefore  the  imaginary  kx  axis  maps  into  the  two  regions  in  the  a  plane 
which  do  not  contain  the  real  axis  mapping. 

A. 3  Hyperbolic  Branch  Cuts  for  k 

The  branch  cuts  for  ±k  in  the  kx  plane  can  be  shown  to  "open  up"  so 

that  no  branch  cuts  for  8  are  required  in  the  a  plane.  The  multiplicity 

of  function  values  for  3  as  k  migrates  between  the  two  sheets  of  the  k 

x  °  x 
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plane  is  manifested  in  the  periodicity  in  Re(a)  under  the  mapping. 

However  the  mappings  of  the  branch  cuts  on  the  a  plane  delimit  the  boun¬ 
daries  of  the  top  sheet  of  the  function  */K2-k2  . 

r  x 

The  branch  of  6  =  /k^-k2  is  given  by  Im(8)  <  0.  Therefore  the 
branch  cut  lies  on  Im(8)  =  0. 

Im(S)  =  lm(/k2(l  -  sin2(a))  )  (A. 3.1) 

*  -k  sin(u) sinh( v)  -  k  cos(u)cosh(v)  =  0 
r  i 

Hence 

-(k./k  )  =  tan(u)tanh (v)  (A. 3. 2) 

i  r 

After  some  simple  manipulation  we  find  that  the  curves  in  the  ot  plane 
which  correspond  to  the  branch  cuts  from  -  k  in  the  k^  plane  are  given  by 

tanh(v)  *  (k^/k^)  tan(u  -  it/ 2 ).  (A. 3. 3) 

These  lines  correspond  to  the  mapping  of  the  real  kx  axis  onto  the  Ct 
plane  but  offset  by  -(tt/2)  along  the  u  axis.  This  region  corresponds  to 
the  top  sheet  of  the  functions  involving  8.  Therefore  these  lines  toge¬ 
ther  with  the  mappings  of  the  real  and  imaginary  k  axes  divide  the  top 
sheet  region  of  the  a  plane  into  the  mappings  of  the  four  quadrants  in 
the  kx  plane. 

A. 4  Hyperbolic  Branch  Cut  for  (-kQp) 

The  integrands  in  the  field  integrals  arising  in  the  interface 
problem  will  involve  the  factorizat ion  function  G_(kx)  which  contains 
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functions  with  the  original  hyperbolic  branch  cut  for  BQp  extending  from 
kx  ■  -kop .  (We  have  chosen  to  use  k  as  the  wave  number  of  the  observa¬ 
tion  medium,  therefore  kQp  is  the  wave  number  for  the  opposing  medium). 
Hence  it  is  necessary  to  determine  the  position  of  the  mapping  of  this 
branch  cut . 

As  in  the  case  of  the  branch  cuts  for  B,  the  8op  branch  cuts  lie 
along  the  curve  Im(Bop)  =  ImC^k^p-  k£)  =  0.  So  S^p  must  be  positive  and 
real  for  points  on  the  cut.  Therefore 

ImCB^p)  =  Im(k£p-  k2  sin2(ct))  =  0 

yields  the  defining  equation  for  the  locus  of  points  on  the  S0p  branch 
cut , 

(k  ) (k  . )  +  (k2-k2) (sin  u cos u  sinh  v  cosh  v)  * 
op,r  op, i  r  i 

(k  kJ)(sin2u  cosh2v  -  cos2u  sinh2v)  (A. 4.1) 

r  i 

where  kop  =  kop>r-jkop>i  and  k  *  kr-jki. 

After  algebraic  ma-ipul at  ions ,  we  arrive  at 

a2  -  2ab  sin (2u) sinh (2v)  +  b2sin2 (2u)sinh2 (2v)  *  cos2 (2u)cosh2 (2v) 
[b2sin2(2u)  -  cos2 (2u) ] sinh2 (2v)  -  [2ab  sin(2u) ] sinh(2v)  +  [a2-  cos2(2u)] 

(A.4.2) 

where  a  *  1  -  2  (kop , r‘ kop > j.)/( kr* kt) , 
and  b  *  (k^-  k^) /( 2* kr* k^ ) . 


Applying  the  quadratic  formula  and  simplifying: 
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sinh(2v) 


[(bsin(2u))z  -  cosz(2u)]  ]_ 


ab  sin(2u) 


±|cos(2u)!  /(a2-!)  +  (b2+l) sin2 (2u) 


where  the  correct  sign  for  the  square  root  must  be  chosen.  The 
for  choosing  the  appropriate  sign  of  (A. 4. 3)  is  straightforward 
results  are  given  below. 

For  | a |  >  1, 

sinh(2v)  ■  - - -  * 


[(bsin(2u))2  -  cos2(2u)] 


ab  sin(2u)  -  sgn(^r  +  u)  *sgn(-^r-  +  u) 


x  |cos(2u) |/Taz-1)  +  (b 2+1) sin 2 (2u) 


for  (-  T+  tan  (kVk^))  <  u  <  u 


(k  — k  f 

x  op 


nd  v  <  0. 


where  sgn  is  the  signum  function. 


For  a |  <  1 , 


cos(2u) 


a-b  sin(2u) sinh(2v) 
cosh  (2v) 


for  v  >  v  and  — ( tt / 2 )  <  u  <  0, 

—  (k  =-k  ; 

x  op 


(A. 4. 3) 


analysis 
and  the 


(A. 4. 4) 


(A. 4 . 5  .a) 


where 
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sin(2u)  - - - -  x 

[cosh2(2v)  +  b2sinh2(2v)] 

r  _ 7 

|ab  sinh(2v)  -  cosh(2v)  /-(a^+bz)  +  (bz+l)cosh2  (2v)  t 

(A.4.5.b) 

For  |a|  *  1  there  are  two  cases:  a  =  ±1.  If  a  =  -1,  then 

(k  k  . )  «  (k  k. )  and  the  branch  cut  for  6  will  coincide  with  at 
op ,  r  op ,  l  r  i  op 

least  a  portion  of  the  locus  of  points  on  the  8  cut  in  the  ct  plane  from 
a  m  — (  it/  2 ) .  In  fact,  either  (A. 4. 4)  or  (A. 4. 5. a)  may  be  used  recogniz¬ 
ing  that  (A. 4. 4)  gives  the  portion  of  the  branch  cut  for  v  <  0  and 
(A. 4. 5. a)  will  give  the  portion  for  v  going  to  +  °°.  On  the  other  hand 

if  a  =  1,  then  k  „  .  ■  0  and  the  8  cut  must  run  along  the  projection 
op ,  l  op 

of  the  positive  imaginary  kx  axis,  hence  we  must  use  (A. 4. 5. a). 

A. 5  Steepest  Descent  Path 

The  characteristic  curves  discussed  in  previous  sections  of  this 
appendix  are  fixed  for  a  given  problem  once  the  media  are  specified.  As 
stated  in  Section  4.1.2,  the  saddle  point  for  the  field  integral  is  found 
as  a  root  of  the  derivative  of  the  exponent  in  the  integrand  and  is  found 
to  be  given  as  At  the  saddle  point  the  exponent  is 

7T 

-jkp  sin(as-<j>)  ■  -jkp  sinCy)  *  -jkp.  (A. 5.1) 

Now  the  steepest  descent  path  is  determined  such  that  the  imaginary  part 
of  the  exponential  for  all  points  on  the  path  is  equal  to  the  saddle 
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point  value,  i.e. 

Im(-jk  sin(a-<f))  *  -(kr  sin(u-$)cosh(v)+k^  cos(u-4>)sinh(v))  *  -k 
or  equivalently 

1  -  sin(u-<J>)cosh(v)  =  (k^/k^cosCu-^sinhCv) .  (A. 5. 2) 

This  expression  may  be  reduced  to  an  explicit  formula  for  the  steepest 
descent  path. 

By  manipulating  (A. 5.2)  and  using  the  quadratic  formula,  one  arrives 
at 


tanh(v)  *  (-B±/B2-AD  )/A  ,  for  v  <  0, 
where  A  ■  1  +  (k^cos ( u-<j>) /k^ ) 2  , 

B  -  y  sin(2(u-$))/2, 

D  *  -  cos2(u-4>) , 


(A.5.3) 


We  observe  here  that  as  v  "*■  +  *>,  the  asymptotic  form  for  the  original 
expression  is: 


-  tan(u-0)  'v  (ki/k^ )tanh(v)  n.  (k^/k^).  (A. 5. 4) 

Using  the  trigonometric  identity  tan(u-u  +tt/2)  *  -l/tan(u-u  )  where  u 

s  s  s 

is  the  position  of  the  saddle  point  (u  ■  <i>+  t/2),  we  arrive  at  the 

s 

following: 

tan(u-u  )  ■  l/(k./k  )  -*■  u-u  ■  ir/ 2~tan  ^(k  /k  ). 

s  i  r  s  i  r 


(A. 5. 5) 
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So  as  v  +  +«  we  see  chat  u  -*■  ug  +  (TT/2-tan_^(k^/kr) ) .  Similarly  as  v  -*■ 

u  ug-(^/2-tan  (k^/k^)).  Using  these  values  it  can  be  seen  that  when 

the  saddle  point  lies  at  u  ■  ~^/2,  the  asymptotes  of  the  steepest  descent 

path  are  the  same  as  the  asymptotes  of  the  respective  boundaries  of  the 

top  sheet  (mappings  of  the  hyperbolic  branch  cuts  from  the  points 

k  *  ±k ) . 
x 

One  final  point  of  interest,  which  is  important  in  Che  evaluation 

of  the  steepest  descent  integral  through  the  saddle  point,  is  the  value 

of  the  slope  of  the  path  at  the  saddle  point.  If  we  consider  points  in 

1  0 

the  vicinity  of  the  saddle  point  (  a  *  ug  +  peJ  ),  then  the  exponent  of 

the  integrand  of  the  field  integral  will  be  given  by 

(A. 5. 6) 


g  *  -jk  sin(a-<())  ^  -jk  +( jk/2! )(a-ag)2+. . . 

■  -j (k^-jk^+i  j/2)  (kr~jki)p2(cos(  20)+j  sin(29))  +  ... 

■  (~kj+( p2/2) (k^  cos(28)-kr  sin(20)>) 

+j (-kr+( p2/2) (kf  cos( 29)+kisin( 20) ) ) 


Now  since  ( Im(g)  )  ga{j^g  *  -kr  »  Chen  we  must  require 
point 


(k^  cos(29)+k^  sin(20))  ■  0. 


(A. 5. 7) 


Hence : 


(tan(29))  ,,  ■  -(k  /k.)  (A. 5. 8) 

saddle  r  i 

point 

which  is  the  slope  of  the  steepest  descent  path  at  the  saddle  point. 
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A. 6  Path  of  Integration 

The  field  integral  involves  an  integration  path  which  lies  in  the 
Wiener-Hopf  strip.  The  path  will  lie  in  a  region  bounded  by  the  singu¬ 
larity  at  k^  *  -k^  cos(8)  which  represents  the  incoming  wave  in  the  upper 
half  space.  Therefore  we  can  take  this  path  to  lie  along  the  real  k^ 
axis  for  0  <  0  <tt / 2  using  the  aperture  formulation.  For  tt / 2  <  0  <  tt  ,  we 
can  take  the  path  to  lie  on  a  line  with  the  imaginary  value  of  k^  given 
by  T  *  Im(-k^  cos(8))  and  indented  counterclockwise  below  the  pole.  For 
this  case  we  can  use  the  expression 

Im(-k^  cos  0)  *  k^cosS  *  Im(k  sin  a) 

*  k  cos  u  sinh  v  -  k.  sin  u  cosh  v.  (A. 6.1) 

r  i 

However  this  is  not  convenient.  The  main  significance  of  this  curve  is 
the  following  observation.  If  we  consider  the  integration  path  to  be  the 
real  axis  in  the  ot  plane  for  the  sake  of  simplicity,  then  as  the  geometri¬ 
cal  optics  pole  moves  to  the  right  of  the  real  kx  mapping  (9  >  ^ / 2 ),  it 
will  only  be  crossed  in  the  deformation  of  the  original  integration  path 
when  the  steepest  descent  path  is  to  the  left  of  the  geometrical  optics 
pole.  This  point  is  made  visually  in  Figure  A.l.  Put  another  way,  the 
only  possible  singularity  of  the  integrand  which  can  lie  between  the  true 
integration  path  and  the  real  k^  axis  mapping  is  the  geometrical  optics 
pole . 

A. 7  Locus  of  the  Geometrical  Optics  Pole 


A  singularity  whose  residue  contributes  to  the  far  field  and  is  of 
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particular  interest  is  the  root  of  (k  sin  a+  k^cosS),  where  k  is  the  wave 
number  of  the  observation  medium,  k^  is  the  wave  number  of  the  upper 
medium  and  6  is  the  angle  of  the  incident  plane  wave  (0  <  9  <tt  ) .  This 
contribution  constitutes  the  so-called  geometrical  optics  field.  Clearly, 
if  k  *  (that  is,  the  observation  point  is  in  the  upper  medium),  then 
ar  *  0  -  it/ 2  and  the  locus  of  the  geometrical  optics  pole  (GOP)  in  the 

a  plane  is  the  real  a  axis  between  u  *  — (tt / 2 )  and  u  *  +  (tr/2).  If  k  * 
(lower  medium  wave  number),  then  the  locus  of  the  GOP  leaves  the  real  a 
axis,  crossing  it  at  a  *  0. 

To  compute  the  locus  of  Che  GOP  for  k  *  k2 ,  we  expand  the  root  ar 
in  the  a  plane 

k2sin(ar)  *  -k^cosCS).  (A. 7.1) 

Substitute  k^  ■  klr-jkli*  k2  *  k2r  -jk2i  and  a  =>  u+jv  and  after  some  alge¬ 
braic  manipulations  arrive  at 


tan(u) 


klrk2r  +  klik2i 


klrk2i  _  klik2r 


tanh(v)  for  k^rk2i  ^  klik2r' 


Note  that  if  ^r*^2i  *  kli'k2r‘  tlTen  ^1  ^es  011  line  in  the  k^  plane 
which  runs  through  both  k2  and  the  origin.  In  this  case  the  GOP  locus 
will  run  between  -  ci^  ^  and  +  ^  along  the  line  segments  shown  in 

Figure  A. 2. 

Clearly  for  a  specific  value  9  we  can  use  the  standard  mapping  for¬ 


mula  (A. 1.2)  to  find  the  location  of  the  GOP  and  it  will  lie  on  the  curve 
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specified  by  (A. 7.2).  We  note  in  passing  that  it  is 


to  specify  a  curve  which  moves  from  oi  onto  the 

V 


possible  for  (A. 7. 2) 
lower  sheet  of  the  a 


plane  and  running  to  ot  =  -  tt  rather  than  to  a  *  0  as  one  might  expect. 

This  corresponds  to  the  situation  in  the  k  plane  where  the  locus  of  the 

x 

GOP  (a  straight  line  from  -k^  to  +k^)  crosses  the  branch  cut.  In  ef¬ 
fect,  a  pole  moving  from  k  =  -k,  to  k  =  0  is  seen  to  move  onto  the  lower 
sheet  projection  and  simultaneously,  the  root  from  the  lower  sheet  is  seen 
to  emerge  onto  the  upper  sheet  projection. 

It  is  useful  to  establish  the  slope  and  curvature  characteristics 
of  the  GOP  locus  in  the  ct  plane  (only  the  projection  of  the  second  quad¬ 
rant  of  the  kx  plane  is  of  interest).  Before  proceeding  it  is  necessary 
to  consider  a  classification  scheme  for  the  second  quadrant  of  the  kx 
plane  (and  its  image  in  the  a  plane)  which  is  based  on  the  properties  of 
the  mapping  k^  *  k^sina.  The  following  equivalences  hold  in  all  cases: 


kx  plane 


a  plane 


1.  The  straight  line  segment 
between  =  0  and  kjj  *  -k2  . 


1.  The  segment  of  the  real 
axis  between  u  =  0  and 
u  -  -(tt/2). 


2.  The  extension  of  the  above 
line  segment  from  kx  *  -1^ 
to  00 . 


2.  The  line  segment  for 
u  *  —  (tr /2)  and  v  <  0. 


3.  The  straight  line  segment 
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for  (k  )  *  -k.  and 

x  r  2r 

(k  ) .  >  k- . . 
x  i  —  2i 

4.  The  extension  of  the  line 
segment  in  3  for  (kx)^  <  k^^ 


3.  The  SDP  through 
a  =  -(tt/2). 

4.  The  SAP  through  a  =  -(tt/2) 
(only  the  portion  on  the  top 
sheet  is  of  interest). 


5.  The  extension  of  the 

branch  cut  hyperbola 

from  k  ■  -k  (asymptotic 
x  2 

to  the  negative  real 

k  axis), 
x 


The  locus  for  v  <  0 

represented  by  the 

mapping  of  the 

imaginary  k^  axis 

moved  to  intersect  u  =*  — ( rr/ 2 ) . 


These  images  are  shown  in  Figure  A. 3.  As  mentioned  these  particular 
loci  will  always  map  as  indicated.  Now  we  classify  the  second  quadrant 
of  the  k^  plane  and  its  image  in  the  a  plane  as  shown  in  Figure  A. 4. 

For  the  purposes  of  this  discussion  we  will  treat  regions  2A  and  2B  as  a 
combined  domain,  region  2. 

First  we  form  the  derivative  (dv/du)  of  the  GOP  locus  equation 
(A. 7. 2)  as 


dvl 

du) GOP 
locus 


cosh2 (v) 

K  cos2(u) 


(A.7.3) 


where 


< 


klrk2r 

+  kUk2i 

klrk2i 

"  klik2r 

|—  k  . 
r  _i±)  _ 

k.  -j 

f-  -ii) 

. 

J 

(  k,  ' 

(  k 

-  Klr 

2r-> 

Y/A  Region  of  no  interest 

(a)  (b) 

Figure  A. 3.  Prominent  lines  in  the  k^  plane  and 
their  projections  in  the  a  plane. 

a.  Lines  in  the  k^  plane. 

b.  Curves  in  the  a  plane. 


Figure  A. 4.  Regions  in  the  k  plane  and  their  projections  in 
the  a  plane. 

a.  Regions  in  the  k^  plane. 

b.  Regions  in  the  a  plane. 
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Clearly  <  >  0  if  the  line  in  the  plane  from  -k.,  to  the  origin  lies 


above  the  GOP  locus  (that  is,  <  — )  in  region  2);  otherwise  <  <  0. 
Forming  the  second  derivative. 


(  A  2  1 
d  v 

—*  — 

2(1+K2)crshu  (v) 

sin(u)cos(u) 

du2J 

GOP 

locus 

_  <2COS4(u)  j 

(r 

[  K 

(A. 7.4) 


For  in  region  1,  <  <  0  and  -(ff/2)  <  u  <  0  so  (d2v/du2)  _>  0.  For  a. 
in  region  2,  <  >  0  so  (d2v/du2)  ^  0.  For  in  region  3  we  see  from 
Figure  A. 4. a  that  the  GOP  line  will  cross  the  boundary  between  regions  1 
and  3.  Now  since  a^  is  in  region  3,  <  <  0;  but  for  the  region  1  portion 
of  the  GOP  locus,  -(tt/2)  <  u  <  0  while  for  the  region  3  portion  of  the 
GOP  locus,  -T  <  u  <  -(tt/2).  These  results  are  collected  in  table  A.  1 
for  the  GOP  locus  in  the  a  plane  and  shown  in  Figure  A. 3. 


For  a,  in 
b 

d v 
du 

d2v 

du2 

Region  1 

>0 

Region  2 

>0 

<0 

Region  3  (as  Reg  1) 

£° 

>0 

Region  3  (as  Reg  1) 

<  0 

_ 

<_Q 

Table  A.  1.  Curvature  characteristics 
of  the  Locus  of  the  Geometrical  Optics 
Pole. 
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A. 8  Locus  of  the  Connection  from  -k^  to  -k^ 

In  Section  3.4,  which  deals  with  the  factorization  of  the  function 
G21  we  saw  that  the  branch  cuts  for  and  S 2  can  be  deformed  along  the 
connection  lines  between  the  branch  points  and  will  cancel  on  the  semi¬ 
infinite  portion  of  the  cuts  resulting  in  finite  cuts  between  the  branch 
points.  (See  Figure  3.3).  It  is  also  demonstrated  in  Section  3.4  that 
the  function  (and  hence  G)  will  not  have  roots  on  the  top  sheet  using 
these  finite  branch  cuts.  It  is  important  to  the  discussion  in  Appendix 
C,  which  examines  the  exposure  of  these  roots  in  the  process  of  deforming 
the  S-^  branch  cut,  that  the  locus  of  the  connection  between  -k^  and  - k ^ 
in  the  a  plane  be  established. 

In  this  section  the  formula  for  the  connection  locus  from  -k^  to  -k2 
in  the  a  plane  is  given  for  observation  angles  in  the  lower  medium,  k£. 

A  similar  analysis  can  be  performed  when  the  observation  angle  lies  in 
k^.  We  start  with  the  parametric  description  of  the  connection  in  the  kx 
plane : 

kx  “  k2  +  (ki  -k^t,  for  t  e  [O.ll.  (A. 8.1) 

Substituting  kx  3  k2  sin(ci)  into  (A. 8.1),  the  expression  can  be  reduced 
to  the  equality  of  the  real  and  imaginary  parts. 

1  +  cos  (2u)cosh(2v)  =  -(<p-2)t  (A. 8. 2) 

-sin(2u) sinh(2v)  =  -ft  (A. 8. 3) 

where  a  *  u+jv  and  4>,  V  are  given  by 
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*2>  -  ttlr“2r  +  »uk21)  +  Mklrk21  -  kuk2r)  -  A  +  JB  (A.8.4) 


;(kl'lt2>1.  2(A2-B.;).  +j_iAB_  .  *  +  j 


(A. 8. 5) 


After  manipulations  the  quadratic  formula  can  be  applied  to  arrive  at 


sinh(2v) 

where  c  =  cos(2u)  and  s  ■ 


c  2-s 2 

4>-2 


[s  ±  C  /l- (c2- s2 )  ] 


(A. 8. 6) 


sin (2u) . 


It  is  necessary  to  establish  which  sign  is  appropriate  and  in  the  proc¬ 
ess  the  behavior  of  (A. 8. 6)  must  be  examined.  The  result  of  this  analy¬ 
sis  for  the  path  between  kx=  -k  ^  and  -k2  is 

.  ^  -  3g”^  "  !u+f!)  |.  ./rrrzrrzrr 


sinh(2v) 


c2_s2 


c 


£or  v(-kj)  *  0 


and  either  -  —  <_  u  _<  u^_k  ) 


or  u(_ki>  1  »  1  -  2 


where 


c  »  cos(2u) 


s  *  sin(2u) 


$,4*:  Defined  by  (A.8.4)  and  (A. 8. 5) 
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An  exceptional  case  for  (A. 8. 7)  occurs  if  '•>  =  0,  then  s  00  .  This  will 
occur  if 


klrk2i  *  klik2r 


kli  .  k2i 
klr  k2r 


(A. 8. 8) 


Hence  when  this  case  occurs,  -k^  and  lie  on  the  same  ray  from  the 

origin  in  the  plane.  If  k^r  <  k2r  then  we  have  that  v  =  0  on  the 
connection  path  in  the  a  plane  and  the  locus  lies  on  the  real  a  axis. 

If  k£r  <  k^r  then  the  connection  must  lie  on  the  line  u  =  -(tt/2),  0. 

We  now  examine  the  slope  and  curvature  characteristics  of  the 
connection  line  in  the  a  plane  using  the  same  classification  system  as 
in  Section  A. 7.  By  forming  the  derivatives  of  the  locus  equation  for 
the  connection  line  between  -kj_  and  -k2  and  manipulating  the  results,  it 
can  be  shown  that  the  curvature  of  the  locus  is  as  given  in  table  A. 2. 


For  a  in 

d2u 

dv7 

Region  1 

>  0 

n 

Region  2 

Region  3 

>0 

B 

Tablr  A.  2.  Curvature  Characteristics  of  Connection 
from  -k^  to  "^2” 
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These  results  imply  that  if  (kx  =  -k^)  is  in  region  1  or  2,  then 

the  entire  locus  is  in  that  region.  It  is  possible  for  to  be  in 
region  3  to  have  the  locus  cross  the  SDP  boundary  and  its  image  appear 
in  region  1.  This  case  is  considered  in  more  detail  in  Section  C.4. 

In  concluding  this  section  on  the  connection  line  between  -k^  and 
-k^  »  it  is  interesting  to  note  that  in  the  kx  plane  the  curve  is  actu¬ 
ally  a  portion  of  a  hyperbola  which  passes  through  kx  =  -k^  and  -k£  and 

has  asymptotes  which  are  perpendicular.  This  is  due  to  the  fact  that 

.  .  .  .  2 
the  connection  line  lies  on  a  straight  line  in  the  kx  plane.  This  is  a 

.  2 

property  of  the  mapping  kx  and  accounts  for  the  fact  that  the  hyperbolas 
used  to  define  the  original  branch  cuts  for  and  &2  will  map  into 
parallel  straight  lines  in  the  kx  plane. 


APPENDIX  B 


ASYMPTOTIC  FAR  FIELD  RESULTS 
FOR  THE  INTERFACE  PROBLEM 


B.l  Introduction 

The  purpose  of  this  appendix  is  to  present  explicitly  the  asympto¬ 
tic  contributions  to  the  far  fields  in  the  interface  problem.  The  solu¬ 
tion  in  the  transmitted  region  (y  <  0)  is  derived  in  this  appendix.  Simi¬ 
lar  steps  can  be  followed  in  the  reflected  region  and  the  generalized  re¬ 
sults  are  given  in  Chapter  4.  Throughout  this  appendix  we  follow  the  ap¬ 
proach  of  Felsen  and  Marcuvitz  (1973a). 

B.2  Asymptotic  Far  Field  Constituents  in  the  Transmission  Region 


The  integral  to  be  evaluated  is  shown  below: 


ez(x,y')  -5 


k^cosa 


jk2psin(a-4>) 


G  (k2sina) (k^sina  +  k^cosS) 


da 


(B.2.1) 


for  y  <  0 


where  £  *  constant  » 


-jEQ  sin0 

irri1G+(-k1cose)  ’ 


Eg  *  magnitude  of  incident  field 

9  *  angle  of  propagation  of  incident 

plane  wave 

<J>  ■  observation  angle  (-ir«fc<0) 


G+,G_ 


factorizations 
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p  *  radial  distance  from  the  edge  of 
the  PEC  screen,  normalized 

*  normalized  values 

T  ■  path  of  integration  in  the  a  plane 

The  singularities,  which  may  be  crossed  (thereby  contributing  explicitly 
to  the  asymptotic  field),  are  as  follows: 

1.  Geometrical  optics  pole  due  to  the  root  of 
(k^  sin  a  +  k1  cos  0). 

2.  A  pole  due  to  the  root  of  G_  in  the  mapping  of 
the  second  quadrant  of  the  kx  plane  (or  on  the  adja¬ 
cent  lower  sheet ) , 

3.  Branch  point  of  *  /k|  -k£  in  the  mapping  of 
the  second  quadrant  of  the  k ^  plane.  Due  to 

G_ (k^sin  a) . 

Each  of  these  contributors  is  considered  below. 

B.3  Steepest  Descent  Integral 

After  having  deformed  the  integration  path  of  the  radiation  inte¬ 
gral  (B.2.1)  to  the  steepest  descent  path  through  the  saddle  point, 

/2 ,  one  is  in  the  position  to  evaluate  an  asymptotic  approxima¬ 
tion  to  the  saddle  point  integral. 

The  first  two  terms  of  the  asymptotic  expression  for  the  steepest 
descent  integral  are  known  to  be  (Felsen  and  Marcuvitz  (1973a)) 


i 


^ez^saddle 

point 


1 


) 

r 

f 

2  *  J 

!  f  (v  i 

f(a  )  d2f (a) 

S  l  S 

•i  k 

J  k2 

P*5  j  2  k2  p3/2 

l  4  d°2 

1 

'i 

X 


for  <p  <  0  , 
where 


e“Jk2p 


(B.3.1) 


f  (a) 


E,  cos  a 

G_(k2  sin  a)  (k^  sin  a  +  k^  cos  9) 


(B.3.2) 


as  *  tr/2  +  <J>  and  p  is  the  normalized  distance  from  the  edge  of  the 
PEC  half-plane  (P  *  2^  distance  in  wavelengths  in  free  space).  As  long 
as  <f>  is  not  in  the  neighborhood  of  <P  ■  0°  or  -180°,  the  lead  term  in  the 
series  will  dominate  for  observation  distances  greater  than  a  wavelength 
in  free  space  from  the  edge  of  the  PEC  half-plane.  As  the  observation 
point  approaches  the  interface  between  the  media,  however,  the  (cosa) 
factor  in  (B.3.2)  causes  the  lead  term  in  (B.3.1)  to  vanish  and  the 
second  term  dominates  the  expansion.  In  principle  this  evaluation  can 
be  performed  but  this  avenue  of  work  has  not  been  pursued  since  the  ex¬ 
pressions  are  tedious  and  of  little  practical  value.  Instead,  the  fol¬ 
lowing  observations  are  made  concerning  the  significance  of  the  second 
term  in  the  asymptotic  expansion. 

For  observation  radii  exceeding  one  wavelength  in  free  space,  the 

lead  term  dominates  the  first  part  of  the  second  term  in  (B.3.1).  It  is 

appropriate,  therefore,  to  compare  the  behavior  of  the  lead  term  with 

the  second  part  of  the  second  terra  in  (B.3.1).  For  <t>  ^  -  180°,  the 

second  derivative  of  f ( ot  )  behaves  as 

s 
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„  -  2*u  sin  a 

f"  (a  )  % - i - c - § - 

s  Uo*8,  •  (k-sin  a  )4G  (k_sin  a  )*(k„sin  at  +  k,  cos  0) 
2  12  s-2  s2  si 


(B. 3.3) 


If  the  modulus  of  the  ratio  of  the  lead  term  to  the  second  term  in 
(B.3.1)  is  evaluated  in  the  vicinity  of  ct  *  -  tt/2  (that  is,  for 
<p  ■  -  180°)  using  (B.3.3),  we  see  that 


first  term  uf 

(B. 3.1) 

p 

U  *8. (k  sin  a  ) 

2  12  s 

second  term  of 

(B.3.1) 

tan  a  *sin  a 
s  s 

ui 

(B. 3.4) 

As  long  as  the  contrast  of  the  two  media  is  substantial,  we  see  that  an 
estimate  of  the  relative  significance  of  the  two  terms  is  given  by 

(B.3.5) 

cos  <p 

We  observe  that  if  p  is  taken  as  one  wavelength  in  free  space,  then  R  *  1 
for  $  *  171°  and  if  p  is  taken  as  five  wavelengths,  then  R  =  1  for 
<p  »  178°. 

In  effect,  we  assert  that  the  "propagation  field"  (which  behaves 
as  P  v  dominates  the  saddle  point  field  contribution  except  in  the  vici¬ 
nity  of  the  interface  plane  where  a  lateral  wave  type  of  field  exists  and 
preserves  the  phase  and  amplitude  propagation  behavior  of  the  total  field 
across  the  interface  aperture.  The  lead  term  in  (B.3.1)  is  given  by 
(4.1.8)  when  generalized  to  observation  points  in  either  medium. 

B.4  Geometrical  Optics  Pole 

This  pole  occurs  at  a  such  that  k»  sina  +  k.  cos0  =  0. 
r  gc  2  go  1 
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Since  the  integration  path  T  runs  from  v  =  -»  to  v  =  +  *>(where  ^  *  u+jv), 

we  see  that  the  pole  of  the  integrand  at  a  will  be  encircled  either 

go 

positively  or  negatively  when  crossed  depending  on  whether  the  geometrical 
optics  pole  lies  to  the  left  or  right  of  the  path  of  integration.  All 
that  remains  is  to  evaluate  the  residue  of  the  pole, 


Residue (a  ) 

go 


(a-a  ) 


-jk2psin(a-<t>) 
(£k2cosa  e  ) 


go '  G  (k2sina) (k2sina  +  k^cosS) 


-Ja«a 


go 


-jk  psin(ct~4>) 

(e  )  -  j  2tt  Res  (a  )  -  ^  e  Z 

z  go  J  go  G  (-k^cosS) 


(B.4.1) 


2E„sin3  -jk  psin(a  -<t>) 

,  x  „  0 _  2  go 

go  niG_(-k1cos6) 


(B.4.2) 


B.5  Root  of  G2(k2sin  a) 

It  is  noted  in  Section  3.4  that  G+(a)  *  G_(-a),  G(a)  *  G+(a)*G_(a) 
and  that  G+  and  G_  each  have  a  single  root  which  must  be  located  at  the 
same  position  as  the  roots  of  G.  It  was  shown  that  G  does  not  have 
roots  on  the  top  sheet  of  the  analytic  continuation  of  G2  which  has  the 
finite  branch  cuts  connecting  the  branch  points.  (Here  we  refer  to  the 
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partition  G q_  of  G  given  in  (3.4,2)).  However  when  the  branch  cut  for 
3^  is  deformed  in  the  process  of  evaluating  the  asymptotic  contribution 
for  the  branch  cut  integral,  this  root  may  be  brought  up  to  the  top  sheet 
and  be  crossed  as  the  integration  path  deforms  to  the  SDP.  In  that  case 
this  pole  contribution  must  be  included  and  is  discussed  in  this  section. 
The  question  of  when  the  root  is  exposed  and  crossed  is  discussed  in  Ap¬ 
pendix  C.  The  residue  at  the  root  of  G_  is  given  as 
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A 
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Residue  (a  ) 
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(a-a  ) 


C  k^cosa  e 


-jk2psin(a-<)>) 


r  G_(k2sina) (k2sina  +  k^cosG) 


Ja»a 


C  k  cosa  e 
2  r 


-jk2psin(ar-<J>) 


k  sina  +  k. cos9 
2  r  1 


Lim 

a-+a 


(a-oO 


G_(k2sina) 


Y  •  Lim 
a-*a 


(a-ar) 


G_Tk2sina) 


F  (a-ar) ' 

Y  •  G  (k,sina  )  Lim  -<7r - — r 

+  2  r  w  |J(k2slna) 
r  *- 


y  •  G  (k.sina  ) 

•h  2 _ r_ _ 

G'  (k2sina^)  (k^osa^) 


(B.5.1) 


Now  we  observe 
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I 


G'  = 


dk 


1 

r_  — 

S!  +  V 

k  1  1  +  _L_ 

U2 

L 

(B.5.2) 


However  at  a  *  ar,  G(ar)  *  0  implies 


S1  *  “(62;I1)/U2  at  a=Ctr 


CB.5.3) 


Substituting  this  into  (B.5.2)  yields 


G' (k2sinar) 


•k  • 

X 

i 

.  .  1 

1 

r62yi 

U2B2 

W  y2 

- 

, 

k  =k„sina 
x  2  r 


y2  -  y2 

P2  P1 


VlU  2 


tana 


(B.5.4) 


B.6  Branch  Cut  of  8^  for  =  -k^ 

The  last  singularity  contributor  that  we  need  to  evaluate  is  the 

branch  cut  integral  due  to  the  presence  of  G_(kx)  in  the  integrand  of 

(B.2.1).  As  we  deform  from  the  integration  path  T  to  the  SDP,  there 

will  be  a  critical  angle  $  such  that  the  SDP  will  cross  the  branch 

cr 

point  for  all  <p  such  that  |<J>|  _>  |$  |.  In  order  to  deal  with  the 

contribution  of  the  branch  point  in  a  systematic  manner  which  will  be 
amenable  to  asymptotic  evaluation,  we  choose  to  follow  an  approach  which 
is  analogous  to  conventional  steepest  descent  analysis. 


| 


A 


1. _ 
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The  essential  point  in  steepest  descent  analysis  is  to  deform  the 
integration  contour  to  lie  along  a  path  of  constant  phase  in  the  expo¬ 
nential  factor.  As  a  result,  we  are  guaranteed  that  the  magnitude  of  the 
exponential  factor  will  change  most  rapidly  along  such  a  path.  Because 
this  path  is  chosen  to  pass  through  the  saddle  point  of  the  exponential, 
there  are,  in  fact ,  four  possible  paths  along  which  we  may  proceed.  Two 
of  the  paths  are  colinear  at  the  saddle  point  and  are  referred  to  as  the 
steepest  descent  path  (SDP)  while  the  remaining  two  colinear  paths  make 
up  the  steepest  ascent  path  (SAP).  However  if  we  choose  some  other 
point  in  the  complex  plane  and  examine  the  contours  of  constant  phase  of 
the  exponential,  we  find  that  there  is  only  one  set  of  colinear  paths 
from  the  chosen  point  (see  Figure  B.l).  One  direction  will  follow  a 
path  which  is  asymptotic  to  the  SDP  while  the  other  direction  will  cause 
one  to  move  on  a  path  which  is  asymptotic  to  the  SAP. 

Consider  the  following  procedure.  We  specify  a  given  observation 
angle  <p  such  that  1 4>  |  ^  |i>cr|  and  therefore  the  integration  path  must  cap¬ 
ture  the  branch  point  which  characterizes  -  -k^  as  it  deforms  to  the 
SDP.  We  deform  the  branch  cut  for  to  the  steepest  descent  path  through 
ab  (SDP^)  which  is  the  portion  of  the  constant  phase  contour  that  is  asymp¬ 
totic  to  the  SDP.  Then  we  deform  the  integration  path  to  the  SDP.  In 
the  process  of  this  deformation  we  must  enclose  the  entire  branch  cut 
and  the  integration  paths  for  this  branch  cut  integral  will  lie  along  the 
steepest  descent  path  from  3^.  The  specifics  of  this  deformation  are  con¬ 
sidered  in  Appendix  C.  For  our  purposes  here,  we  assume  that  such  a  de¬ 
formation  has  occurred  and  the  branch  cut  for  8,  lies  along  the  SDP  from 

1  D 

ab  * 


We  are  now  interested  in  evaluating  the  branch  cut  integral  and  fol- 


Figure  B.l. 


An  example  of  the  steepest 
loci  for  an  arbitrary  point 


descent/ascent 

V 


AO  A 1 38  968  SCA1IERING  OF  AN  ELF.CI ROMAGNF  1 1C  PLANE  WAVE  FROM  A 

PERFECTLY  ELECTRICALLY ..( U )  MISSISSIPPI  UN1V  UNIVERSITY 
DEPT  OE  ELECTRICAL  ENGINEERING  R  D  COBLIN  ET  AL . 
UNCLASSIFIED  FEB  83  N00014-81 -K-0256  F/G  20/14 


184 


low  the  analysis  discussed  in  Section  4.8  in  Felsen  and  Marcu.vitz  (1973a). 
We  want  to  characterize  the  portion  of  the  integrand  in  (B.2.1)  exclusive 
of  the  exponential  in  the  vicinity  of  o^.  Hence, 


5  l^00801 


*  G  (l^sinct)  •  (k2sino  +  k^cos0) 
jl;  *k2  *  cosa*  G+ (k2  s  inct) 


k2Si.n0.  +  k^cos0 


u(a) 


G(k2sino) 


k^cosa 


("A,  -  k_  sina 


u(a) 


y^tl+^2Sina 


k2Cosa 


k2cosa 

2 

k^-k^sin2 (of 

h  -1 

/k,  -  k„sina 


A.i+k2sina 


(B.6.1) 


For  a  near  o^,  we  have 

n 

[f(oO] 


4  near 


1  U^C-k^" 

-  £  ^2  u2  G+(-k1) 

k^cosS-l) 

■X 

Ul  ^nqi^tcose-D 

/(a-o^) 


(B.6.2) 


where  /  refers  to  the  analytic  continuation  of  the  original  branch 
of  viTJ+iT^TinO  in  the  vicinity  of  o^Cdiscussed  below).  /-*v 

is  the 

principal  value  with  the  branch  cut  along  the  negative  real  axis. 
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Before  proceeding  we  need  to  clarify  the  meaning  of  /a  -  .  To 

do  this  we  begin  with  /k  ^+k  2  sin  a.  The  original  specification  of  this 
function  is  given  in  Section  2.2  (  specifically  (2.2.2)  and  Figure  2.3). 
We  see  there  that  the  quantity  (  kj+k^  *  kx-(-k^)  can  consi-dered 
as  a  vector  from  (-k^)  to  kx>  If  this  is  done  then  the  following 
equations  hold  for  the  branch  definition: 


If  (kj*  kx)  *  [kj*  kje-59  .  then  /(k^  kx)  -  Ik^  kj**  e^2. 


(B.6.3) 


It  is  shown  in  Section  C.3  that  the  SDP^  branch  cut  for  never  crosses 

the  a  plane  locus  of  the  straight  line  from  k  ■  -k,  to  k  ■  0  in  the  k 

XXX  X 

Q 

plane  (locus  of  the  geometrical  optics  pole).  Hence  the  function  ^kj+k 
will  retain  its  original  value  from  the  hyperbolic  branch  cut  definition 
for  all  points  on  this  line.  If  we  examine  points  on  this  line  (referred 
to  as  a  )  in  the  a  plane  near  a^,  we  find  that 


(-k2)t  -  (k^+k^sin  a^)  for  te[0,l] 

~  "ki  +  vv^w 

[S2(k1)]*(a  -ab)-(l-t)k1  . 


(B.6.4) 


We  have : 


Arg  02 (k^)  +  Arg  (a  ~a^)  -  Arg  ^ 


Arg(/k1+kx0)  -  4  Arg  k^Argf  /B^k^ 


l2(kj_)  •  -oth)9]  -  H  Arg  B2(k. 


)  +  h  Arg  (a  "V 


Hence 
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So 


ArgCa^-ct^)  -  0  -  Arg  ^  -  Arg  B2(k1) 


and 


^S2(k1)*(ag"ab)9  "  /S2(kl)I>V  *  IV^b^  e  8 


JV2 


7—  Q 

Having  specified  values  on  the  branch  of  /  near  c^,  we  can  analytical¬ 
ly  continue  them.  So, 

/02(k1)*(a-ocb)0  -  /S^k^  •  la-c^l*5  e^2'  (B.6.6) 

where  a  is  near  and  0  assumes  a  value  by  rotation  from  0^.  Consider,  • 
in  particular,  the  ray  from  ct^ along  the  SDP^  which  results  from  the  coun¬ 
terclockwise  rotation  from  9^  and  call  the  angle  9+.  Then 

f(a+)  ■  y  +  wv'da-o^;  .  (B.6.7) 

where  Y  and  m  are  given  in  (B.6.2)  and  ^a+  -  »  |a+  -  e  + 

with  0^  <  0+  <  9g+2TT  and  0+  is  the  angle  of  the  SDP^  from  a^.  Turning  to 
the  analysis  of  Felsen  and  Marcuvitz,  we  see  that  the  branch  cut  contri¬ 
bution  for  a  branch  integral,  which  has  the  integration  path  running  coun¬ 
terclockwise  around  the  branch  point,  is  given  by 


mJj  mWWZ> 

IkjPco.Co^-*) |3^2 


CB -6.8) 


Generalizing  the  result  for  the  reflection  half  space  gives  (4.1.16). 


APPENDIX  C 


DEFORMATION  OF 
THE  BRANCH  CUT 

C.I  Introduction 

The  case  to  be  considered  in  this  appendix  is  that  of  the  interface 
problem  where  the  field  integral  is  to  be  evaluated  in  the  transmitted 
region,  i.e.  the  same  case  as  that  in  Appendix  B  ( specif ical ly(B. 2 . 1 )) . 

As  discussed  in  Section  B.6,  it  is  useful  to  deform  the  branch  cut  of  8^ 
from  kx  *  -k^  arising  in  the  factor  G_(k2sin  a)  in  (B.2.1).  In  partic¬ 
ular,  we  deform  to  the  steepest  descent  path  (SDPfe)  from  the  point 
where  **k^*  k^inc^*  Section  C.2  examines  the  constant  phase  contours 
specifying  SDP,  for  a  given  saddle  point  location  (a  )  and  develops  the 
relevant  equations.  In  Section  C.3,  we  examine  the  effect  of  the  deform¬ 
ing  branch  cut  integral  on  the  various  singularity  contributors.  Final¬ 
ly,  in  Section  C.4  we  consider  the  exposure  of  the  root  of  G_Ck2sinCt). 

C.2  Constant  Phase  Contours 

The  steepest  descent  path  (and  steepest  ascent  path)  from  the  saddle 
point  c*8  is  actually  a  path  in  the  a  plane  along  which  the  phase  of  the 
exponential  in  (B.2.1)  is  constant.  Therefore 

Iaf-Jkj  3in(a-$)]  -  Im[-jk2cos(a-aa)]  .  (C.2.1) 
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where  ot  »  u  +j  v,  k2  ■  k2r-j  k2i’  k2r*  k2i  —  °* 

So  Che  phase  is  given  by 

Y  ■  k2^sin(u-us)sinh(v)-k2rcos (u-ug)cosh(v) . 

Along  the  SDP  (with  respect  to  ctg),  y  *  ~k2r'  c*ear  the 

discussion  in  Section  4.1.3  (and  Figure  4.4),  that  if  ag  ■  0,  then  the 
region  between  the  two  shifted  mappings  of  the  imaginary  k^  axis  (the 
zero  phase  lines)  will  have  Y  <  0.  It  is  useful  to  consider  the 
contours  of  constant  phase  which  lie  to  the  right  of  the  SDP  (  and  hence 
characterize  points  which  will  be  crossed  by  the  SDP  for  -it  <  <p  <  -tt/2). 
See  Figure  C.l.  Observe  that  for  -k  <  Y  0,  the  contours  run  more  or 
less  vertically  and  lie  between  the  SDP/SAP  and  the  right  zero  phase 
line.  For  Y  <  “k2r  tlie  contours  curve  downward  and  lie  between  the  SDP 
and  SAP.  We  observe  another  point:  namely  it  is  clear  that  the  portion 
of  the  contour  which  is  asymptotic  to  the  SDP  is  the  SDP^  from  the  point 
ot^.  Hence  if  Y  >  -k2r’  t*ie  ^**b  runs  UP  an<*  to  v  *  +°°  while  if  Y  <  ~k2r’ 
the  SDP^  curves  to  the  left  and  down  to  v  ■  -  ®.  Therefore  as  the 
saddle  point  moves  to  -tt/2,  the  SDP^  deforms  continuously  as  shown  in 
Figure  C.2.  We  notice  in  particular  that  when  lies  on  the  SAP,  the 
SDP^  may  be  taken  as  running  to  either  v  *  -  00  or  v  *  +  00 .  We  can  view 
the  phase  contour  plot  of  Figure  C.l  as  fixed  to  otg  and  visualize  as 
moving  to  the  right  on  the  phase  plot  as  a  moves  from  0  to  -tt/2. 

We  have  seen  in  the  above  discussion  that  the  SDP^  alters  radically 
as  ctg  moves  to  -tt/2.  However,  note  that  the  SDP^  will  always  lie  to  the 
left  of  the  SDP  through  cig  which  passes  through  a^.  This  can  be  seen  by 
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observing  that  if  the  derivative  (dv/du)  is  formed  from  (C.2.1),  then  in 

all  cases  the  SDP^  will  have  a  tangent  from  which  moves  into  the 

region  to  the  left  of  the  SDP  through  and  remains  in  that  region. 

Also  we  know  that  the  SDP^  is  bounded  on  the  left  by  the  SDP  through  as. 

So  for  any  and  ag  we  have  a  band  between  the  two  SDP  loci  in  which  a 

root  of  G  (k^sina)  may  be  exposed  (and  obviously  crossed  by  the  SDP 

through  oi  ).  Viewed  the  other  way  if  the  root  of  G  lies  outside  this 
s  - 

band  then  it  cannot  lie  on  the  top  sheet  and  be  crossed  by  the  SDP  so  it 

can  be  ignored.  To  use  this  information  it  is  helpful  to  know  the  value 

of  a  (referred  to  as  a  ,  )  such  that  the  SDP  runs  through  a,  .  To  deter- 
s  sb  b 

mine  this  we  set  Y  ■  -k„  and  a  ■  a,  in  (C.2.1)  and  solve  for  a  ,  .  The 

2r  b  sb 

result  is 


cos(u  J  —  [k.  A±/B*(A*+B*-kO]  for  v  >  0  .  (C.2.2) 

sb  a2+b*  2t  2t 

where  a*  u+jv, 

A  -  Re [^(c^)]  -  k2j.cos(ub)cosh(vb)  -  k2isin(ub)sinh(vb) 

B  -  Relk^c^)]  •  k2rsin(ub)cosh(vb>  +  k2icos(ub)sinh(vb) 
and  -  ■=■  <  u  <  0  . 


Now  consider  an  arbitrary  point  ab(which  must  lie  in  the  second 
quadrant  of  the  kx  plane)  and  an  arbitrary  saddle  point  ag(between  — tt/ 2 
and  0).  We  wish  to  derive  an  explicit  equation  for  points  on  the  SDPb> 
To  begin  we  solve  (C.2.1)  for  a.  Then  we  invert  (C.2.1)  for  to  arrive 


at  the  following  expressions: 


If  Y  >  -k2r: 


sin(u-u  ) 


3  A2+B2 


-B  +  A/A2+B2-l 


for  (u-u  )  >0,  all  v,(C.2.3) 
s 


where 


a  •  u  +  j  v  and  A 


k^cosh (v) 


W1 


B 


k2±sinh(v) 


If  y±-k2r: 


sinh(v)  -  — - —  [B-sgn§  ~  u  +  u)  /A* (Bz-Az+1) ] 


where 


A2+B2 


k2r*cos(u) 


k2i*sin(u) 


for  -0^  <  (u-us>  <  8^ ,  v  <  0, 


9^  *  tan 


-1 


2r 


lk2ij 


(C.2.4) 

(C.2.4.a) 

(C.2.4.b) 

(C.2.4.c) 


on 


Note  that  in  the  limit  of  Y  ■  -k2r,  these  equations  describe  points 

the  SDP  and  SAP  which  define  the  bounds  of  their  respective  regions 
> 

(Y  <  ~k2r).  However  we  recognize  the  following  limitations  for  the  appl: 

cation  of  (C.2.3)  and  (C.2.4)  when  Y  •  -k2r: 

(C.2.3)  gives  u  as  a  function  of  v  for:  a.  SDP  for  v  >  0 

b.  SAP  for  v  <  0 
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(C.2.4)  gives  v  as  a  function  of  u  for:  a.  SAP  for-  u  >  ug 

b.  SDP  for  u  <  ug 

These  results  can  easily  be  extended  by  recognizing  that  the  SDP  and  SAP 
are  symmetric  with  respect  to  the  saddle  point  ct  . 

C.3  Effect  of  Deformation  on  the  Contributing  Singularities 

As  discussed  in  Section  B.6,  it  is  desirable  to  deform  the  8^  branch 

cut  from  k  ■  -kp  which  arises  from  G_(kx)  in  the  denominator  of  the 

field  integral,  to  lie  along  the  SDP  in  preparation  for  the  deformation 

b 

of  the  path  of  integration  to  the  SDP  through  the  saddle  point.  This 
permits  a  formulation  of  the  branch  cut  contribution  which  is  amenable 
to  asymptotic  analysis.  As  is  recognized  in  Sections  B.5  and  B.6,  this 
deformation  of  the  8^  branch  cut  has  two  effects  on  the  integrand 
resulting  from  the  fact  that  the  steepest  descent  analysis  is  performed 
on  an  analytic  continuation  of  the  integrand.  First,  the  question 
arises  concerning  value  of  8^(  that  is,  +  or  -  the  original  hyperbolic 
definition  of  8^).  Since  we  only  evaluate  the  integrand  at  the  positions 
of  the  singularities  in  the  complex  a  plane,  the  question  reduces  to 
determining  if  the  specific  singularity  is  crossed  in  the  process  of 
deforming  the  8^  branch  cut.  The  second  issue  arising  from  the  deforma¬ 
tion  is  the  determination  of  whether  a  root  of  G_  is  brought  onto  the 
top  sheet  of  the  integrand  (and  is  crossed  by  the  integration  path  as  we 
move  to  the  SDP  of  integration).  This  question  is  rather  involved  and 
is  considered  in  Section  C.4. 

We  begin  by  recalling  the  field  integral  to  be  evaluated 
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®+Jc 


j (6-y-k  *) 


-®+jc 


G_ • (k^+k^cosS)  e 


dk 


?JC__VV 

J  cog-tt^ 


j (B,y-k x) 

cose!  e  dkx  fory<0.  cc.3.1) 


-°°+j  c 


We  note  that  the  second  form  displays  explici  r  the  B^  and  B^  branch 
cuts  in  the  second  quadrant  of  the  k^  plane  t  v  _,h  the  function  G(kx>. 
The  branch  cut  for  B^  affects  only  the  portion  of  the  integrand  outside 
of  the  exponential,  hence  deformation  of  the  branch  cut  does  not 
affect  convergence  of  CC.3.1).  In  the  second  form  given  in  (C.3.1)  the 
value  of  G+(kx>  is  unaffected  by  the  deformation  of  the  branch  cut  from 
kx»  -k^  and  can  be  evaluated  using  the  factorization  integral  without 
alteration. 

We  utilize  the  classification  scheme  of  the  Ot  plane  given  in 
Section  A. 7  and  add  a  level  of  refinement.  We  now  break  region  2  into 
2A  and  2B  as  shown  in  Figure  A. 4. 

Having  stated  this  classification,  we  begin  by  considering  the 
saddle  point  contribution.  Consider  in  region  1.  It  is  clear  from 
the  phase  contour  plots  (Figure  C.l)  that  the  SDP^  will  have  v  0. 
However  it  is  also  clear  from  Figure  A. 4. a  that  the  original  hyperbolic 
branch  cut  for  in  region  1  will  lie  in  region  1  and  not  cross  the 
connection  from  -k^  to  the  origin.  Therefore  we  have  demonstrated  that 
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if  the  branch  point  is  in  region  1  then  the  deforming  branch  cut  3^ 

will  not  cross  the  real  a  axis  and  so  the  saddle  point  contribution  will 

use  the  original  (3  definition.  Now  consider  a,  in  region  3.  Similar 
1  b 

arguments  to  those  for  in  region  1  apply  and  the  conclusion  is  that, 
again,  for  in  region  3  the  original  value  of  8^  is  correct  for  the 
saddle  point  contribution. 

We  now  consider  in  region  2  and  must  consider  regions  2A  and  2B 

separately.  Consider  first  ct,  in  region  2B.  When  we  examine  Figure 

b 

A. 4. a,  it  is  clear  that  the  original  hyperbola  will  cross  the  ray  from 

the  origin  through  -k^  at  some  point  beyond  -k2 .  Therefore  in  the  a 

plane  it  must  cross  the  line  u  =  — ( tt / 2 )  for  some  point  v  <  0  and  curve 

down  into  region  3.  At  this  point  we  consider  the  SDP  for  a  j.n  region 

b  b 

2B.  Referring  back  to  Section  C.2,  if  we  consider  the  derivative 
(dv/du)  of  (C.2.1)  and  set  it  to  zero,  we  have 


const 


phase 


contour 


k,rtan(u-ug)  +  k^tanhCv) 
k2rtanh(v)  -  k2itan(u-us) 


=  0 


which  gives 


^2rtan ^U-Us^  ™  ~k2^tanh(v) 


(C. 3 . 2a) 


(C. 3.2b) 


We  note  that  this  equation  is  identical  to  (A. 2. 2)  with  the  exception 

that  the  intersection  with  the  real  a  axis  occurs  at  u  *  u  .  Hence  the 

s 

locus  of  zero  tangent  lines  to  the  constant  phase  contours  lies  on  the 

image  of  the  imaginary  kx  axis  mapping  shifted  to  u  =  ug  and  when 

u  *  — < tt/2 )  this  locus  can  be  shown  to  map  back  into  the  k  plane  as  the 
s  x 
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extension  of  the  hyperbolic  branch  cut  through  -k^  which  is  asymptotic 
to  the  negative  real  k^  axis.  This  implies  that  the  boundary  between 
regions  2A  and  2B  in  the  a  plane  (which  is  the  extension  of  the  hyper¬ 
bola  going  through  -k2  in  the  kx  plane)  is  the  locus  of  zero  slope  for 
all  the  phase  contours  with  y  <  “k2r  and  ug  ■  -(tt/2).  Hence  we  are  guar¬ 
anteed  that  the  SDP^  will  go  into  region  3  and  that  in  the  course  of  de¬ 
forming  to  it  we  will  not  cross  the  saddle  point.  So  for  a^in  region  2B 
we  also  have  that  the  original  8^  definition  is  correct. 

Finally  we  must  consider  in  region  2A,  a  more  difficult  case. 

We  note  that  the  original  hyperbolic  branch  cut  will  now  cross  the 
saddle  point  locus  (Figure  A. 4. a)  in  all  cases.  Specifically  we  will 
assume  that  the  original  8^  hyperbola  crosses  the  real  a  axis  at  a  *  u^ 

where  -(tt/2)  <  u  <0.  Consider  the  behavior  of  the  SDP,  as  u  moves 
cr  —  b  s 

from  u  *  0  to  u  *  -(tt/2).  We  recognize  that  gi^en  a  there  is  a  unique 

D 

value  u  *  u  specified  by  (C.2.2)  for  the  saddle  point  at  which  the  SDP 
sb 

through  the  saddle  point  u  ,  crosses  the  branch  point  ct  ,  it  is  clear 

sb  b 

that  for  u  >  max(u  ,u  ,  ),  the  saddle  point  moves  on  the  original  top 
s  cr  sb 

sheet  of  8^.  Now,  if  ug*  us^>  the  SDP^  must  lie  on  the  SDP  through 
and  run  to  v  ■  For  simplicity,  we  will  deform  the  8^  cut  to  this  lo¬ 

cus  initially,  recognizing  that  in  the  process  all  values  of  u  <  u^  on 
the  real  a  axis  enforce  the  lower  sheet  value  of  the  original  hyperbolic 

definition  of  ft  •  Now  as  u  proceeds  from  u  ■  0  to  u  ■  -(tt/2),  it  will 
1  s 

run  into  either  u  or  u  ,  .  If  it  encounters  u  first  there  is  no  ef- 
cr  sb  sb 

feet  on  g  since  the  SDP  must  curve  below  the  real  a  axis  until  the  SAP 
1  b 

crosses  a^.  Therefore  only  top  sheet  values  of  8^  will  be  achieved  until 


l 
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ug  reaches  either  uCf  or  the  value  of  ug  at  which  the  SAP  crosses  o^, 

referred  to  as  u .  Assume  that  u <  u ,  then  for  u„_  <  u„  <  u , 

sap  sap  cr  sap  s  cr 

the  lower  sheet  value  of  Bi  is  correct.  Since  for  all  u„  <  u„„  the  SDP, 

1  s  sap  b 

must  lie  to  the  right  of  the  saddle  point,  we  must  use  the  top  sheet  value 

of  as  the  deformation  from  the  original  B^  branch  cut  to  the  SDP^  will  not 

cross  Che  saddle  point.  Therefore  for  the  case  where  u  <  u  ,  we  have 

sap  cr 

the  following: 


f+(B, )  j  ,  u  <  u 
1  lorig  cr  —  s 


-((3,)  .  ,  u.,B  <  u  <u  where  -  -z  <  u  <  0  (C.3.3) 

1  orig  SAP  —  s  cr  2  —  s 


t+(6l>oris  '  us  <  USAT 


Now  we  assert  that  u^,.,^  <  u„„  for  all  cases.  This  can  be  seen  from 

sap  cr 

a  somewhat  heuristic  argument  based  on  Figure  A. 3.  Clearly  the  possible 
exceptional  points  must  occur  in  the  triangle  in  Figure  A. 3. a  bounded  by 
the  lines  1,  4  and  the  real  k„axis,  since  branch  points  outside  this 
region  (but  in  region  2A)  will  not  have  u  >  -(tt/2).  But  observe  that 
the  SAP  in  A.3.b  maintains  a  constant  acute  angle  from  the  saddle  point 
with  respect  to  line  1.  Since  the  mapping  (kx»  k2sin  a)  is  conformal 
except  at  a  ■  — C  it/ 2) ,  this  angle  is  maintained  in  the  k^  plane  with 
respect  to  line  1.  However  the  B]_  hyperbola  clearly  has  a  tangent  line 
at  the  intersection  with  line  1  which  maintains  a  greater  angle  (usually 
obtuse)  with  respect  to  line  1.  This  angle  must  be  maintained  in  the  a 
plane.  But  both  the  original  B^  hyperbola  and  the  SDP ^  must  originate 
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at  cl  ,  so  clearly  u  <  u 

sap  cr 

Summarizing  the  results  for  the  value  of  to  be  used  for  the 
saddle  point  contribution,  when  the  effect  of  the  deformation  of  the  8^ 
branch  cut  to  the  SDP^  is  included,  we  have: 

For  in  regions  1,2B  and  3:  Use  *  +  ^l^original* 

(C. 3.4) 

For  in  region  2A:  Use  6^*  (C.3.3). 

Now  we  consider  the  effect  of  the  B^  deformation  on  the  geometrical 
optics  pole  (GOP).  We  recognize  that  the  locus  of  the  pole  is  a  straight 
line  in  the  k^  plane  connecting  -k^  with  the  origin.  Consider  the  case 
where  is  in  region  1.  Examining  Figure  A. 4. a,  we  see  that  the 
hyperbolic  branch  cut  and  the  GOP  locus  intersect  only  at  the  branch 
point,  -k^.  Considering  Section  A. 7  and  examining  Figure  A. 5  (which  de¬ 
picts  the  curvature  of  the  locus  of  the  geometrical  optics  pole  in  the 

a  plane),  we  see  that  the  same  will  be  true  for  the  SDP  ,  which  must  run 

0 

up  from  to  v  "  +00.  Therefore  the  GOP  locus  will  not  be  crossed  during 
the  deformation  of  the  8^  branch  cut  to  the  SDP^  when  is  in  region  1 
and  hence  the  original  top  sheet  definition  of  8^  is  correct.  A  similar 
analysis  is  seen  to  apply  for  ct^  in  region  2.  Finally  we  consider  in 
region  3.  As  mentioned  above  all  GOP  locus  lines  will  cross  the  boundary 
line  between  regions  1  and  3.  Also  we  know  that  for  a  in  region  3  the 

D 

deformed  SDP^  cut  must  stay  in  region  3.  Regarding  the  original  hyper¬ 
bolic  branch  cut  for  8^,  we  see  from  Figure  A. 4. a  that  it  does  not  cross 
the  GOP  locus  and  in  deforming  to  the  SDP^(which  lies  totally  in  region 


3)  that  it  will  not  cross  the  locus.  Therefore  for  in  region  3,  the 
original  top  sheet  value  of  is  correct.  Collecting  these  results  for  the 
GOP  contribution,  we  see  that  the  original  value  for  8^  will  always  apply 
Having  established  these  rules  for  choosing  the  value  of  6^  when 
evaluating  the  asymptotic  contribution  of  each  singularity,  we  see  that 
they  are  amenable  to  implimentat ion  during  numerical  evaluation  of  the 
far  fields  in  a  problem. 

\  C.4  Effect  of  the  8.  Branch  Cut  Deformation  on  the  Hoot  of  G  (k  ) 

We  begin  by  determining  under  what  situations  the  root  of  G  (k^ ) , 

is  exposed  in  the  deformation  of  the  8.  branch  cut  to  the  SDP,  and  is 

x  b 

crossed  by  the  deformation  of  the  integration  path  to  the  SDP  through 
the  saddle  point.  We  utilize  the  same  division  of  the  a  plane  as  pre¬ 
sented  in  Section  A. 7  (that  is,  into  regions  1,2  and  3).  In  addition 
we  need  the  following  terminology: 

I 

a  ■  u+jv 

■  ub+jvb  is  the  branch  point  in  the  a  plane  such  that 

r  -kx-  ^sir^. 

a  or  u  is  the  value  of  the  saddle  point  ( — (tt /2 )<  u  <  0). 

ss  —  s — 

i  otg  *  Ug+jVg  is  the  root  of  G_(kx>. 

a  ,  or  u  ,  is  the  value  of  the  saddle  point  such  that  the  SDP 
sb  .  sb 

i  intersects  a^. 

i 

a  „  or  u  .  is  the  value  of  the  saddle  point  such  that  the  SDP 
sO  .  su 

intersects 

a  *  u  +jv  refers  to  points  on  the  connecting  line  between 
conn  conn  conn 

k^  -  -k^  and  -k£  (discussed  in  Section  A. 8). 


I 
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Y(ot,ug)  refers  to  the  value  of  the  phase  function  (C.2.1) 


We  recognize  that  the  branch  cut  for  6 ^  *-n  the  second  quadrant  of 
the  kx  plane  exists  in  both  the  exponential  factor  of  the  integrand  of 
the  field  integral  and  in  the  denominator  factor  G_(kx>,  while  the  B^ 
branch  cut  exists  only  in  G_(kx).  Therefore  we  can  view  the  B^  branch 
cut  as  deforming  from  the  original  hyperbolic  locus  to  overlay  the  hyper¬ 
bola  as  depicted  in  Figure  C.3.a.  In  the  a  plane  the  configuration  will 
appear  as  in  Figure  C.3.b  (shown  for  in  region  3).  We  will  take  this 
as  our  initial  6^  branch  cut  configuration  since  we  are  guaranteed  by  the 
analysis  of  Section  3.4  that  no  roots  of  G  will  occur  on  the  top  sheet. 


For  all  values  of  the  saddle  point  such  that  u  .<  u  ,  there  will  be  no 

sb  s 

roots  crossed  in  deforming  to  the  SDP.  We  know  that  the  root  of  G  lies 
somewhere  in  the  region  Re(kx)  _>  0  and  the  point  at  issue  is  whether  it 
is  exposed  on  the  top  sheet  and  crossed  in  deforming  the  integration  path. 
We  also  know  that  for  a  root  to  be  crossed  by  the  integration  path, 


— ( tt/2)  <  u  <  U  -. 

S  SO 

Consider  in  region  1.  When  ug  <  ug^,  then  the  SDP^  must  run  up 
to  v  ■  +  ».  Since  the  initial  branch  cut  for  8^  also  runs  to  v  ■  +  * 
along  the  SDP  through  a  ■  -(tt/2),  it  is  clear  that  any  root  with 
Vq  <_  0  will  never  be  brought  to  the  top  sheet  and  can  therefore  be  ig¬ 
nored.  We  assume  now  that  Vq  >  0.  We  now  observe  that  the  root  01q  must 


lie  to  the  left  of  the  SDP  through  and  o^. 


Since  all  SDP,  contours 

D 


lie  to  the  left  of  this  boundary  line,  any  root  which  may  be  exposed  in 


the  branch  cut  deformation  must  lie  to  the  left  of  this  curve.  Hence  we 


t 

*a« 

C.3.  Initial  deformation  of  the  branch  cut  for  0,  from 
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have  that  usQ<usb  i-n  order  for  a^to  be  exposed  at  all.  The  one  excep¬ 
tion  to  this  is  if  the  connection  line  from  -k^  to  -k2  of  the  initial, 
root-free  top  sheet  configuration  is  indented  below  the  SDP  through  a 
and  the  root  0Lq  lies  in  this  indention.  (Note  that  the  connection  line 
has  positive  curvature  and  positive  slope).  If  this  is  the  case  then 
for  all  ug  <  ugg  the  root  will  be  both  exposed  and  crossed  by  the  SDP 
through  ug.  See  Figure  C.4.  This  is  due  to  the  fact  that  for  ug  <  ugg, 
the  branch  cut  for  0^  which  initially  coincides  with  the  connection  line 
must  deform  to  coincide  with  the  SDP  (which  goes  through  a^) .  Further 
movement  of  ug  toward  u  *  -(tt/2)  will  force  the  branch  cut  to  expose  the 
root  which  is  clearly  crossed  by  the  integration  path  deforming  to  the 

SDP  through  u  since  u  <  u  .  Once  u  <  u  . ,  the  SDP,  deforms  up 
s  s  sO  s  sb  b  r 

toward  v  *  +<»  and  the  root  will  continue  to  be  exposed  as  ug  progresses 
to  -(tt/2). 

Continuing  the  analysis  for  in  region  1,  we  consider  that 
us0  <  usb  s0  *****  t*ie  root  °>ay  be  exposed.  If  ct^  lies  to  the  left  of 
the  curve  defined  by  the  connection  from  -k^  to  -k2  and  the  SDP^  running 
from  to  v  *  +®(See  Figure  C.5),  then  the  root  has  been  exposed  on 
the  top  sheet  (and  crossed  by  the  integration  path  since  ug  <  ugg). 
Otherwise  the  root  is  not  exposed  and  does  not  contribute. 

Consider  cx^  in  region  2.  We  again  note  that  unless  ugg  <  u^  the 
root  cannot  be  exposed  and  there  is  no  exception  this  time  as  the 
connection  line  has  positive  curvature  and  negative  slope.  See  Figure 
C.6.  Now  for  in  region  2  it  is  possible  for  the  SDP^  to  go  to 
v  *  ±  «  depending  on  whether  y^  ■  y(a^),  the  phase  function  for  the  phase 
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Figure  C.4.  Possible  region  of  exposed  roocs  of  G_ 
in  Che  a  plane. 


Figure  C.5.  Possible  region  of  exposed  roots  of  G 
in  the  a  plane. 
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Figure  C.6.  Possible  region  of  exposed  roots  of  G 
in  the  a  plane. 
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contour  through  a^,  is  ^  ”^2r*  Assume  first  that  Yfc  <  ”^2r  s0  that  the 
SDP^  curves  down  to  v  a  -«.  Clearly  if  Vq  0,  then  the  root  is  never 
exposed.  Assume  Vq  <  0.  Then  a  possible  situation  is  shown  in  Figure 
C.7.  If  aQ  lies  between  the  curves  defined  by  the  lower  leg  of  the  SDP 
through  a  -  -(tt/2),  the  connection  line  from  -(tt/2)  to  and  the  SDP^, 

l 

then  <Xq  is  exposed  (and  crossed  since  ug  <  usq). 

Consider  the  same  situation  as  that  discussed  in  the  paragraph 
above  but  with  ~^2r  <  ^b  30  ^at  c^e  SDP^  runs  UP  to  v  *  +°°.  Then 
we  may  encounter  into  the  situation  shown  in  Figure  C.8.  As  for  the 
case  when  y^  <  -k^,  we  that  ^  ao  *ies  between  the  upper  leg  of 

the  SDP  through  a  *  -(ir/2),  the  connection  line  and  the  SDP^  ,  then  ct^ 
is  exposed. 

Finally  we  consider  in  region  3.  We  saw  in  Section  A. 8  that  the 
slope  and  curvature  of  the  connection  line  was  given  by  Table  A. 2.  If 
we  consider  sets  of  these  curves  we  arrive  at  Figure  C.9.  We  observe, 
in  particular,  that  it  is  possible  for  to  lie  in  region  3  and  have 
its  connection  line  appear  in  region  1.  This  is  of  particular  signif¬ 
icance  since  the  SDP^  loci  must  lie  completely  in  region  3.  Therefore 
if  a  root  appears  in  region  1  between  the  connection  line  and  the  SDP 
through  -(tt/2),  it  is  exposed  and  crossed  for  ug  <  usq.  Otherwise,  if 
Vq>0,  root  is  never  exposed.  Now  assume  that  the  root  has  Vq  <  0. 

Assuming  that  u  <  u  ,  so  that  the  branch  cut  has  deformed  from  the 
s  sb 

initial  position  to  the  SDP^  and  exposed  a  portion  of  the  lower  sheet  of 
G  ,  then  any  root  in  the  region  bounded  by  the  lower  portion  of  the  SDP 
through  -(tt/2),  the  connection  between  -(tt/2)  and  ,  and  the  SDP^  will 
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Figure  C.7.  Possible  region  of  exposed  roots  of  G 
in  the  a  plane. 
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Figure  C.9.  Loci  of  connection  lines  from  -k^  to  -k^ 
in  the  a  plane. 
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be  exposed  and  crossed  (since  ug  <  usq)« 

We  can  summarize  the  results  of  this  analysis  in  a  manner  amenable 
to  implementation  on  a  computer.  We  accept  the  following  definitions: 


ug  is  the  given  saddle  point  value 

ugQ  is  the  saddle  point  value  for  which  the  SDP  runs  through  o the 
root  of  interest. 

(u^,v^):  ci^  *  u^+jv^,  ^b  *"S  C^e  *3ranc^  poi-nt  for  -k^ . 

(Uq.Vq):  ctQ  *  u0+jv0,  aQ  is  the  root  of  interest, 
u^  is  the  saddle  point  value  for  the  SDP  running  through  . 

is  the  value  of  the  phase  function  for  ct^  and  ug. 


voo„(u0); 


[if  ub>  -  Or/2)]  - 


SDP 


,  SDP  through  u  »  -  >  for  u^  <_  - 


conn 


,  value  on  connection  curve;  for  -  —  <u_< 


(-§)  to  ^ 


2<U0<Ub 


[if  Uj,  <  -  (TT/2)] 


VSDP’  SDP  throu8h  for  \  £  uq 


vSDp»  SDP  through  u  *  -  j;  for  -  -|<u0 


conn 


;  for  ^  <  uQ  <  -  - 


VSDP»  SDP  through  ab;  for  uQ  <_  \x^ 


on  SDPb;  for  vfe  <  vQ 


u^,  fixed;  for  vQ  < 


(vsdp,  SDP  through  V  ub  £  u0 


pk  <V  = 

D 


v  __  ,  on  SDP.;  for 
SDP^  o 


sb 


-«  ;  for  u_  < 


sb  2 


im 

con 


(O  -  -V  (-TT  -  u  ) 
'  0  con  □ 


vcOh<U0) 


for 


%i"o 


!St 


(uq) 
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We  enumerate  the  following  inequalities  which  will  be  use  in  the 
logical  constructions  of  the  flowchart  derived  from  the  earlier 
analysis : 


a. 

-!<‘s 

g- 

V0  >  Vcon 

m. 

_  im 

vn  >  V 

0  con 

b. 

u  <  u  n 
s  sO 

h. 

U0  <  USDP, 

D 

n. 

u  <  u 

S  Sb 

c. 

v,  >  0 

D  — 

i. 

Yu  >  k. 

'b  2r 

o. 

v„  >  V 

0  test 

i 

|  d. 

i 

TT 

2 

j. 

v  <  v  <  V 

con  o  SDP. 

D 

| 

!  e. 

vo-° 

k. 

V  <  V  <  V 

SDPV  0  con 

D 

f. 

uso  -  Usb 

1. 

U0  >  USDP, 

0 

We 

can  then  organize 

the 

logic  of  this  section 

into  the 

flowchart 

shown  in  Figure  C.10.  The  question  concerning  the  value  of  8^  to  be 
used  in  the  residue  contribution  of  the  root  of  G  does  not  arise 
because  the  form  of  the  contribution  (Section  B.5)  implicitly  includes 
the  correct  value  of  g^  due  to  (B.5.3). 


I 

i 

) 

I 

1 

i 

i 

i _ 

I  i&h.'  r  '  ;  .  . 


APPENDIX  D 


BEHAVIOR  OF  THE  GREEN'S  FUNCTION  IN  THE  SLAB  PROBLEM 


In  this  appendix  the  behavior  of  the  Green's  function  for  the  slab 
problem  is  considered  when  the  media  are  lossless,  have  the  same  permea¬ 
bility  U^and  the  slab  has  a  greater  permittivity  than  the  surrounding 
medium  (£  ^  >  ^he  Green's  function  is  given  in  (5.2.6)  and  we 

consider  here  the  functional  form,  F(k^),  given  in  (5.2.7).  In  normal¬ 
ized  coordinates,  we  have 


F(k;) 


r  ■je2t’ 

<B{-B’)e  * 

-  (S{+6')e 

2  j02tT 

-  (B|+8p2e 

■\ 


(D.l) 


F(kx) 


2S^t' 

—  — 

JB{t'  +  Bjt,cot(BJt') 

jC(B^t')2  +  (3£t')2]  +  2(6{t’)(B’t’)cot(8’t') 

» 

* 

(D.2) 


The  slab  is  taken  to  have  thickness  t'  (expressed  in  normalized  coordi¬ 
nates,  t  ■  t'  ‘kg).  Also  the  normalized  permeability  is  equal  to  1,  so 
k^  "  Only  the  roots  and  poles  of  the  expression  in  brackets  need  be 
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discussed  in  detail,  so  the  factor  3^  will  be  ignored.  We  follow  the 
analysis  of  Barone  (1956)  and  give  the  results  for  our  formulation  in 
the  present  notation. 

Consider  the  roots  of  the  numerator  of  (D.2), 

jB^t  +  S2t  cotS2t  ■  0  .  (D.3) 

Letting  0^t  *  w  *  y+jv  and  £2t  *  z  *  x+jy,  then  (D.3)  becomes 

uj  •  jz  cotz  .  (D.4) 

We  recognize  that 

z2-w2  -  (k2-k2)t2  -  (£2-£1)t2  -  a2  (D.5) 

where  a  ^  0.  Analyzing  (D.4)  and  (D.5)  together  allows  for  the  locus  of 

Che  roots  as  a  function  of  (a)  to  be  determined.  It  can  be  shown  that 

all  the  roots  reside  on  the  lower  sheet  except  for  a  finite  number  which 

must  lie  between  k,  and  k„  on  Che  real  k  axis  (and  their  images  which 
12  x 

lie  between  -k2  and  -k^  on  the  negative  real  axis).  Determining  the 
location  of  the  roots  requires  solving  the  equations, 

x2+  v2  -  a2  (D. 6) 

and 

v  *  ±  a  cos(x)  such  that  v  <  0,  simultaneously.  (D.7) 

This  can  be  thought  of  graphically  as  shown  in  Figure  D.l,  where  the 
locus  defined  by  (D.6)  is  shown  as  a  circle  of  radius  a  and  (D.7)  is 


Figure  D.l.  Graphic  determination  of  the  roots 
of  the  numerator  of  the  function  F 
In  the  slab  problem.  Roots  of 
x  ■  t  a  sin  x 
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shown  for  v  <  0.  Alternatively,  the  roots  of 
x  *  ±  a  sin(x) 


(D.8) 


can  be  determined  (clearly  |x|  £  a)  and  then  (D.7)  gives  v.  It  is 
evident  from  Figure  D.l  that  no  roots  occur  until  a£  (tt/2).  This 
being  the  case,  there  will  be  n  roots  on  the  interval  (k^ ,  ) ,  where  n 

is  the  integer  portion  of  (a/7r  +  1/2). 

The  roots  of  (D.3)  are  seen  to  correspond  to  the  surface  waves  of  a 
dielectric  slab  waveguide  of  thickness  t  with  one  surface  covered  by  a 
PEC  screen  (that  is,  the  so-called  short  circuit  bisection  problem).  It 
is  shown  in  the  asymptotic  analysis  for  the  transmitted  region  (y  <  -t) 
that  the  roots  give  rise  to  surface  waves  propagating  in  the  positive  x 
direction  from  the  edge  of  the  PEC  screen.  It  is  no  surprise  that 
these  surface  wave  roots  are  the  same  as  the  surface  wave  roots  of  the 
slab  structure  since  for  x  >  0,  the  configuration  is  coincident  with 
that  of  the  short  circuit  bisection  problem. 


Consider  now  the  roots  of  the  denominator  of  F(k  ), 

x 

-jB-t  jS,t 

(Bj^-Bj)  2e  -  (B^B^e  «  0  • 

After  some  manipulation  this  reduces  to  two  families  of  roots, 


(D.9) 


V  .  S2C  „ 

—  "  S  T  c°t 


a2t 


(D.lO.a) 


and 


B, t  B,t 
~Y~  — j-§-  tan 


S2t 


(D. 10. b) 
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Comparing  (D.lO.a)  with  (D.3),  it  is  clear  that  an  identical  analysis 
applies  using  (t/2)  rather  than  t.  In  particular,  from  (D.5)  we  see 
that  a  becomes  b  *  (a/2).  Therefore  we  have  no  roots  until  a  *  tt  and 
then  n  roots  such  that  n  is  the  integer  portion  of  (b/tr  +  1/2).  These 
roots  solve  (D.6)  and  (D.7)  or  (D.8)  with  the  appropriate  adjustment  in 
t  and  a.  From  the  earlier  discussion,  it  is  clear  that  these  poles 
correspond  to  the  short  circuit  bisection  problem  for  a  slab  of  thick¬ 
ness  (t/2).  Thus  these  waves  have  an  odd  functional  behavior  with 
respect  to  the  center  of  the  slab. 

Examining  (D.lO.b),  it  can  be  shown  that  using  (B^  t/2)  *  w  and 
(82C/2)*z,  the  equation  can  be  reduced  to 


and 


w  ■  -J z tan(z) 


(D.ll.a) 


2  2 
Z  -U) 


(VE1> 


(D.U.b) 


Hence  the  respective  real  and  imaginary  parts  of  the  roots  satisfy 

(D. 12 .a) 


and 


or 


x2  ♦  v2  -  b2 


v  ■  ±  b  sin(x)  , 


x  ■  +  b  cos(x) 


(D. 12. b) 


(D. 13) 


I 
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The  graphical  depiction  of  (D.12)  is  shown  in  Figure  D.2.  Observe  that 
this  pole  does  not  suffer  a  cutoff  for  b  <  ( "n-/ 2 )  unlike  the  short 

circuit  waves  discussed  previously.  These  waves  represent  the  even 
functional  waves  characterized  by  the  surface  waves  of  a  dielectric  slab 
of  thickness  t/2  with  one  surface  covered  by  a  perfectly  magnetically 
conducting  screen  (the  open  circuit  bisection  problem  for  a  slab  of 
thickness  t/2). 

By  combining  the  functions  shown  in  Figures  D.l  and  D.2,  it  can  be 
shown  that  as  a  (which  equals  /e2  -  t  )  increases  from  zero,  the 
following  sequence  of  roots  and  poles  appears  on  the  interval  (k. ,  k2) 
of  the  real  k^  axis  (and  by  symmetry,  on  the  interval  (-k^,  -k^)). 

First,  a  pole  emerges  from  k^  and  moves  toward  k^  as  a  increases  to  (tt/2). 
For  a  >  (tt/2),  a  root  then  emerges  and  follows  the  preceding  pole 

toward  k2.  Then  when  a  reaches  T,  a  second  pole  emerges  from  k-^  and 
moves  toward  k2«  This  alternation  between  poles  and  zeros  is  seen  to 
continue  as  a  increases.  This  progression  of  roots  is  shown  in  Figure 
5.2  which  describes  the  behavior  of  the  logarithm  of  the  Green's  func¬ 
tion  of  the  slab  problem. 


■vturtEs  w  - 


[*-(+  b) — ^(+  b)  — (-  b)  — (-  b)  — >] 

-  -  -  Roots  of  x  ■  ±  b  cos  x 

-  Roots  of  x  *  ±  b  sin  x 

Figure  D.2.  Graphic  determination  of  the  roots 
of  the  denominator  of  the  function 
F  in  the  slab  problem. 
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